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Chapter 1 



Preface 



The aim of this work is to present a role of the reggeized gluon states in perturbative Quantum 
Chromodynamics. The reggeized gluon states, also called Reggeons, appear in the Regge limit 
where the square of the total energy s is large while the transfer of four-momentum t is low and 
fixed. In this limit the leading contribution to the scattering amplitude of hadrons is dominated 
by the exchange of intermediate particles, Reggeons, which are the compound states of gluons 

mca nam 

Even for the Regge limit in the generalized leading logarithm approximation [3 HI El this 
problem is technically very complicated due to the non-abelian structure of QCD. Therefore, 
in order to simplify colour factors the 't Hooft's multi-colour limit [HI HH E21 is performed, 
in which a number of colours N c — > oo. This causes that the Reggeon wave-functions become 
the eigenstates of a Hamiltonian which is equivalent to the Hamiltonian for the non-compact 
Heisenberg SL(2, C) spin magnet. Because of this symmetry the multi-Reggeon system with N 
particles is completely solvable [IHllIlj. Thus, it possesses a complete set of integrals of motion, 
conformal charges (q 2 , q 2 , (fe, q 3 , ■ ■ ■ , (?jv>?jv)- The eigenvalues of the SX(2,C) Hamiltonian are 
also called the energies of the Reggeons. The Schrodinger equation for the lowest non-trivial 
case, i.e. for N = 2 Reggeons, was formulated and solved by Balitsky, Fadin, Kureav and 
Lipatov [H3 EH Sj. They calculated the energy of the Pomeron state with N = 2 Reggeons. 
An integral equation for three and more Reggeons was formulated in Refs. [HI Ej in 1980. 
However, it took almost twenty years to obtain the solution for iV = 3, which corresponds to 
the QCD odderon [THJ [Tj3 [20] • Finally, the solutions for higher N — 4, ... ,8 recently in series 
of papers [21112211231 written in collaboration with S.E. Derkachov, G.P. Korchemsky and A.N. 
Manashov. 

Description of scattering amplitudes in terms of the reggeized gluon states is most frequently 
used in two different kinematical regions. The first region appears in calculation of the elastic 
scattering amplitude of two heavy hadrons whose masses are comparable. Here the amplitude 
is equal to a sum over the Regge poles. They give behaviour in the Regge limit, s — > oo and 
t = const, like s a ^ where a(0) is called the intercept and its value is close to 1. On the other 
hand, the intercept is related to the minimum of the Reggeon energy defined by the SL(2, C) 
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Hamiltonian. Thus, evaluating the spectrum of the SL(2, C) XXX Heisenberg model we can 
calculate the behaviour of the hadron scattering amplitudes. 

The second kinematical region is the region of deep inelastic scattering (DIS) where a 
virtual photon 7*(Q 2 ) is scattered off a (polarized) hadron with the mass M 2 = p 2 , such 
that Aq CD <C M 2 <C Q 2 with x = Q 2 /2(p ■ q) and q 2 = —Q 2 . In this case the moments of 
the structure function F 2 (x, Q 2 ) can be expanded (OPE) in a power series in 1/Q 2 where the 
exponents of this expansion are related to the twist of Reggeons and anomalous dimensions 
of QCD. Analysing the spectrum of the SL(2, C) Hamiltonian one can find a twist related to 
a given Reggeon wave-function and calculate the corresponding anomalous dimension. This 
was done for the leading twist n = 2 in Ref. |2I| and for other twists of two-Reggeon states in 
Ref. [2H]. The anomalous dimension and twists for more than two Reggeons were discussed in 
Ref. |25| ■ 

In this work we present the methods of constructing the Reggeon eigenstates. Moreover, 
we calculate the reach spectrum of the the Reggeon energy and the conformal charges {qk,Qk\- 
Finally, we calculate anomalous dimensions of QCD and twists in operator product expansion 
(OPE) which are also provided by the Reggeon states. In the first three Chapters we discuss 
the current state of knowledge. The main results of the present work are collected in Chapters 
5-7 and in Chapter 9. For completeness in Chapter 8 we added description of another method 
presented in Refs. [271128] . 

Thus, in Chapter 2 we explain when the reggeization of the gluon appears [21] • Next, we 
perform the multi-colour limit [101 HU H2] , discuss the properties of the SL(2, C) symmetry and 
construct invariants of this symmetry, the conformal charges. In next Chapter we introduce a 
Baxter Q— operator method j^Dl with Baxter equations which allows us to solve the Reggeon 
system, completely. Chapter 4 contains the quasi-classical solution for the Baxter equations |21j. 
It gives us the WKB approximation for the spectra of the Reggeon energy and the conformal 
charges. It also explains the structure of these spectra. 

In Chapter 5 we show construction of the Reggeon eigenfunctions which consists in solving 
the the eigenequations of the conformal charges. We systematize the knowledge about the 
ansatzes for eigenstates of {qk, q k } with an arbitrary number of Reggeons, N, extend the calcu- 
lations to an arbitrary complex spins s and derive differential eigenequations for the conformal 
charges with N = 3,4. Moreover, we show solutions to the differential eigenequations with 
N = 3 and s = and resum obtained series solutions for the q 3 = case. In next Chap- 
ter we present an exact solution to the Baxter equations. It consists in rewriting the Baxter 
equation into the differential equation, which may be solved by a series method, and finding 
the quantization conditions for {qk,%} which comes from single- valuedness of the Reggeon 
wave-function and analytical properties of the Baxter functions. The numerical results of this 
method for N — 2, . . . , 8 Reggeons are shown in Chapter 7 [221 123) • I n the A" = 3 case we 
present quantized values of q 3 for different Lorentz spins = 0, . . . , 3 as well as corrections 
to the WKB approximation. For A" = 4 we show the resemblant and winding structure of the 



14 



q4 and spectrum and also corrections to the WKB approximation when g 3 = 0. At the end 
we discuss the ground states with N = 2, . . . , 8 Reggeons and describe their properties. In 
Chapter 8 we review another approach to the Baxter equation which was presented in Refs. 
(23 [2B]. We discuss the advantages and disadvantages of this method and compare it with the 
results presented in previous Chapters. Finally, in Chapter 9 we concentrate on deep inelas- 
tic scattering processes. We calculate anomalous dimensions of QCD and twists coming from 
Reggeized gluon states [26] . One can obtain them by performing an analytical continuation of 
the Reggeon energy into the complex space of the scaling dimension Uh and analysing the pole 
structure of the energy. At the end we make final conclusions. 
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Chapter 2 



Introduction 



2.1 Hamiltonian for the iV-Reggeon states 
2.1.1 The gluon reggeization 

The Regge limit can be conveniently illustrated by considering the elastic scattering amplitude 
of two hadrons 




P2 

in the kinematical region: 



with s = (pi +p 2 ) , and t = (p 1 - pi) 2 , (2.1) 



s — > oo, t — const . (2.2) 

In this case the largest contribution to the scattering amplitude is provided by a multi-gluon 
exchange. It was shown in the papers |31 IU El El UH HH1 EH] that in the limit (I2.2jl so called 
gluon reggeization occurs. It means that the contribution to the scattering amplitude may be 
written as a sum of ladder diagrams, where compound states of the reggeized gluons, called 
Reggeons, are exchanged in the t— channel and interact with each other. 




a) Born approx. b) composite virtual and Born diagrams c) gluon emission 



Figure 2.1: Some diagrams of the order g A and g 6 which give contribution to the total cross 
section cr 2 ^2- 

In order to illustrate what the gluon reggeization is, we show in Figure 1 some diagrams of 
the order g 4 and g 6 which give contribution to the total cross section cr 2 ^2 for two scattered 
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gluons. The Born approximation diagram, Figure l.a, appears in the lowest order of the 
perturbative calculation. Next diagrams, Figure l.b and l.c, are of the order g 6 . They include 
virtual contributions and real gluon emissions. 

The non-triviality of the virtual diagrams due to the colour factor causes modification of 
the gluon propagator. This phenomenon is called gluon reggeization. A leading contribution in 
Ins coming from the diagrams of the l.b type, i.e. of the order g 6 , leads effectively to exchange 
of the gluon propagator 

1 u(t)his , f d 2 l 1 , , 

where u(t) = 3a s t \ ;9/ ^_ (2.3) 



t t w J (27r) 2 /2( v ^t -if 

Note that u(t) in 1)2.3)) is infrared divergent. Here we perform regularization adding a small 



gluon mass (HES). Summation of the leading terms obtained from the virtual diagrams gives 

(2.4) 

Similarly the real gluon diagrams from Figure l.c lead to the modification of vertices re- 
sponsible for the gluon emission. The j-th vertex give contribution 

kj Bj dj 

sPit^tenOi+uq)), . T (2.5) 

where F(qj + \,qj) is a function defined for example in Ref. [29j, where T is the gluon colour 
matrix. These vertices with the emitted gluons serve to build the ladder rungs. 

In the ladder diagrams there are also Reggeon couplings with two scattering particles. For 
gluons they look like 

in+l 





2Pi^a 1 a 1 T- 1 ^KyXjT (2.6) 

where 5\\> means conservation of helicity. 

One can see that the absorptive part of the gluon scattering amplitude for two gluons may 
be written in the form of the ladder diagrams summed over the intermediate n-gluon states 

lmA(s,t) = ^X] / dT n+2A* 2 ^ n+2 (p' 1 ,P2)A 2 ^ n+2 (pi,P2) , (2.7) 

n 

where <ir n+2 is an element of the phase space for n + 2 intermediate particles, whereas A 2 ^2+n 
is an amplitude represented by an effective diagram in Figure 2. 

In this diagram, momenta Pi and k , k n+ i and also helicities Aj and A^, for % = 1,2 char- 
acterize two initial and two final states. The remaining intermediate gluons with momenta ki 
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Pi Ai ai 



P2 A 2 a 2 k n+ i A' 2 a 2 

Figure 2.2: Diagram symbolizing the scattering amplitude A 2 ^ n+ 2 

numbered from 1 to n are characterized by polarization vectors ej, where j = 1, ...,n. Simi- 
larly for the colour indices, the two initial gluons and the two final ones have the colour indices 
ai, a\ where % = 1,2, whereas the indices of the produced gluons are denoted as dj where 
j = 1, . . . ,n. The Reggeons transfer squares of 4-momenta tj = q 2 and possess colour indices 
ij, where j = 1, . . . , n + 1. 

The total contribution of the diagram from the Figure 2 to the scattering amplitude has 
the form 

A 2 ^ n+ 2 = 2sg n+1 5 XiX[ 5 X2 y 2 

E ^a^TfX^F^q^^Tf^F^q,)) . . . -f^-T^ 1 , (2.8) 



x 



n+1 



where s it j = (ki + kj) 2 . 

One can decompose the momenta ki in a basis of the scattered particles momenta 

h = otiVx + PiP2 + k ix . (2.9) 

In the limit (|2.2J) the largest contribution comes from the kinematical region where Oj-i 3> 
(3j-\ <C /3j and small gj ~ g| . Taking square of the diagram in Figure 2 one can perform 
integration over longitudinal degrees of freedom. All non-trivial dynamics is then concealed in 
the 2-dimensional subspace of transverse variables. 

In Eq. p.7jl one has to sum over polarizations and colours of intermediate particles. The 
sum over polarizations gives 

^^(^.ft)^^^) = -F,{q l+l , qi )g^FM + i^[) = 2K q ( gi ,q i+1 ), (2.10) 

i 

where 



K q (quq i+ i) = q r- 2 (2.11) 

K i± 



with t = q 2 . In order to perform the sum over colours it is convenient to introduce projectors 
P( R ) onto irreducible representation (R) composed of two colour octets in the t— channel: 

(8) ® (8) = (1) © (8) 5 © (8) A © (10) © (10) © (27). (2.12) 
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The scattering interaction of two reggeized gluons may be then written as 

A^oiCM) = T, P £ia ia[ A {R) (s,t), (2.13) 

(R) 

where A^ R \s,t) is the amplitude with specified colour state in the t channel. For the i-th 
intermediate gluon we obtain a colour part in the form 

di R 

where \r is the colour factor. 

Summarizing, the absorptive part of the scattering amplitude for two gluons with conser- 
vation of helicities is given as jHS [33] 



n n n 

x n g.+o^n - ^ n i^d • ( 2 - 15 ) 

i=l i=0 i=l 

The reggeization may only occur in the non-abelian theories. In contrast to the quantum 
chromodynamics, the reggeization of particles does not appear in such theories as quantum 
electrodynamics as well as scalar field theories. These theories do possess neither a colour 
structure nor three-gauge-field vertices, which are responsible for gluon emission, and as a 
result in these theories u(t) = 0. 



2.1.2 Integral equation for the scattering amplitude 

Performing the Mellin transformation on (j2. 15^ . the sum over n and introducing new 2- 
dimensional variables in the transverse subspace: 

h = ?i, h = -q[ = q-qi, 

l[ = q n +i, = -q'n+i = ? ~ 1n+l, (2-16) 
t { = -I 2 < 0, t = -q 2 < for i = 1, 2, 

we obtain 

where the function m q '(li) describes the gluon ladder with two external reggeized gluons and 
with the fixed total angular momentum j and the transferred momentum q. This function 
satisfies the integral BFKL equation [TH1 [TJl 0] 

{j-l-uj{t 1 )-u{t 2 ))^\h) = l + 2a s \ R J ^V q (hJ[)C(l[), (2.18) 

20 



where a s = g 2 /4:7i, and the interaction potential between two Reggeons is 

Equation ()2.18|) may be illustrated using the diagrammatic form shown in Figure 3. Notice 
that Eq. ()2.18|) does not contain the propagators of out-going gluons related to the momenta 
h and l 2 . 

(i-i) - - + 




Figure 2.3: Illustration of the BFKL equation by means of Feynman diagrams. 

Equation (I2.18j) contains infrared singularities. The Reggeons in the t— channel may form 
different colour representations ()2.12|) . For the singlet A^ = i = 3. In this case, for q ^ all 
infrared singularities cancel. For q = 0, the integrand kernel in Eq. ()2.18|) contains singularities 
at l[ = 1' 2 = 0. Equation (|2.18j) describes the process with a gluon of momentum p 2 being a 
target. In reality we are not interested in the scattering on a gluon, but on a physical colourless 
particle. In this case thanks to inhomogeneity in Equation ()2.18|) . which describes the Reggeon 
interaction with a target, all infrared singularities disappear. 

Equation (I2.18|) has the form of the inhomogeneous integral Schrodinger equation 

(H 2 -E)®(j) = f, where H 2 = w(ti) + u(t 2 ) + 2a s X 1 V 12 , 8 = j - 1, (2.20) 

whereas <S> means convolution of the Hamiltonian with the wave-function (p. 

Equation (|2.20p describes exchange of two reggeized gluons. It is possible to generalize it 
to more gluons. The kinetic energy of the i-th Reggeon is equal to ui{ti). The potential of the 
interaction between two Reggeons has the form 6a s Vi2, where V\ 2 (I2.19J) is an integral operator. 
For N Reggeons we have 

(H N - S) <8> = /, (2.21) 
where the Hamiltonian, which was firstly derived in the works [3 HI EH takes a form 

N N 

H N =y]o;fti) - 2a s ^ T i T kV ik 

(2 _ 22) 

= -i £ T^M*) +6«A) = -~ E TiT k H ik . 

i<k=2 i<k=2 

The interaction potential V ik between i-th and k-th reggeized gluons is defined by Eq. (|2.19|) . 
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2.1.3 Impact parameter representation 

The Hamiltonian ()2.22|) was defined in two-dimensional space of transverse Reggeon momenta. 
It appears that it is more convenient to perform calculations using variables in the impact 
parameter space using the transverse spatial coordinates of Reggeons (xj,yj). We pass to this 
representation performing Fourier transform of Eq. (I2.22J) . After introducing holomorphic and 
anti-holomorphic complex coordinates for the j-th reggeized gluon 

Zj = Xj + ii/j, Zj = Xj — ii/j (2.23) 



the Hamiltonian for a pair of Reggeons H ik defined in Eq. (|2.22j) separates into two independent 
parts: 

H lk = -^(H{z h z k )+H{z h z k )). (2.24) 

Z7T 

Now, one can see that (|2.21j) was reduced to two one-dimensional Schrodinger equations. 
The holomorphic part and analogously the anti-holomorphic one has the form 

H( Zj , z k ) = -Pr 1 \n.(z jk )Pj - P, 1 Hz jk )P k - HPiPk) + 2^(1), (2.25) 

where Pj = —id Zj , and Zj k = Zj — z k . 

The Hamiltonian (12.251) is SL(2, C) conformal invariant. It means that is invariant under 
translation, scaling and inversion operation [2H]. An eigenstate of the full Hamiltonian (|2.22|) 
may be written as the bilinear form <P = ^ x \P. In this work we demand that this function is 
single- valued and normalized. 

The colour factor T{F k in the Hamiltonian (I2.22J) may be calculated exactly for = 2 (i.e. 
the leading order of Pomeron exchange) and for = 3 (i.e. the leading order of odderon and 
the next to leading order of Pomeron exchange). For > 3 we use an approximation N c — > oo. 
Finally we obtain the Hamiltonian of the nearest neighbour interaction 

— N 

H N = H N +t N = ^Yl ^+0 + H(zk, Zh+i)] (2-26) 

fc=i 

with a s = a s N c /ii where N c is the number of colours, and z\ and z^+x correspond to the same 
Reggeon. 

The total cross section for the exchange of two reggeized gluons in the singlet representation 
(Pomeron) is proportional to s w W. This behaviour is in contradiction with unitarity of the S 
matrix. There is a conjecture [3 02] that in order to conserve unitarity one has to sum up 
diagrams for all N- Reggeon exchanges. Additionally, as discussed in the papers [SSHMlElj the 
multiple Pomeron exchanges, so called the fan diagrams, are more likely to unitarize the total 
hadronic cross section [33|. In the latter diagrams the number of the Reggeons varies in the 
t— channel. Those diagrams are beyond the scope of this work. 
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2.2 System with SL(2, C) symmetry 
2.2.1 Symmetry SL(2, C) 

The Hamiltonian (|2.26|) is invariant under the coordinate transformation of the SL(2, C) group 

z , = az k + b _, = az k + b 

Zk czk + d' Zk cz k + d [ ' 1 

with k = 1, . . . , N and ad — be = ad — be = 1. 

Now one may associate with each particle generators of this transformation [38J. This 
generators are a pair of mutually commuting holomorphic and anti-holomorphic spin operators, 

ie aPl 5 kn S\ k) and 



similarly for S a . The generators act on the quantum space of the k-th particle, V^ Sk,Sk ^ as 
differential operators 

#0 = z kd Zk + s k , = -d Zk , = z 2 k d Zh + 2s k z k , ^ 

Sq = z k c\ k + s k , = -c\ k , = ~z~ldz k + 2s k z k , 

where S± = iiS^ while the complex parameters, s k and s^, are called the complex spins. 
Thus, the Casimir operator reads 

2 

J2(Sf f = {S {k) f + (S {k) S {k) + S {k) S {k) )/2 = s k (s k - 1) (2.29) 

3=0 

f k \ 

and similarly for the anti-holomorphic operator (S ) . The eigenstates of 5X(2,C) invariant 
system transform as ERT] 

*{z k ,z k ) -> tt'(z fc ,z fc ) = (c2 fc + rf)- 2 ^(cz fc + d)- 25 ^(4,4). (2.30) 

Due to the invariance (|2.27jl of the system we can rewrite the Hamiltonian (|2.26j) 

Hn = H N + H N , [H N , H N ] = (2-31) 

in terms of the conformal spins ()2.28|) 

N N 

H N = J2 H (Jk,k+i), Hn = J2 H ( J w)> H(J) = i/>(1- J) 2i>(l) (2.32) 
fc=i k=i 

with = dha.T(x)/dx being the Euler function and Jn,n+i — Jn,i- Here operators, J k)k +i 
and Jfc,fc+i, are defined through the Casimir operators for the sum of the spins of the neigh- 
bouring Reggeons 

J k>k+ i(J k , k+ i ~ 1) = (S W + 5 (fc+1) ) 2 (2.33) 
with Sl iV+1 ' > = «9q 1 , and J k ,k+i is defined similarly. 
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In statistical physics (j2.32j) is called the Hamiltonian of the non-compact SL(2, C) XXX 
Heisenberg spin magnets. It describes the nearest neighbour interaction between N non- 
compact SL{2, C) spins attached to the particles with periodic boundary conditions. 

For the homogeneous spin chain we have to take s k = s and St = s. In QCD values of 
(s, s) depend on a chosen scalar product in the space of the wave-functions (|2.30j) and they are 
usually equal to (0, 1) or (0,0) jEEHZZj. 

2.2.2 Scalar product 

In order to find the high energy behaviour of the scattering amplitude we have to solve the 
Schrodinger equation 

Wft = °'^ 1, *(*i> ■ ■ ■ > ?n) = En*&> *2, • • • , z N ) (2.34) 

with the eigenstate ty(zi,z 2 , . . . ,z N ) being single-valued function on the plane z — (z,z), nor- 
malizable with respect to the SL(2, C) invariant scalar product 

H^H 2 = = J d 2 z 1 d 2 z 2 ...d 2 z N \V(z 1 ,z 2 ,...,z N )\ 2 , (2.35) 

where d 2 Zi = dxidy-i = dzid~Zi/2 with Zi = z*. 

It is possible to use other scalar products [2S]. Let us consider the amplitude for the 
scattering of two colourless objects A and B. In the Regge limit, the contribution to the 
scattering amplitude from N— gluon exchange in the t— channel takes the form 

A(s,t) = is^2(ia s ) N A N (s,t). (2.36) 

N 

Using the SL(2, C) scalar product ()2.35|) we have 

A N (s,t) = s J d 2 ^ e^■P(<l A (f )|e-^ y ^/ 4 (9^..4) _1 |<l B (0)), (2.37) 

where = d/dz k , the rapidity Y = Ins and . . .dj^j are gluon propagators omitted in 
(j2.18j) . Here the Hamiltonian Hn is related to ()2.26|) as Hn = —Q-s "Hn/^, so it is given by the 
sum of N BFKL kernels corresponding to nearest neighbour interaction between N reggeized 
gluons. The wave-functions \&a(b)(z6)) = $A{B){zi ~ %o) describe the coupling of N— gluons 
to the scattered particles. The z — integration fixes the momentum transfer t = —ft 2 whereas 
the operators l/<9f stand for two-dimensional transverse propagators. 
Defining the functions <3>(z) as 

$(z) = (-i) N d zl d Z2 . . . d ZN §{z) (2.38) 

the scalar product in the amplitude (|2.37|) can be rewritten as 

($ A (f )|e-^ y ^ /4 (St . . .^) _1 |$ B (0)> = (<S> A (Zo)\e-« sYn » =0 '" 1)/4 |$b(0)> • (2-39) 
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The Hamiltonians, 1~Ln and Hn are related to each other as 

n (s=o,s=i) = ^ ^ ^ ^ g N y\ (2 40) 

From the point of view of the SL(2, C) spin chain 

H N = H N =0 '° =0) . (2.41) 

Indeed, the transformation S a — ► (<9i . . . 9a?) 5q, (9i . . . 9a?) 1 maps the SL(2, C) generators of 
the spin s = into those with the spin s = 1. One concludes from (j2.41j) that the Hamiltonian 
^(s-o,s-i) j g advantageous with respect to Ti n as it has the quantum numbers of the principal 
series representation of the SL(2, C) group. 

However, one can also use H N ~ 0,S ~°" > [22j or even H%T 1 ' a ~ 1) (2H]. Then the factor (di . . . d N ) Tl 
or (pi . . . <9at) t1 may be included into the scalar product, i.e. for (s = 1, s — 1) we have 

||^|| 2 = J dh 1 d 2 z 2 ...d 2 z N \(d 1 ...d N y 1 ^(z 1 ,z 2 ,...,z N )\ 2 . (2.42) 

Here the scalar product is no longer in the principal series representation of the SL(2, C) group. 
All these Hamiltonians with the corresponding scalar products are equivalent up to the zero 
modes of the (di . . . <9tv) t1 and (pi . . . d^) Tl operators. 

2.2.3 Conformal charges qk and the conformal spins 

The Hamiltonian (|2.31|) possesses a complete set of the integrals of motion {p, q k } where q k = 
{Qk: Qk} with k = 2, . . . N are called conformal charges while p = {p,p} is the total momentum 
of the system. In order to construct them we introduce the Lax operators [1H H21 EC3 S3 in 
holomorphic and anti-holomorphic sectors: 



L k {u) = u + i((T-S {k) ) 



u + iS^ iS<L h) 



L k {u) = u + i{a-S )= {k) _ {k) ) (2.43) 

V ib + u — tb 

with u and u being arbitrary complex parameters called the spectral parameters and a a being 
Pauli matrices. 

To identify the total set of the integrals of motion of the model, one constructs the auxiliary 
holomorphic monodromy matrix 

T N (u) = Li(u)L 2 (u)...L N (u) (2.44) 

and similarly for the anti-holomorphic monodromy operator Tn(u). Taking the trace of the 
monodromy matrix we define the auxiliary transfer matrix (spectral invariants) 

i N (u) = tiT N (u) = 2u N + q 2 u N - 2 + ... + q N (2.45) 
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and similarly for We see from 1)2.451) the advantage of using the transfer matrix tjv(u): 

that is a polynomial in u with coefficients given in terms of conformal charges % and q k , which 
are expressed as linear combinations of the product of k spin operators: 



TV 2 



& = - 2 E E K* 1 ^) (2-46) 
«2>«i=i ii=o 

N N 2 

* = - E E E ( W) (4"^) (2-47) 

i2>«l=l «4>i3 = l ilJ2=0 
^ AT AT AT 2 

<?6 = — g ^ E E £ »1*2»3*4*5»6 (^ji ^jl 



i 2 >ii=l i 4 >i 3 =l i 6 >i 5 =l ju' 2 ,j 3 =0 



x 4 3) <^ 4) ^; 5) 4 6) > ( 2 - 48 ) 



where Ei x i 2 ...i k is completely anti-symmetric tensor and £i 1 i 2 ...i fc = 1 for Zi < z 2 < • • • < 4- So for 
even /c we have a formula for conformal charges 



Qk 



(k/2)\ 



N 

E 

i 2 > i\ = 1 

«4 > «3 = 1 



E 



fc/2 



'U«2---*fc 



n a 



jl J2vJfe/2=0 



n=l 



jt(l2n-l) C<(*2n) 



(2.49) 



^ ^n— 1 



For odd fc's we have 



N 

<73 = 2 ^11213 3 ^ 

H,«2,«3 = l 

AT 2 

= 2 ^ ^ e iii2j3 e iu2j3 {^ji^j^^hj ' (2.50) 

«3>«2>«1=1 jl J2 J3=0 

AT AT 2 

& = -2 E E E^^^^f^fOHr^f) ' ^ 2 - 51 ) 

i3,«2,«l=l 15>M = 1 J4=0 

N N N 2 

A V V V V * f q(<0 ofeO qfaA 

V7 — / v / y / v / y fc U«2i3«4«5i6*7 I °0 °1 °2 J 

«3,«2,«1 = 1 «5>«4 = 1 i7>«6 = l j4,J5=0 

x (^jf) (4 6)5 1i 7) ) ( 2 -52) 
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and the general expression for an odd number of the conformal spins is 

N N 2 



q k = 2 y~ 1 > sr sr sr .... fcfrOcteOcfco 



EST e ( c(*a) eft 

( t hi2-i k \ o Q J 1 2 

V 2 J' i 3 ,i2,U=l . . Jlj2,->i(fe-3)/2=0 

«5 > «4 = 1 



Z 7 > 2 6 = 1 

i k > 4-i = 1 

(fc-3)/2 

x n (^n 2 " +2) ^i 2n+3) ) • ( 2 - 53 ) 

n=l 

In the above formulae we have two basic blocks [Sj^S^j and ^Sq^S^S^J whose 
products are summed with antisymmetric tensor 6^...^. 

2.2.4 Two-dimensional Lorentz spin and the scaling dimension 

The Hamiltonian (|2.32|l is a function of two-particle Casimir operators [HE!, an d> therefore, it 
commutes with the operators of the total spin S a = and S a = ^2 k S^\ acting on the 

quantum space of the system V N = V^ Sl ' Sl ^ ® ^( S2 > S2 ) <g> . . . <g> ]/( s iv,sjv)_ ^his implies that the 
eigenstates can be classified according to the irreducible representations of the SL{2, C) group, 
V^ h ' h \ parameterized by spins (h,h) [38J. 

The Hamiltonian depends on differences of particle coordinates so eigenfunctions can be 
written as 

tf^zi, z 2 ,...,z N ) = J d 2 z e iS °-^{z x - z , z 2 -z ,...,z N - z ) . (2.54) 

The eigenstates ^(5*1, z 2 , ■ ■ ■ , z N ) belonging to V^ h ' h > are labelled by the centre-of-mass coordi- 
nate z and can be chosen to have the following the SL(2, C) transformation properties 

" 4'}) = (cz + d) 2h (cz + rf) 27r (j[( cz k + d) 2s *(cz k + d) 2 ^ *{{z k - z }) (2.55) 

with Zq and z transforming in the same way as z k and z k , 1)2. 27J) . As a consequence, they 
diagonalize the Casimir operators: 

(S 2 - h{h - l))*0?i, z 2 ,...,z N ) = (2.56) 

corresponding to the total spin of the system, 

N 2 N 

s2 = £ £ ] = - £ ^ - !) • ( 2 - 57 ) 

12,11=1 j=0 fc=l 

The complex parameters (s k ,s k ) and (/i, h) parameterize the irreducible representations of 
the SL(2, C) group. For principal series representation they satisfy the conditions 

s k -s k = n Sk , s k + (s k )* = l (2.58) 
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and have the following form 

1 + n s , _ 1 — n s , . . 

s k = — y~ + lv *k> Sk = — Y^ + W s k ( 2 - 59 ) 

with i/ Sfc being real and n Sk being integer or half-integer. The spins (h,h) are given by similar 
expressions with n Sfc and v Sk replaced by rih and z^, respectively 

h = — Viv h , h = — Y%v h . (2.60) 

The parameter n Sk has the meaning of the two-dimensional Lorentz spin of the particle, whereas 
v Sk defines its scaling dimension. To see this one can perform a 27r-rotation of the particle on 
the plane, and find from eigenstates transformations (I2.55jl that the wave-function acquires a 
phase. Indeed 

z k -> z k exp(27ri) and z k -> z k exp(-2ni) gives ^(z k ,z k ) (-l) 2ns k^(z k ,z k ). 

(2.61) 

For half-integer n Sk it changes the sign and the corresponding representation is spinorial. Sim- 
ilarly, to define scaling dimension, s + s = 1 + 2iu Sk one performs the transformation 

z —> \z and z -> Xz giving ty(z k , z k ) -» \ l+2iUs k^(z k , z k ) . (2.62) 

Because the scale product for the wave-functions is invariant under SL(2, C) transformations, 
(I2.27L the parameter v Sk is real. 

We notice that the holomorphic and anti-holomorphic spin generators as well as Casimir 
operators ()2.33|) are conjugated to each other with respect to the scalar product l|2.35j) : 

[SL k) ] ] = -S ( " } , [J k ]1 = l-J k . (2.63) 

Moreover, because of the transformation law (|2.63jh h* = 1 — h 1 . This implies that = H n 
and, as a consequence, the Hamiltonian is hermitian on the space of the functions endowed 
with the SL(2, C) scalar product, W N = TLn- 

2.3 Conformal charges qk as a differential operators 

We noticed in the previous paragraphs that the conformal charge operators q k are given by 
invariant sum of linear combinations of the product of k spin operators. They can be rewritten 
as fc-th order differential operators acting on (anti) holomorphic coordinates (z,z). 
Two particle spin square can be written as 

Y,4 l] 4 2) = -^(^-^)X^ 2 +(^-^)(^^n+^^ 2 ) + s ^. (2.64) 
ii=o 

1 * - denotes complex conjugation 
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For homogeneous spins s = s\ = s 2 = ■ ■ ■ = sn } what is also the QCD case, we have 

N / 2 \ N 

& = - 2 E E^: 2) = E (fe 1 ) 2(l - s) ^ 2 ^ i (^ 1 ) 2s +2 S ( S -i)), 

«2>U = 1 \jl=0 / t2>U=l 

JV AT 

11,12,13=1 i3>«2>U=l 

+ SZi x i 2 (Zi 2 i 3 — z i 3 i-i)d Zi2 d Zii + SZj 2 j 3 (^igjj — z iii 2 )dzi 3 dzi 2 

~^~ sz i3il (^13H Z i\i-i)^Zi,J^z i2 2s Z^^d^ 2s Z{ 2 i 3 d Zi ^ 2s ZjgjjC^) , (2.65) 

where 2y = — Zj. Similar relations hold for the anti-holomorphic sector. 

In that way one can also construct operators for the higher conformal charges. They have 
a particularly simple form for the SL(2, C) spins s = 

q k = i k Yl z hh ■ ■ ■ z 3 h -i,h z h,h d z h ■ ■ ■ d Zjk i d Zjk (2.66) 

l<jl<32<-<jk<N 

as well as for s = 1 

Vk = i k d *h ■ ■ ■ dz h-i dz >u z hn ■ ■ ■ z h-i,h z h,h • ( 2 - 67 ) 

i<ii<j2<...<j fc <Af 

2.4 Other symmetries 

The states (I2.54B have additional symmetries [38J: 

W q>7! (zi, z 2 ,..., z N ) = V qrq {z 2 , z 3 ,..., zi) = e i6 »^ q , q (zi, z 2 ,..., Z N ) , 
M^ r± (fi, Z2,..., z N ) = ■V ± (z N ,z N „ 1 ,...,z 1 ) = ±'V ± (z h z 2 ,...,z N ) 



(2.68) 



so called cyclic and mirror permutation where the conformal charges are denoted by q = 
(q 2 , q 3 , ... , q n ) and q = (q 2 ,q 3 , . . . ,q n ). They generators P and M, respectively, commute with 
the Hamiltonian 7i but they do not commute with each other. They satisfy relations 

P W = M 2 = 1, pt = p-i = pJV-i ) M f = M, PM = MP _1 = MP JV_1 . (2.69) 

The phase 6^{q) which is connected with eigenvalues of P is called the quasimomentum. It 
takes the following values 

9 N (q, q) = 2tt A fo r k = 0, 1, . . . , JV - 1 . (2.70) 

The eigenstates of the conformal charges qu diagonalize TC and P. 

The transfer matrices (|2.45j) are invariant under the cyclic permutations t N (u)F = t N (u) 
whereas they transform under the mirror transformation as 

Mfjv(u)M= (-l) N i N (-u). (2.71) 
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Substituting (|2.45j) into (|2.71jl one derives a transformation law of the conformal charges q k 

^q k F = q k Mq k M=(-l) k q k (2.72) 



and similarly for the anti-holomorphic charges. Since the Hamiltonian Ij2.32j) is invariant under 
the mirror permutation, it has to satisfy 

H(q k j k ) =MH(q k j k )M = H(Mq k M,Mq k M) = H((-l) k q k , (-1)%) . (2.73) 

This implies that the eigenstates of the Hamiltonian (I2.32J) corresponding to two different sets 
of the quantum number {q k ,q k } and {(—l) k q k , {—l) k q k } have the same energy 

E N (q k ,q k ) = E N ({-l) k q k , (-l) k q k ) . (2.74) 

Similarly one can derive relation for quasimomentum 

d N (q k ,q k ) = -9 N ((-l) k q k , {-!)%) . (2.75) 

The cyclic and mirror permutation symmetries come from the Bose symmetry and they 
appear after performing the multi-colour limit [10] • Physical states should possess both sym- 
metries. 

2.5 Pomeron and odderon 

In the previous Sections we introduced the wave-functions of the Reggeized gluon states. Now, 
we answer how one can distinguish Pomeron states from odderon states |45| . Since the Pomeron 
states have conjugation charge parity C = 1 and the odderon states have C = —1, we have 
to check the C— parity of studied states. Let us consider for simplicity purely perturbative 
gluon fields A M (x) = A^(x)t a , where t a are the generators of the gauge group SU(N C ) and x are 
coordinates of the gluon. Thus, under a charge conjugation transformation this field transforms 
as 

A,(x) - -Aj(rr) , (2.76) 

where T denotes matrix transposition. 

For iV = 2 we have only one possibility to form a colour singlet state from two gluon 
operators 

V, v (x, y) = tr(A M (x), A„(y)) = ^5 ab Al(x)A b v (y) . (2.77) 

This state is invariant under charge conjugation (|2.76j) so that it has C = +1. Therefore, for 
two gluons we have only the Pomeron states. 

In the N = 3 case they are two ways of constructing a colour singlet state. Firstly, one can 
build 

V^ p (x,y,z) = -itr([A^x),A v (y)]A v (y)) = l -f abc Al(x)A b MA c p (z) (2.78) 
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with the total antisymmetric structure constants f a b c defined via the Lie algebra of SU(N C ): 

[t a ,t b ] =if abc t c . (2.79) 

Using (|2.76L we find that V IJiU p(x,y, z) have also C = +1 so for N = 3 the state (I2.78j) gives a 
subleading correction to the N = 2 Pomeron state. 

The other possibility to form colour singlet states is to use the totally symmetric constant 

d abc = 2 [tr (t a t b t c ) + tr (t c t b t a )] . (2.80) 

In this way we obtain a state 

0^ p (x,y,z) = tr({A^(x), A„(y)} A v (y)) = \d^A#x)Al{y)A%z) . (2.81) 

Similarly, applying (|2.76jl we find that (I2.81j) is the odd under C— parity. Thus, it is the leading 
contribution to the odderon state. This state only appears for the gauge groups SU (N c ) with 
iV c >2. 

In order to find the eigenstates of the Hamiltonian 1)2. 22j) we perform the multi-colour limit 
(j2.26jl . In this limit the colour factor of the Hamiltonian is reduced to a trivial delta- function. 
Thus, a solution to the Schrodinger equation (|2.34|) does not have a colour factor. On the 
other hand, the Reggeon wave-function should be invariant under Bose symmetry. Adding a 
proper colour factor into the Reggeon wave-functions we are able to restore Bose symmetry. 
In the multi-colour limit the Bose symmetry reduce into the cycle and mirror permutation 
symmetries (12.681) . Thus, we have states which under the mirror permutations are either odd 
or even. For example in the N = 3 case it turns out that adding colour factor: corresponding 
to the antisymmetric constant f a b c for the odd-mirror states and the symmetric one d a b c for the 
even-mirror states, we are able to restore Bose symmetry. As we have said in (I2.78J) and (|2.81j) 
f a bc corresponds to the Pomeron states while d a bc is related to the odderon states. One may 
expect similar relationship for iV > 3 states. Therefore, in order to check C— parity of a given 
state we need to study its parity under the mirror permutation M. The states ^ satisfying 
M\& = — $ are the Pomeron states whereas states for which M\l/ = are the odderon states. 
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Chapter 3 
Baxter Q-operator 



3.1 Definition of the Baxter Q-operator 



The Schrodinger equation (|2.34|) may be solved applying the powerful method of the Baxter 
Q-operator [30J. This operator depends on two complex spectral parameters u, u and in the 
following will be denoted as Q(u,u). This operator has to satisfy the following relations 



Commutativity: 
Q — t relation: 



(u,u),Q(v,v)\ = 0, 
[i N (u),Q(u,u)) = t N (u),Q(u,u) 







Baxter equation: 

i N (u)Q(u, u) = {u + is) N Q(u + i,u) + (u- is) N Q(u - i, u) 



(3.1) 
(3.2) 

(3.3) 



t N (u)Q(u, u) = (u + is) N Q(u, u + i) + (u- is) N Q(u, u-i), (3.4) 

where i N (u) and tjy(u) are the auxiliary transfer matrices (|2.45j) . According to ()3.2|l the Baxter 
Q(u, w)-operator and the auxiliary transfer matrices as well as the Hamiltonian ()2.26|) share 
the common set of the eigenfunctions 

Q(u, z 2 , ■ ■ ■ , z N ) = Q q> q{u, u)^ q z(zi, z 2 ,..., z N ) . (3.5) 

The eigenvalues of the Q-operator satisfy the same Baxter equation (I3.3JI and (I3.4jl with the 
auxiliary transfer matrices replaced by their corresponding eigenvalues. 

In the paper [HHJ the Q-operator was constructed as an N— fold integral operator 

Q(u,u)^(z 1 ,z 2 ,...,z N ) = J d 2 w 1 J d 2 w 2 ... 

...J d 2 w N Q u ,u(zi, z 2 ,..., z N \wi, w 2 ,..., w N )^(w h w 2 ,..., w N ) , (3.6) 
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where the integrations are performed over two-dimensional planes. The integral kernel in 
()3.6|) takes two different forms: 



N r _ n 1 — 2s 

=a(2-2s,s + iu,s- ]J - ^ (3.7) 

t=l l W k - M [Wk ~ Zk+l\ 



and 



«S<*M = n 7 ^^r^ f= • (3-8) 

fc =1 [ z k ~ Wk\ [Zk ~ Wk+l\ 

which appear to be equivalent [38j. In Eqs. (13.711 and 1)3.8)1 a(. . .) factorizes as 

j / 1 — \ 

a(a, /?,...) = a(a)a(/?) . . . and a(a) = — — (3.9) 

and a is an anti-holomorphic partner of a satisfying a — aE Z. Moreover, the two-dimensional 
propagators are defined as 

[z k - w k ]~ a = (z k - w k y a (z k - w k y a . (3.10) 

In order for the Baxter Q— operators to be well defined, ()3.7|) and 1)3 .8)1 . should be single- valued 
functions. In this way we can find that the spectral parameters u and u have to satisfy the 
condition 

i{u — u)=n (3-H) 

with n being an integer. 

The Baxter Q-operator has a defined pole structure. For Q^ + \u,u) we have an infinite set 
of poles of the order not higher than N situated at 

{u+=i(s-m),u+=i(s-m)}; { u~ = -i(s - m),u^ = -i(s - m)} (3.12) 

with m,m = 1,2,.... The behaviour of Q qj q(u,u), i.e. an eigenvalue of QW(u,«), in the 
vicinity of the pole at m = m = 1 can be parameterized as 



Qq,g( u i +e,uf + e) = ^(q, q) 



1 ^E ± (q,q) 

e N + € N-1 



(3.13) 



The functions R^i^q^q) fix an overall normalization of the Baxter operator, while the residue 
functions E ± (q,q) define the energy of the system (see Eqs. 1)3.21)1 an d 1)3.24)1 below). It has 
also specified asymptotic behaviour. For |ImA| < 1/2 and ReA — > oo 

Q q>7j {\ - in/2, A + in/2) ~ e ie ^ X h+h-N( s -s) + e ~ie h ( q ,q) x i-h+i-h-N{s-s) ? ^ u ^ 

where Oh is a phase which should not be confused with quasimomentum ^(g, <?). 
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3.2 Observables 



The Hamiltonian (|2.31|) may be written in terms of the Baxter Q-operator 



Ti-N = e N + i— In (Q|( + ) (u + is,u + is) 
du 



— ( i—]nQ( + \u — is,u — is) 



u=Q 



where the additive normalization constant is given as 

e N = 2N Re[^(2s) + ip(2 - 2s) - 2^(1)] . 
Applying to (|3.15j) the eigenstate we obtain the energy 

E N(q, q) = £n + i-j- In [Q q ,q(u + is,u + is) {Q q4 {u -is,u- is))*] 
or equivalently 



u=0 



E N (q,q) 



Im In 

du 



u 2N Q q ,q(u + i(l-s),u + i{l - s)) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



q,-q 



[u + i(l — s),u + i(l - s)) 



u=0 



where Q q q(u,u) = Q q + ^(u,u) is eigenvalue of the Q+(u, u) operator, while 

±q = (q 2 ,±q3,...,(±) N q N ) 

are the conformal charges. 

It is also possible to rewrite the quasimomentum operator in terms of Q+(u, u): 

IJisJs) 



(3.19) 







N 



-i In P = % In 



(3.20) 

\l+\—is, -is) 

Moreover, using the mirror permutation ()2.72|) one finds the following parity relations for the 
residue functions R + (q,q) defined in (|3.13jl : 

R + (q,q)/R + (-q, -q) = e 2i9 ^ 



and for the eigenvalues of the Baxter operator: 



Q 



-u, —u) 



Q-q,-q(u, U) 



(3.21) 



(3.22) 



where — q = (q 2 , —q3, ■ ■ ■ , (—l) n q n ) and similarly for q. Examining the behaviour of (|3.22|) 
around the pole at u = uf and u = uf and making use of Eq. (|3.13jl one gets 

R ± (q,q) = (_l)"e**fe'*> R*(-q, -q) , E ± (q,q) = -E*(-q, -q) . (3.23) 

To obtain the expression for the energy E N (q,q), we apply (|3.17|) and replace the function 
Q q ,q(u i i(l — s),tii i(l — s)) by its pole expansion (|3.13l) . Then, applying the second relation 
in (I3.23L one finds 

E N (q,q) = E + (-q, -q) + (E + (q, q))* = Re [E + (-q, -q) + E + (q, q)] , (3.24) 

where the last relation follows from hermiticity of the Hamiltonian (I2.32jl . We conclude from 
Eqs. ()3.24|) and ()3.13|) . that in order to find the energy E N (q,q) 7 one has to calculate the residue 
of Qg t q(u,u) at the (N — l)-th order pole at u = i(s — 1) and u = i(s — 1). 
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3.3 Construction of the eigenfunction 



The Hamiltonian eigenstate ^ p,{ q ,q}{z) is a common state of the total set of the integrals of 
motion, p and {q, q} as well as the Baxter Q-operator. Thanks to the method of the Separation 
of Variables (SoV) developed by Sklyanin |HJ EE] we can write the eigenstate using separated 
coordinates x — (xi, . . . , x^-i) as 



^P,{q,q}(z) = d x/j,(x 1 ,...,x N - 1 )U^ u ...^ N _ 1 (z 1 ,... } z N )(^ q ^(x 1} ...,x N ^ 1 ))*, (3.25) 



where is the kernel of the unitary operator while 



N-l 



. . . , x N ^y = e ^^)/ 2 n f — r(s + tXk)T{ * - tXk) 



k=l 



s + ix k )T(l - s - ix k ) 



N 



Qq,q(x k ,X k ). (3.26) 



The functions Q qt q(x k , Xk) are eigenstates of the Baxter Q-operator. In contrast to the Zi 
(zi,Zi) - coordinates, the allowed values of separated coordinates are 



m k _ in k 

x k — v k — , Xk — Vk H — ^~ 



(3.27) 



with n k integer and v k real. Integration over the space of separated variables implies summation 
over integer n k and integration over continuous v k 



N-l / oo 
fc=l \n k =-oo^ 00 



dvi. 



2tt 



_At2 N-l 



j>k 



\Xf£ Xj 



3\ ' 



(3.28) 



where \x k — Xj\ 2 = [y k — Uj) 2 + (n h — rij) 2 /A. 
The integral kernel can be written as 



U^fa, ...,z N ) = c N {x){p 2 ) {N -^l 2 I d 2 w N e 2 ^U^{z x , ...,z N ; w N ) 

where 2p • wn = V w n + V wn, 

U$(zx,..., z N ; w N ) = Aj'* lj(ft) AjJ;_2',(s? a *y • • • A i,(L-i) ( : 1 : A I ,r - x 5 
for even N, and 



(3.29) 



(3.30) 



U*(zi 



gj^l, • • • , X.JV; 



z N ;w N 



\( S , S ) aC 1 -^ 1 "*) AC 1 - 5 . 1 -*) 

ly N-l,(x 1 ) 1)L N~2,(x2) ■ ■ ■ lv l,{S N -i) 



IZi, . . .,Z N \W N 



(3.31) 



for odd N. Here the convolution involves the product of (N — 1) functions A N _ k ^ k ^ with 
alternating spins (s,s) and (1 — s, 1 — s). They are defined as 



A 



iV-n,(a-)(^> • • • ' ZNWn+1, ■ ■ ■ , Vn) = [Zl ~ J/ 2 

(N-l 
k=n+l 



1 —x+iu 



I [ [z k - y k ] x lu [z k - y k+1 ] x+m ) [z N - y N ] x tu , (3.32) 
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where the convolution [A N -k t (£ k )A N _ k+1 ^ k _^] contains (N — k) two-dimensional integrals. The 
coefficient cjy(x) is given for N > 3 

[(AT-l)/2] [JV/2-1] 

c N (x) = Y\ ( a ( s + ix 2k,s - ix 2 k)) N ~ k ]| (a(s + ix 2k +i,s - ix 2k+1 )) k (3.33) 

k=l k=l 

while the sums goes over integer numbers lower than upper limit. For iV = 2 we have c 2 (xi) = 1. 
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Chapter 4 

WKB approximation 



4.1 WKB approximation for Reggeons 

The WKB approach was presented in a series of papers [101 H3 EI] • We will use this method to 
find an approximation set of conformal charges as well as to understand the structure of the q^s 
values. In the representation of the separated coordinates (J3.25)) - 1)3.27)) the Baxter equations 
(13.3)1 . (13.41) play the role of one-dimensional Schrodinger equations where the eigenfunctions 
are the eigenvalues of the Baxter operator 1)3.26)1 . These equations may be solved in the quasi- 
classical limit [21] which corresponds to large values of the conformal charges, and q k , 1)2.45)1 . 

In order to apply the WKB limit we introduce an auxiliary parameter rj and rescale simul- 
taneously the coordinates, x — > x/r], and the charges. Thus we define 

t N (x) = r] N t N (x/ri) = 2x N + q 2 x N ~ 2 + ... + q N x N (4.1) 

with qk = qkV k = 0{rf) as rj — > 0. The small parameter r] appears in the Baxter equation 

(x - ir]s) N Q((x - ir))/rj) + (x + ir]s) N Q((x + irj)/rj) = Q(x/rj) t N (x) (4.2) 

which allows us to perform the WKB expansion 

Qix/rj) = exp ( - ! dxS'{x)\ with S(x) = S (x) + rjS^x) + r] 2 S 2 (x) + C(r/ 3 ), (4.3) 

where S'(x) = dS(x)/dx and Xo is an arbitrary reference point. Substituting 1)4.3)1 into the 
Baxter equation (14. 2 j) we obtain relations 

2 cosh ^(x) = S[(x) = ^(lnsinh^(a;)) , + (4.4) 

In this way one can obtain further relations for higher Sk(x). 
The first relation in (14. 4j) can be rewritten as 

y 2 (x)=t N \x)-4x 2N , (4.5) 
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where 

y(x) = 2x sinhp x with p x = Sq(x) (4.6) 
is a hyperelliptic curve. Solving this relation ()4.4|) we obtain the leading term of ()4.3|) 



r<e dx 

S (x) = I dxp x = I ^r(Nt N (x) - zAn(x)) + xp. 



(4.7) 

x 



The hyperelliptic curve, y(x), as well as S'(x), are double- valued functions on the complex 
x-plane. Two different branches of S' (x) will be denoted by (±), so we have S' ±(x). To specify 
them one makes cuts on the x-plane in an arbitrary way between the 2(N — 1) branching points 
Oj where yioj) = 0. Taking (|4.4j) and (14.71) we have the first non-leading correction 



x — cr,- x 



which is a single- valued function on the complex x— plane. 

Combining together (|4.7I) and (I4.8J) we obtain two different solutions to the holomorphic 
Baxter equation (14. 2 j) 

Q±{x/rj) = exp ( - / dxS' ± (x)) , where ^ = S ,± + r/^(x) + C(r/ 2 ) (4.9) 
W Jx J 

with asymptotics at x — > oo: 

Q+(x/^) ~ x 1 -^, Q-(x/rj) ~ x""^. (4.10) 
Similarly, one gets solutions in the anti-holomorphic sector with 

Q + (x/ri) ~ x 1 "^, g_(x/r^) ~ x 7 ^, (4.11) 

as x — > oo. 

4.2 Quantization conditions with WKB 

The quasi-classical solution of the Baxter equation (|4.2j) is given as a bilinear combination of 
the chiral solutions Q±(x/r}) and Q±(x/rj): 

Q{x/r),x/rj) = c + Q + (x/7])Q + (x/r]) + c-Q-(x/rj)Q_(x/r}). (4-12) 

The cross-terms Q±Q T do not enter (I4.12|) since they do not satisfy ()3.14j) . 

The function ()4.12|) should not depend on reference point x . This gives condition 

( i f x '° i r'° _— ' _ \ 

c ±( x 'o) — c±exp - / dxS'Ax) H — / dxSAx) (4.13) 

\VJxo VJxo J 

with x"o = Xq. Additionally, taking into account that 

c + (a j )=c4a j ) with j = l,2,...,2(iV-l) (4.14) 
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one can derive [21] the quantization conditions 

Re * dxS' (x)/7] = irizk-i + 0(v), Re <j> dxS' (x)/r] = Txt 2k + 0(r}), (4.15) 

where k = 1, . . . , N — 2 and £ = (Zi, . . . , I2N-4) being integer and the integrals are performed 
over cycles around branching points CTj, Fig. l4.11 The relations (j4.15j) give the WKB quantization 
conditions for the integrals of motion, q k and q k . 




Figure 4.1: The Riemann surface Tn with integration cycles, ctj and f3j. The dashed lines 
represent the contours on the lower sheet [21]. The contour lp-p+ used in ()4.16|) . goes from 
the point x = on the lower sheet to x = on the upper sheet. 



4.3 Lattice structure of the conformal charges spectrum 

Additionally [21], one can obtain the quasimomentum, (|3.2()l) as 

f 2tt 
N = -2Ref dxS' (x)/r] = —t (mod27r) (4.16) 

and the energy (|3.17jl in the WKB approximation 

E^ s) =41n2 + 2Re (^0-- s -i\ k ) +ip(s -i\ k ) -2ip(l)) 

lmx k >o 

+ 2Re &(l-s + iX k )+ip(s + iX k ) -2^(1)) , (4.17) 

ImA fc <o 

where t N (X k ) = for = 1, . . . , N. 
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The further analysis of (|4.15|) and (|4.16l) where we fix rj = (gAr/4)~ 1//Ar , which is equivalent 
to expansion for large conformal charges, gives the system of (N — 2) equations for the integrals 
of motion 



UN ■ N B Q, i) - {u 2 u N )* ■ B (± ^) = vr £ e^^ 1 ^ „ 4 , 
(u N+ i- m u N ) -bQ,^)- (wiv+i-mWiv)* • B Q, = 71 ^ e-^- 1 



(4.18) 



N 

)m/N 



with m — 2, . . . , N — 2 and 



Un = 4 V 4 



(fp forn = 2,...,iV-2, «,= (f)^. (4.19) 



The Euler beta-function B(x, y) = T(x)T(y)/T(x + ?/) and are integer numbers related to 
the quasimomentum ()4.16ft by 



A' 



= -^n k (modiV). (4.20) 



k=l 



The equations (I4.18JI give only (N — 1) conditions so the system is underdetermined. 

Solving (|4.18j) one can obtain the quantized values of the highest conformal charge q^ as 

<lT = -^^ 1 Q(n) (l + f j^cot(n/N)\Q(n)r + 0(Q(n)n^ , (4.21) 



where 

N 



Q(n) = ^n k e -^-i)/N _ ^ ^ 

k=l 

with n = {m, . . . ,Un} integer forming a lattice structure on the spectrum of qx (14.21^ . 
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Chapter 5 

Eigenfunctions for the conformal charges 



In this Chapter we present the general method of constructing the ansatzes for eigenfunction 
of the conformal charges %■ These ansatzes have been commonly used before |4*5l l4T? | IHU | I5T]. 
Here we systematize the knowledge about this ansatzes and extend our calculations to an 
arbitrary complex spin s. Moreover, we derive differential eigenequations for the conformal 
charges. Next, we show solutions to the differential equation for N = 3 and s = using the 
series method and present derivation of the quantization conditions for the integrals of motion 
qk- Finally, we resum obtained series solutions for the q 3 = case. 

5.1 Solution of the q_i eigenproblem 

In this Section we solve the q 2 eigenproblem. This gives us not only the solution for the N = 2 
Pomeron problem but also provides us with a structure of the ansatz for the eigenfunctions of 
the conformal charges, qk for N > 2. 

As one can see from Eq. (I2.54JL the eigenfunctions (|2.55j) of the Hamiltonian (|2.31j) depend 
on differences of coordinates, z\ — z . They may contain two types of factors. The first ones are 
arbitrary functions of (N — 2) variables, Xi = ; which are SL(2, C)— invariant. 

They are defined in Appendix A. The other ones are some products of the coordinate differences 
that satisfy the transformation law (|2.55jl . Let us consider the general product 

N 




(5.1) 



*2>«1=0 



where z i2il = z i2 — z ix and k ili2 is arbitrary. 

It is convenient to parameterize the eigenvalues of the lowest conformal charge as 



N 




(5.2) 
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where h and Sj are the conformal spins l|2.56j) . Now one can act with the q 2 operator (j2.65j) on 
the product of holomorphic coordinates 1)5. lj) : 

N N 

Q2 IJ ( z nh) kili2 = 12 IJ (^u) fcl2n - (5-3) 

«2>H=0 i2>H=0 

In this way we obtain the equation for which can be easily solved. 

For each N we obtain two sets of solutions: the first ones, with the z coordinate, and the 
second ones, without the z coordinate. The latter solutions may be calculated by integrating 
the first ones over zq with p = with the help of 1)2.54)1 . They don't transform as (12.551) so they 
will not be further considered. The integration with an arbitrary p 1)2.541) for two Reggeons was 
performed in [H2]. In the first set we obtain two groups of solutions related by the symmetry 
h — > 1 — h, which comes from that the form of the q 2 eigenvalue ()5.2|) . More detailed studies of 
this problem are included in Appendix B. 

As a result we obtain an ansatz for the eigenfunction. For an arbitrary N it has a form 

h— 81— S2 — s N 



{z w )^(z 20 y^ . . . {z N0 )* s " \z 10 z 30 J 



#(zio, z 20 , z N0 ) = -_ ( ^— ) F(x 1 ,x 2 , x N „ 2 ) . 



(5.4) 

For N particles we have N — 2 invariant independent variables, Xi. For the homogeneous spin, 
i.e. with si = s 2 = ■ ■ ■ = sn = s, chains we have 



1 ( ^_\ 



h-Ns 



ty(z 10 ,z 20} ... ,z NQ ) = — F(x 1 ,x 2 ,...,x N - 2 ) ■ (5.5) 

{z 10 Z 20 . • -Z m y s \Z W Z 3 qJ 

5.2 Various ansatzes 

Using the above ansatzes we can construct other ansatzes which are equivalent to the original 
ones. For example, for N = 3 [51j, 



h—ai— S2-S3 



(zio, z 20 , z 30 ) = - — ^— — — — — F(x) , (5.6) 

(Zl0) 2si {Z 20 ) 2s 2 (Z 3 0) 2S3 V ^10^30 / 

{ (x-l) 2 \ V3-si/3-s 2 /3-s 3 /3 

where x = x 2 we substitute F(x) = G(x) ( v _J J what gives an ansatz 

(h _ £1 _ £2 _ £3 
Z 3 1 Z\ 2 Z 23 \ 3 3 3 3 

(^ioA^oA^o) 2 / 

or equivalently 



s 3 

^2 \ -f 



*( Zl0 ,z 2Q ,z 3Q ) = ( (^ 2 ^°) 2 V ( (^o) 2 (^o) 2 y ( (z 20 n Zw ) 

" \{z w ) 4 z 12 z 23 z 31 J \{z 20 ) 4 z 12 z 23 z 31 ) \{z 30 ) 4 z 12 z 23 z 31 



x{ uJ(Zn^ ] G(x) - (5i 
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As we can see we obtained a totally symmetric ansatz which is equivalent to original one. For 
the homogeneous spin chains we have 



^ Zm) = few Uo?Sw) ' G{x) ■ <5 ' 9) 

Both ansatzes have advantages and disadvantages. The symmetric one 1)5.8)1 is appropriate 
if we want to deal with the particle symmetries. The original one 1)5.6)1 has a simpler structure, 
it contains powers of h (not h/3), so we can use it when we want to construct proper single- 
valuedness conditions. One can easily notice that z h z h is single- valued because h — h = rih € Z. 

Similarly, we can go from the original ansatz (15.6)1 to an ansatz with different permutation 
of particles 

tf(*io, z 20 , z m ) = ( Z , n )2 ai ( z j)2s 2(z . mr js 3 (i) F(x) = 

h-s 1 -s 2 -s 3 

H(x) (5.10) 





231 \ 




^10^30 / 


i 


[ Z12 


(*10) 2S 1(Z2C)P2(*30) 2S 3 


\ Z1QZ20 


1 


( Z32 


(«10) 2S 1(Z20)^2( Z30 )^3 


\Z20Z3O 



h— SI— S2— S3 

I •/(*) • 

Now we can generalize our ansatz for different number of particles, N. In this case we have 

(zm) 2 (z2 ) 2 ...(z NO ) 2 \^ ( (z 10 f(z2 y 2 ...(z N0 ) 2 \% 



V{Z W ,Z 20 ,...,Z m ) - ^ {Z10 )2» Zl2Z23 ...z N1 ) ^ 



(2l0) 2JV 2i2223---^JVl / \ (22C)) 2JV Z12223---2JV1 



' (Z 10 ) 2 (Z2Q) 2 ...(Z N0 ) 2 \ N ( Z 12 Z 23 ...Z N1 \ " r( \ 

(zno) 2n z 12 Z2 3 ...z n1 J \(z w ) 2 (z2o) 2 -(z NO ) 2 J ^^1)^2, • • • , • 

(5.11) 

And for the homogeneous spin chains, Sj = s, it has a form 



i, x 2) • • • , a;Ar-2j • (o.lzj 

{z 12 z 2 3 ■ ■ ■ z N1 y \{ziqY{z 2 qY . . . {z NO yj 

To sum up, we have many equivalent ansatzes. Here we presented two of them, 1)5.6)1 and 
(15.8)1 . They are most frequently used due to their simplicity and symmetry properties. 



5.3 Solutions of the eigenproblem for N = 3 

In this Section we solve the ^3 eigenproblem for three Reggeons. We derive differential equations 
for the integral of motion ^3 for an arbitrary complex spin s and solve them for s = s = making 
use of the series method. We also construct the quantization conditions for q 3 and in the end 
we resum the series for zero-modes of g 3 . We obtain the known solution 1)5.44)1 as well as some 
others with Log(x)— terms (15.511) . 



5.3.1 Derivation of differential equations for N = 3 

In the previous Section we solved the eigenproblem for q 2 . However, for N particles there are 
N — 1 conformal charges which commute with the Hamiltonian. Thus, for iV = 3 we have q 2 
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and q 3 . In order to to solve the eigenproblem for q 3 we can act with q 3 on our ansatz 



1 ( z 32 \ 



h-S!-S2-S 3 



Zx ' Zm) = www^ F(x) <5 - 13) 

and derive a differential equation for F(x) with x = x 2 = ^ Zl ~ Z2 ^ Z3 ~ z °\ . Thus we obtain 

v ' z (21 -zo) (23 -22) 

iq 3 F(x) = -s 1 (s 1 + s 2 + s 3 — h)(h — 1 — s x — s 2 + s 3 — 2hx + 2(1 + si)x)F(x) + 

+ ((s 1 + s 2 )(l-/i + s 1 + s 2 -s 3 ) + 
-(2 + /i 2 + 7si + 3s 2 + (si + s 2 )(5si + s 2 ) - h( 3 + 6s i + 2s s) + s 3 - s§)x+ 
+ (h 2 + 2(1 + s 2 + s 3 ) - h{3 + 6 Sl + s 2 + s 3 ) + si(7 + 5si + 3s 2 + 3s 3 ))x 2 )F , (x)+ l °' 1 
-(a; - l)z(2 + 2 Sl + 2s 2 - s 3 - (4 + 4 Sl + s 2 + s 3 )x + h{2x - l))F"{x) + 

+ (x-l) 2 x 2 F^(x). 

Using various ansatzes we obtain equivalent differential equations. 
For the homogeneous chain our equation looks like 

iq 3 F(x) = (3s - h)(h - 1 - s)(l - 2x)F(x) + {{(h - 2)(h - l)(x - l)x + s 2 (2 + ll(x - l)x) + 
+s(2 - 2h(l - 2xf + ll(x - l)x))F'(x) + (2 + h - 3s)(l - x)x(2x - l)F"(x) + 

+ (x-l) 2 x 2 F^(x). 

(5.15) 

From the QCD point of view the most interesting cases are for s = 0: 

iq 3 F(x) = (h - l)(h - 2)x{x - l)F'(x) + (h-2)(x- l)x{l - 2x)F"(x) + x 2 (x - l) 2 F^(x) 

(5.16) 

and for s = 1 

iq 3 F(x) = {h - 3)(h - 2) (2a; - l)F(x) + 
+ ((4 -2h-(h- 8)(h - 3)x + (h- 8)(h - 3)x 2 )F'{x) + (5.17) 
+ (5 - h)(x - l)x(2x - l)F"(x) + (x- l) 2 x 2 F^\x) . 

The first such solution was derived and found numerically in [T8] . 



5.3.2 Solutions for s = 

Now, we will solve the above equations by the series method. All these equations have three 
regular singular points at x — 0, x — 1 and x = 00. Let us take an ansatz 1)5. 13j) which for 
s = can be rewritten as 

qr(z 10 ,z 20 ,z 30 ) =w h F(x), (5.18) 



where w is defined similarly to (|A.4j) with z 3 , z 2 and z . Thus, we obtained ()5.1 

iq 3 F(x) = (h- l)(h - 2)x(x - l)F'(x) + (x - l)x(h - 2)(1 - 2x)F"(x) + x 2 (x - l) 2 F^(x) . 

(5.19) 
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This equation can be used to obtain solution around x = + . To generate solutions around 
other singular points we exchange variables. For the case x = l~ we can use a substitution 

x = 1 — y: 

iq 3 G(y) = (h-l)(h- 2)1/(1 - y)G'(y) + (h-2)(l- y)y(l - 2y)G"{y) - y\y - l) 2 G^(y). 

(5.20) 

Moreover, for x = 1 + we can use x — y + 1: 

q 3 G(y) = i (y(y + l)(h- \){h - 2)G'{y) - (y + l)y(h - 2)(1 + 2y)G"(y) + y\y + \fG^{y)) 

(5.21) 

and for x = oo~ we have x = 1/y: 

q 3 G(y) =i((y- l)(h + l)(h- 2y)G'(y) + (y - l)y(2(h + 1) - (h + A)y)G"{y) 

+y { y -iyG^\y)) . (5.22) 

The upper-scripts plus and minus correspond to case where Re[x] is above and below the 
singular point, respectively. 

5.3.3 Wave- function for s = s = around x = x = + 

In order to obtain the full-complex solution containing the holomorphic and anti-holomorphic 
parts we have to glue together solutions from these parts: 

^({Zi}, {Zi}) = %({^}) T • A®(h,h, q 3 , q 3 ) ■ u q ({z t }) , (5.23) 

where we use a (N xN) mixing-matrix, ^4^, [HI which does not depend on particle coordinates 
but only on q = {q2,qs}- From the QCD point of view we have two possibilities of gluing 
solutions: (s = 0,s = 0) [2Ij and (s — 0, s = 1) [3HJ - These two cases are equivalent except 
zero modes of the highest conformal charge q^. Let us consider the first case, s = s = 0. 

The conformal charges are related by conditions h = 1 — h* and q k = q^. The wave-function 
has to be single-valued. This condition defines the structure of the mixing-matrix. 

For h ^ Z and q 3 ^ we have solutions of the following type 

^l(^) % y~] T> _ n Qn.n % > 

U 2 {x) = ^ 1 Y,7=0 a n,T2 Xn > ( 5 - 24 ) 

u 3 (x) = x J2™ =0 b n:r3 x n + x 1 J2™ =0 a n!r2 x n Log(x) 

and similarly in the anti-holomorphic sector. The coefficient recurrence relations for the a n , r . 
are given in Appendix C. One can notice that x a x b is single- valued only if a — b G Z. Moreover 
we have also terms with Log(x) which have to give in a sum Log(xx). So in this case we have 
a mixing matrix of the form 



A^(h,h,q 3 ,q 3 ) 



a 

/3 7 
7 



(5.25) 
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where a, (3, 7 are arbitrary. In the above matrix we have Ai 2 = A% 3 = A 2 \ = A 3 i = in 
order to eliminate multi-valuedness coming from the power-terms, A 23 = A 32 to obtain single- 
valuedness in Log(x)— terms and A 33 = because the term Log(x)Log(x) is not single- valued 
on the x— plane. 

In the case of q 3 = and h G" {0, 1} we don't have any Log(x)— terms so the mixing matrix 
looks like 1 

(3 

A iQ \h,h,q 3 = 0,q 3 = 0) = a 5 . (5.26) 

For q 3 7^ and h G Z we have a solution with all power in x integer and solutions without 
Log, with one-Log and with double-Log. The structure of the matrix is 



A^(heZ,heZ,q 3 ,q 3 ) 



a (3 7 
(3 2 7 
7 



(5.27) 



In the last case for q 3 = and h G {0, 1} we have all powers of x integer and the third 
solution with one-Log term. The matrix has a form 



A®(h G {0,1}, h G {0,1}, q 3 = 0,q 3 = 0) 



a 7 
(3 5 e 
e 



(5.28) 



5.3.4 Wave-function for s = s = around x = x = 1 and x = x = 00 

Around the other singular point we construct the wave-function exactly in the same way ob- 
taining matrices A^- 1 ^(h,h,q 3 ,q 3 ), A^ 1+ \h,h,q 3 ,q 3 ) (around 1) and A^°° \h,h,q 3 ,q 3 ) (around 
00). 

Thus, we have the wave- function similar to (|5.23j) . For h G" Z and q 3 7^ we take solutions 



m(x) = (1 — x) Z)£L an,ri(l-aO n , 

u 2 (x) = (l-^E^o^a-z)", (5.29) 

u 3 {x) = (1 - x)° E: =0 & ffi , r3 (l - x)" + (1 - j) 1 L M =0 V,(l - ^) n Log(l " ^) 

and similarly in the anti-holomorphic sector. Our the wave-functions have to be single-valued. 
Thus, the mixing matrices take the following forms 



A^(h,h,q 3 ,q 3 ) 



a 

(3 7 
7 



A^\h,h,q 3 = 0,q 3 = 0) 



(3 
a 5 
e 7 



(5.30) 



Greek variables in each A— matrix have different numerical values 
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A {1 \h eZ,he Z,q 3 ,q 3 ) 



a /3 7 
(3 2 7 
7 



(5.31) 



A^\he{0,l},he{0,l},q 3 = 0,q 3 = 0) 



a 7 

5 e • (5.32) 
e 

Similarly, we proceed around a; = oo~. For h ^ Z and ^3 7^ we have solutions of type 

= (l/^) Er=oVn(V^)", 

«a0«0 = (l/*) 1 - fc E^ ) a»,™(V*)". 
u 3 (z) = (l/a;)^ 1 £~ =0 & n , r3 x n + {l/xf- h J2™=o a n , r2 x n Log(x) 

and similarly in the anti-holomorphic region. In this case we have the matrix 



(5.33) 



A^\h,h,q 3 ,q 3 ) 



a 

(3 7 
7 



A^\h,h,q 3 = 0,q 3 = 0) 



poo 

a 5 

e 7 



(5.34) 



A^\heZ,heZ,q 3 ,q 3 ) 



a (3 7 
P 2 7 
7 



(5.35) 



A( oo )(/ i G{0,l},/iG{0,l},g3 = 0,g 3 = 0) 



070 
P 5 e 
£ 



(5.36) 



5.3.5 Transition matrices between solutions around different poles 

The above solutions around x = 0,1, 00 have a convergence radius equal to the difference 
between the two singular points: the points around which the solution is defined and the 
nearest of the remaining two. In order to define a global solution which is convergent in the 
entire complex plane we have to glue the solutions defined around different singular points. This 
can be done by expanding one set of solutions in terms of the other solutions in the overlap 
region of the two considered solutions. Thus, in the overlap region we can define the transition 
matrices A, r, where 

u^(x,q) = A(q)u^(x,q) , 
uP)(x,q) = r(q)u {co \x,q). 
Matrices, A and r, are constructed in terms of the ratios of certain determinants 
example, to calculate the matrix A we construct Wrofiskian 



(5.37) 
For 







x; 


q) 


u\ '(x; 


q) 


4 1} 


w = 


u?> 




q) 


r(l) I 
u\'[x 


q) 


uf 




u'f 


(x 


q) 


u'f\x 


q) 





x;q) 
(x;q) 



(5.38) 
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Next we construct determinants Wij which are obtained from W by replacing j-th column by 
the i-th solution around x — 0, i.e. for i — 1 and j = 2 we have 





u i \ 


T ' 

X, 


H3) 


u\ 0) ( 

u l \ 


x\ 


H3) 


a 3 \ 


T • 

X, 


13) 




w 12 = 


uf 




qs) 


uf 




13 ) 


uf 


[x; 


<?s) 


(5.39) 




u'f 


(x 


Is) 


u'f 




13 ) 




(x 


Is) 





The matrix elements A^ are given by 



A 



1.1 



w 



(5.40) 



Matrix A does not depend on x, but only on qj.. In the similar way we can get the matrices r 
and their anti-holomorphic equivalents: A, r. 

Substituting equation (j5.37j) into the wave- function (|5.23j) , one finds the following conditions 
for continuity of the matrix A(q, q): 

A T (q 3 )A^(q 3 ,q 3 )A(q 3 ) = A^(q 3 ,q 3 ) , (5.41) 
T T (q 3 )AW(q 3 ,q 3 )r(q 3 ) = A^(q 3 ,q 3 ). (5.42) 

Each Equation, (|5.41|5.42|) . consists of nine equations. Solving them numerically, we obtain 
values of parameters a, (3, 7, . . . as well as quantized values of the conformal charges, q k and q k . 
Numerical calculations tell us that the spectrum of qk obtained using this method is equivalent 
to the spectrum obtained using the Baxter Q— operator method which is presented in next 
Chapter. 



5.3.6 Additional conditions coming from the particle permutation in- 
variance 

Our states have additional symmetries: the cyclic and mirror permutation (|2.68j) . The confor- 
mal charges commute only with P. Thus, our eigenstates are hardly ever eigenstates of M, so 
they usually have mixed C-parity. 

For q 3 = we can easily resum the series solutions, see Appendix C. Let us take a case for 
h ^ {0,1}. The eigenequation for the cyclic permutation with a quasimomentum 3 (q) gives a 
following condition 

w h w>* [p + ~/(-x) h (-xf + a(x - lf(x - if + 6{-xf(x - If + e(x - lf(-xf} = 
= e^W w h w Tl (a + {3(-xf(-xf + 7(2 - lf(x - if + S(x - if + e(x - if} . 

(5.43) 

Here we have used cyclic transformation laws (|A.6j) . 

Comparing these two lines we obtain conditions: a = e~ td3 ^P, (3 = e~* 93( ^7, 7 = e~ ldz ^a 
and 5 = e = 0. One can derive that exp(3i6> 3 (g)) = 1 so 6 3 (q) = ^p- where k — 0,1,2 (k = 
for physical states). Thus we have an eigenstate of P 

y(z w ,z 20 ,z 30 ) =w h w K (^l + exp(^)(-xf(-xf + exp(i^)(x-l) h (x-lf S J , (5.44) 
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where we have omitted the normalization constant. 

Now we can act with a mirror permutation operator on (|5.44p and test its eigenequation 
(j2.68j) . Using the mirror transformations (|A.7|) . similarly to ()5.43|) . we obtain the following 
relation 



(5.45) 



w h w h (-l) nh (exp(i^) + (-x) h (-x) h + exp(i^)(x - l) h (x - l) h 
= ±w h w Tl (l + ex^{i^f-){-x) h {-xf + exp(i^)(x - l) h (x - if) 



Comparing both sides of (j5.45|) gives (-l) nh exp(z^) = ±1, (-l) nh = ±exp(i^) and 
(— l) nh = ±1 where the SX(2,C) Lorentz spins rih = h — h. These conditions are consis- 
tent with k = 0, |. Only the first case agrees with the cyclic permutation condition. As we can 
see for odd n h we have minus sign, so taking into account colour factors (— 1)^, solution (|5.45j) 
is C-even. For even rih we have plus sign thus solution is C-odd. The last case is unnormalizable 
because when x — > or x — > 1 it does not vanish so the norm, with (s = 0, s — 0), is divergent 

m 

Using the duality symmetry [KU [EEl [EH EEl IEE1 , which corresponds to h — > 1 — h, Bartels, 
Lipatov and Vacca constructed an eigenstate with q 3 = and C — — 1 

*(*io, z 20 , z 30 ) = w h ^x{l - x)x(l - x) (V 2 )(x) - 5 {2 \l -x) + ^5 (2) (-) ) . (5.46) 

This wave-function cannot be constructed using the method described here because it contains 
non-analytical functions, the Dirac delta S(x). 

Now, let us take the second wave-function with five parameters i.e. for q 3 = and h e {0, 1}. 
Choosing h — 1 we have 

u(x) = [1, (-x), (-x)Log(-x) + (x- l)Log(x - 1)] T (5.47) 

and for h = it is 

u(x) = [Log(x - 1), 1, Log(-x)] T . (5.48) 

Combining them we obtain 



^(^10,^20,^30) = w h w h (aLog(a; - 1) + (3 + 7(-x)Log(x - 1) + 5(-x)+ 
+e((—x)Log(—x) + (— a;)Log(— x) + (x — l)Log(a; — 1))) . 

Like in the previous case we write the eigenequation for the cyclic permutation 



(5.49) 



w h(=i) w h(=o) (( a _ 7 _ e )Log(x - 1) + (5 - (3) + a(-x)Log(x - 1) + + 
+£(— x)Log(— x) + a(x — l)Log(— x) + e(— x)Log(— x)+ 
+7Log(— x) + e(x — l)Log(a; — 1)) = 
= e ie3{q) w (aLog(x - 1) + f3 + 7(-a;)Log(x - 1) + 8(-x)+ 
+e((— x)Log(— 2) + (— x)Log(— x) + (x — l)Log(s — 1))) . 

(5.50) 
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Thus we get conditions: a = e~ ie3 ^(a - 7 - e), (3 = e~ ie3 ^(5 - (3), 7 = e~ Wa ^a, 5 = 
e -i6 3 (q)(_p^ = e~ ies ( q \>y - a) and e = e- Ws{q) e. We have two types of solutions. 
The first one with 9 3 (q) = when a = 7 = — e and (3 — S — 0. It has a form 

*(Si„, z 20 , z 30 ) = w ((-x)Log((-x)(-x)) + (x- l)Log((x-l)(x- 1))) . (5.51) 

We obtained in this way solutions with Log(x)— terms which have not been presented before. The 
similar expressions were shown in [27J as asymptotics of the q 3 eigenfunction. Acting with the 
mirror permutation operator on (15.511) we get 

Mw ((-x)Log((-x)(-x)) +{x- l)Log((z - l){x - 1))) = 

_ _ (o.52j 

= — w ((— x)Log((— x)(— x)) + (x — l)Log((x — l)(x — 1))) . 

We obtained minus sign so this state is also symmetric under C parity. 

Other solutions have 6 3 {q) = 27r/3,47r/3. Thus, a = 7 = e = 0and5 = -e~ ie ^ (3. The 
wave-function has a form 

*(zio, r 20 , z 3 o) = ^(1 + ie^«)) . (5.53) 
These solutions are not eigenstates of the M operator. 



5.4 Set of differential equations for N = 4 

Similarly to the case for three particles we can derive a set of differential equations for more 
particles. For N = 4 with the following ansatz 

h— s\— S2— S3— sa 



^(^10) ^20; ^30; z 4o) 



^10 2S1 ^20 2S2 ^30 2S3 ^40 2S4 



(—)' 

\ Z W Z 30/ 



xf( Xi = Z -^,x 2 = Z -^-) (5.54) 

V Z21Z40 ^10^32/ 

we have two equations for q 3 and for q±. Acting with q% we obtain 

t 0>0 F<-°>Q(x u x 2 ) + t lfl F<- 1 &(x 1 ,x 2 ) + t , 1 F^ 1 \x 1 ,x 2 ) + t 2 , F&°\x 1 ,x 2 )+ 
+t 1 , 1 F^ 1 \x 1 ,x 2 ) +t 0t2 F^ 2 \x l7 x 2 ) +t 3fi F^(x u x 2 ) +t 2 , 1 FW(x 1 ,x 2 )+ (5.55) 
+t 1 , 2 F^ 2 \x 1 , x 2 ) + t , 3 F^(x 1 , x 2 ) = , 

where coefficients are defined in Appendix D. 

Moreover, acting with g 4 on our ansatz we have 

fo,oF^°\x u x 2 ) + f lfl F^(x u x 2 ) + fo,iF^( Xl ,x 2 ) + f 2 , F^( Xl ,x 2 ) + 
+fi,iF^ 1 \x 1 ,x 2 ) + f , 2 F^(x u x 2 ) + f 3 , F^(x u x 2 ) + f 2)1 F^( Xl , x 2 )+ 
+h, 2 F^ 2 \x 1 ,x 2 ) + fo, 3 F^ 3 \x u x 2 ) + U, F(°> 2 \x u x 2 ) + f 3 , 1 F^(x 1 ,x 2 )+ { ' } "" 
+ f2,2F^( Xl ,X 2 ) + h, 3 F^ 2 \x u x 2 ) + h A F^\x llX2 ) = . 

These equations are not symmetric in x\, x 2 because the ansatz is also not symmetric. 

Equations, ()5.55|) and (|5.56j) . are very hard to solve even numerically, thus for the states 
with N > 4 we use the Q-Baxter method. 
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Chapter 6 

Quantization conditions in the Q-Baxter 
method 



In Ref. the authors describe a construction of the solution to the Baxter equations, (13. 3 j) 
and (|3.4|) . which satisfies additionally the conditions, (|3.12|) and ()3.14ft . This can be done by 
means of the following integral representation for Q q>q (u, u) 

Q qrq (u,u)= [ fl z - lu z~ m Q(z,z), (6.1) 

where we integrate over the two-dimensional z— plane with z = z* and Q(z, z) depends on 
{q,q}. The advantages of this ansatz are: 

• the functional Baxter equation on Q q<q (u,u) is transformed into the N— th order differ- 
ential equation for the function Q(z,~z) 



N 

N-k 



Z'izdzfz 1 " + z-\zd z ) N z s - 1 - 2(zd z f -J2^ k Qk(zd ; 



k=2 



Q(z,z)=0. (6.2) 



A similar equation holds in the anti-holomorphic sector with s and qu replaced by s = 1— s* 
and q k = q^, respectively. 

• the condition (|3.11|) is automatically satisfied since the z— integral in the r.h.s. of (|6.ip is 
well-defined only for %{u — u) = n. 

• the remaining two conditions for the analytical properties and asymptotic behaviour of 
Q q ^q(u,u), Eqs. (I3.12jl and (|3.14|) . become equivalent to a requirement for Q(z,z = z*) to 
be a single- valued function on the complex z— plane. 

Analogically to the eigenequations for q 3 (|5.15|) . the differential equation (16 .2jl is of Fuchsian 
type. It possesses three regular singular points located at z — 0, z = 1 and z = oo. Moreover, 
it has linearly independent solutions, Q a {z). The anti-holomorphic equation has also N 
independent solutions, Q b (z). 
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Now, similarly to (I5.23j) . we construct the general expression for the function Q(z,z) as 



N 

Q(z,z) = Qa(z)C ab Q b (z), (6.3) 

a,b=l 

where C ab is an arbitrary mixing matrix. The functions Q a (z) and Q b (z) have a nontrivial 
monodromy 1 around three singular points, z, z = 0, 1 and oo. In order to be well-defined on 
the whole plane, functions Q(z,z = z*) should be single-valued and their monodromy should 
cancel in the r.h.s. of (|6.3jh This condition allows us to determine the values of the mixing 
coefficients, C a b, and also to calculate the quantized values of the conformal charges q k . 

The differential equation (|6.2|) is also symmetric under the transformation z — > 1/z and 
qk — > (— l) k qk- This property is related to Eq. (|3.22l) and leads to 

Q qrq (z, z) = e i6 »^ Q^i/g, 1/z) , (6.4) 

where ±q = (q2, ±^3, . . . , (±) N qN) denotes the integrals of motion corresponding to the function 
Q(z,z). The above formula allows us to define the solution Q(z,z) around z = oo from the 
solution at z = 0. Thus, applying ()6.4|) we are able to find Q(z,z) and analytically continue it 
to the whole z— plane. 

6.1 Solution around z = 

We find a solution Q(z) ~ z a by the series method. The indicial equation for the solution of 
Eq. (|6.2jl around z = reads as follows 

(0-1 + 5)^ = 0. (6.5) 

and the solution, a = 1 — s, is iV— times degenerate. This leads to terms ~ Log fc (z) with 
k < N — 1. We define the fundamental set of linearly independent solutions to (16 .2|) around 

z = as 



m—l 



u^Log^iz) + Y J C k m - l u k+1 (z)hog m - k -\z) 
fc=l 



(6.6) 



with 2 < m < N and where for the later convenience c^_ x = (m — l)!/(A;!(m — k — 1)!). The 
functions u m {z) are defined inside the region \z\ < 1 and have a form 

00 

U m (z) = l + J2 ZnU< n l) ( ( l)- ( 6 - 7 ) 
n=l 



1 The monodromy matrix around z — is denned as Q n (ze 2 ") = M n kQ k (z) and similarly for the other 
singular points. 
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Inserting (|6.6|) and ()6.7|) into ()6.2ft . one derives recurrence relations for Un \q)- However, in 
order to save space, we do not show here their explicit form. 

In the anti-holomorphic sector the fundamental set of solutions can be obtained from (I6.6j) 
by substituting s and q k by s = 1 — s* and q k = q%, respectively. Sewing the two sectors we 
obtain the general solution for Q(z,z) around z = as 

TV 

Q(z,z) ]z] =° £ QS(z)C^QS(z). (6.8) 

m,m=l 

The above solution ()6.8|) should be single-valued on the z— plane. Thus, similarly to (|5.23j) we 
find a structure of the mixing matrix which for n + m < N + 1 

TV— n— m+1 ^ 2^fc 

C ™ = (n-l)Km-l)! ^ (6.9) 

v 7 v 7 fc=0 

with cr, ai, . . . , O/v-i being arbitrary complex parameters and cxn = 1- Below the main anti- 
diagonal, that is .for n + m > N + 1, Cnm vanish. 

The mixing matrix depends on N arbitrary complex parameters a and a k . However, 
two parity relations, Eqs. (I3.21jl and (|6.4j) . fix a = exp(i6*jv(g, q)), with (9jv(g, g) being the 
quasimomentum, and lead to the quantization of the quasimomentum. Later, we will use (16 .4j) 
to calculate the eigenvalues of 9]y(q,q) (see Eq. (|6.24ft ). 

The leading asymptotic behaviour of Q(z,z) for z — > can be obtained by substituting 
jEIJ) and ((ESI) into (jHHj). It has a form 



Q q ,q{z,z) = z 1 -^ 1 - 8 e i0N{q ' q) 



Log N -\zz) hog N -\zz) 
{N-l)\ + (N — 2)! aN ~ V 

o 

~1! 



+ L ° g [ Zz) o 2 n, (l + 0(z,z)). (6.10) 



Making use of the integral identity 

r (_-\\n m \ 

+ £>((m - s - in) ) 



rf 2 ^ - r (-l) n n! 



z\<p ZZ 



m — s — lu) 

(6.11) 

with m and m positive integer, we can calculate the contribution of the small— z region to the 
eigenvalue of the Baxter equation (j6.ip . The function Q q ^(u,u) has poles of the order N in the 
points u = i(s — m) and u = i(s — m) what agrees with 1)3. 12|) . For m = m = 1 one finds from 

(jnnD 



p i8 N (q,q) 



Q q ,qi u t + e, u£ + e) = ,. w [1 + ie c^-i + ■ ■ ■ + {ie) N ~ 2 a 2 + (ie)^ 1 a x + 0(e 

(6.12) 

where uf and uf are defined in (13.12)1 . One can see that the integration in (J6.ll) over the 
region of large z with ()6.4|) and ()6.1fl|) gives the second set of poles for Q q ^{u, u) located at 
u = —i(s — m) and u = —i(s — m). 
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Comparing (|6.12|) with one obtains 



7T 



R + (q,q) = -^e i9 »^ ; E+fa-Q = aN ^q) . (6.13) 



i 



N 



Now, we may derive expression for the energy 

E N (q,q) = Re[atjv_i(-g,-g) +q;jv-i(?, q)] • (6.14) 

The arbitrary complex parameters a n , defined in (|6.9j) . will be fixed by the quantization con- 
ditions below. 

In this Section we have obtained following Ref. [22] the expression for the energy spectrum 
E]\f(q,q), as a function of the matrix elements of the mixing matrix Ij6.9j) in the fundamen- 
tal basis (J6.6I) . Moreover, we have defined the solution to the Baxter equation Q(u, u) and 
reproduced the analytical properties of the eigenvalues of the Baxter operator, Eq. (|3.12j) . 

6.2 Solution around z = 1 

Looking for a solution of (I6.2JI around z = 1 in a form Q(z) ~ (z — l) b we obtain the following 
indicial equation 

N-3 

{b + 1 + h - Ns){b + 2 - h - Ns) Yl (6 - *:) = , (6.15) 

fc=0 

where ft, is the total SL(2,C) spin defined in ()5.2j) . Although the solutions b = k with = 
0, . . . , JV — 3 differ from each other by an integer, for h (1 + Uh)/2, no logarithmic terms 
appear. The Log(z)— terms are only needed for Imh = where the additional degeneration 
occurs. 

Thus, we define the fundamental set of solutions to Eq. 1)6. 2j) around z = 1. For Imh ^ it 
has the form 

Q^\z) = z 1 - s (l-z) Ns - h - 1 v 1 (z), 
Q i £\z) = z 1 - s {l-z) Ns+h - 2 v 2 {z) ) 

QW(z) = z 1 - s (l-z) m - 3 v m (z), (6.16) 
with m = 3, . . . , N. The functions Vi(z) (i = 1, 2) and v m (z) given by the power series 

oo oo 

«i(z) = i+x;(i-«r«i (ff), «™w = i+ x; (i-*r^ m) (?), (6.17) 

n=l ra=-/V— m+1 

which converge inside the region |1 — z\ < 1 and where the expansion coefficients Vn and 
satisfy the N— term recurrence relations 2 with respect to the index n. For h — (1 + n h )/2 G 
2Z + 1, one Log(z)— terms appears so for > 0, 



= ^"'(l - 2 )^-K+3)/2 [(! _ 3)"hLog(l - ^) + t5i(*)] , (6.18) 

fe=(l+n h )/2 



2 The factor z 1 s was included in the r.h.s. of Ijfi.lfifl and llfi,18Jl to simplify the form of the recurrence relations. 
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where the function v 2 (z) is the same as before, V\(z) = J^'kLo VkZ h and the coefficients v k satisfy 
the N— term recurrence relations with the boundary condition v nh = 1. For /ifZwe have two 
additional terms: Log(z) and Log 2 (z). 

Similar calculations have to be performed in the anti-holomorphic sector with s and h 
replaced by s = 1 — s* and h = 1 — h*, respectively. A general solution for Q(z,z) for Im(/i) 7^ 
with respect to the single-valuedness can be constructed as 

N 

Q(z,z) le ^ 1 ^Q?(z)Qf\z) + /3 1 ^ h Q?(z)Qf\z) + £ Q% (z) lmm Q%> (z) ■ (6-19) 

m,m=3 

Here the parameters /3h and ^mm build the C^ 1 ) matrix where Q(z,z) = Qmd^Q^-. The 
(3— coefficients depend, in general, on the total spin h (and h = 1 — h*). They are chosen 
in (I6.19J) in such a way that the symmetry of the eigenvalues of the Baxter operator under 
h — > 1 — h becomes manifest. Thus, the mixing matrix C^ x > defined in (I6.19j) depends on 
2 + (N — 2) 2 complex parameters flh, ft\-h and ^ m m which are some functions of the integrals 
of motion (q,q), so, they can be fixed by the quantization conditions. 

For h — (1 + nh)j2 the first two terms in the r.h.s. of (|6.19l) look differently in virtue of 

dnzHD 



Q{z,z) 



h=(l+n h )/2 



Q < t\z)Q^\z)+Q < i\z)Qf{z)\ + p 2 Q { 2 1 \z)Q i 2 \z) + ..., (6.20) 



where ellipses denote the remaining terms. Substituting (I6.19j) into (I6.1j) and performing in- 
tegration over the region of |1 — z\ <C 1, one can find the asymptotic behaviour of Q{u,u) at 
large u. 

Let us consider the duality relation ()6.4|) . Using the function Q(z, z) we evaluate ()6.19|) in 
the limit \z\ — > 1. In this way, we obtain set of relations for the functions Pi(q,q) and •y m m(q,q)- 
The derivation is based on the following property 

N 

(1/z; -q) = Sat Qi 1} (z; g) , (6.21) 

b=l 

with Im(l/,2) > and where the dependence on the integrals of motion was explicitly indicated. 
Here taking limit z — » 1 in (|6.16|) and (|6.17|) and substituting them to (I6.21|) we are able to 
evaluate the S— matrix 

S U = e-^ N *- h -V , S 22 = e -^ N '+ h -V , S k , k+m = (-1) ^ ~ 2 " ~ 1)m (6.22) 

ml 

with (x) m = T(x + m)/r(x), 3 < k < N and < m < N — k. Similar relations hold in the 
anti-holomorphic sector, 

£11 = e-^" 1 ) , S 22 = e MJm*-» , S k , k+m = (_i)*-3 (*- 2g - 1 )"> . (6 .23) 

m! 

The S 1 — matrix does not depend on z because the Q— functions on the both sides of relation 
()6.2H) satisfy the same differential equation (|6.2|) . 
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Now, substituting ([6~TT?|) and (l6~2l^l into (jOjl . we find 

N 

lmm(q,q) = e^,f) S n mlm(-q, -q)Sym ■ (6.24) 

n,n>3 

In this way, similarly to the energy, Eq. (|6.14j) . which was calculated from the mixing matrix 
at z = 0, the eigenvalues of the quasimomentum, ^(g, g), maybe calculated from the mixing 
matrix at z — 1, from the first relation in ()6.24|) . In the special case when q2k+i = <?2fc+i — 
(A; = 1,2.. .), this means fih(q,q) — Ph(—q, —q), the quasimomentum is equal to 

e i8 N (q,q) _ ^_^Nn s +n h 

6.3 Transition matrices 

In the previous Sections we constructed the solutions Q(z, z) to ()6.2|) in the vicinity of z = 
and z — 1. Now, we sew these solutions inside the region |1 — z\ < 1, \z\ < 1 and, then 
analytically continue the resulting expression for Q(z, z) into the whole complex z— plane by 
making use of the duality relation ()6.4|) . 

The sewing procedure is similar to that described by (|5.37p - (|5.4ip . Firstly, we define the 
transition matrices fl(q) and tt(q): 

N N 

Q$>(*) = J2 Q ^(q)Q^(z), Q ( °\z) = £n nni (g)Q^). (6.26) 

m=l m=l 

which are uniquely fixed (|5.40[) . The resulting expressions for the matrices f2(g) and fl(q) take 
the form of infinite series in q and q, respectively. Substituting (|6.26|) into ()6.8|) and matching 
the result into (I6.19JI , we find the following relation 

C^(q,q) = [n(q)] T C^(q,q)m. (6.27) 

The above matrix equation allows us to determine the matrices and and provides the 
quantization conditions for the integrals of motion, q k and q k with k = 3, . . . , N. Therefore, we 
can evaluate the eigenvalues of the Baxter Q— operator, Eq. (|6.1|) . Formula (|6.27ft contains N 2 
equations with: 

• (N — 1) a— parameters inside the matrix 

• 2 + (N — 2) 2 parameters /?i,2 and 'j mm inside the matrix C^, 

• (N — 2) integrals of motion q 3 , . . . , q^ where q k = ql 

Thus, we obtain (2N — 3) nontrivial consistency conditions. 

The solutions to the quantization conditions (|6.27|) will be presented in details in next 
Sections. 
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Chapter 7 



Numerical results 



In this Chapter we present the spectra of the conformal charges obtained by numerical calcula- 
tions [2S1 ED]- To this end we resum solutions (lfi.7j) and (|6.17j) numerically, and using them we 
solve the quantization conditions (|6.27jl . First we discuss the numerical results obtained in Refs. 
[50l 123] for N = 3, 4 and for the ground states for N = 5, . . . , 8. Moreover, we present some 
results which were never published before, i.e. quantized values of g 3 for N = 3 and > 0, 
resemblant and winding spectra of q^, q^ for N = 4 and corrections to the WKB approximation 
for N = 3 and N = 4. 

7.0.1 Trajectories 

Solving quantization conditions (|6.27|) we obtain continuous trajectories in the space of confor- 
mal charges. They are built of points, (q 2 (vh), ■ ■ • , QN^h)) which satisfy (|fj.27jl and depend on 
a continuous real parameter Uh entering q 2 , (15. 2 J) and (I2.fi()j) . In order to label the trajectories 
we introduce the set of the integers 

£ = {£ 1} £ 2} ...,e 2(N - 2) } (7.1) 

which parameterize one specified point on each trajectory for given h. Specific examples in the 
following Sections will further clarify this point. 

Next we calculate the observables along these trajectories, namely the energy (I6.14JI and 
the quasimomentum (|6.24|) . The quasimomentum is constant (I2.70B for all points situated on a 
given trajectory. The minimum of the energy, which means the maximal intercept, for almost 
all trajectories is located at u h = 0. It turns out that the energy behaves around v h = like 

E N (u h ; F round ) = £fT nd + <W + 0[v h 2 ) (7.2) 

Thus, the ground state along its trajectory is gapless and the leading contribution to the 
scattering amplitude around Uh may be rewritten as a series in the strong coupling constant: 

00 s -a B E% ound /4 

A(s,t) ~ -is y2(m s ) N — -——^(t)( B ^(t) , (7.3) 
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where a s = a s N c /n and £x,N(t) are the impact factors corresponding to the overlap between 
the wave-functions of scattered particle with the wave-function of JV— Reggeons, whereas a N 
measures the dispersion of the energy on the the trajectory around = 0. 

On the other hand, the energy along the trajectories grows with Uh and for \uh\ — > oo and 
finally, we have E^iy^t) ~ lni^ 2 . These parts of the trajectory give the lowest contribution 
to the scattering amplitude. 



7.0.2 Symmetries 

The spectrum of quantized charges qz, ■ ■ ■ ,Qn is degenerate. This degeneration is caused by 
two symmetries: 

q k <-> {-l) k q k (7.4) 

which comes from invariance of the Hamiltonian under mirror permutations of particles, (I2.68|) . 
and 

qk <-> q k (7-5) 

which is connected with the symmetry under interchange of the z— and 1— sectors. Therefore, 
the four points, {qk}, {(— l) fe 9fe}j {qk*} and {(— l) k qk*} with k = 2, . . . , JV, are related and all 
of them satisfy the quantization conditions 1)6. 27J) and have the same energy. 



7.0.3 Descendent states 

Let us first discuss the spectrum along the trajectories with the highest conformal charge q^ 
equal zero for arbitrary u h G IR. It turns out [E3 IHU E-l E3 that the wave-functions of these 
states are built of (JV — 1)— particle states. Moreover, their energies are also equal to the 
energy of the ancestor (N — 1)— particle states: 

E N {q2, q 3 ,...,q N =0) = E N -i(q 2 , qN-i) ■ (7.6) 

Thus, we call them the descendent states of the (N — 1)— particle states. 

Generally, for odd N, the descendent state ^/^ N ~ ^ with the minimal energy E^iqN = 0) = 
has for q 2 = 0, i.e. for h — 0, 1, the remaining integrals of motion g 3 = . . . = q^ = as well. 
For h — 1 + ii>h, i-e. q 2 ^ 0, the odd conformal charges q 2 k+i = with k — 1, . . . , (N — l)/2 
while the even ones q 2 k ^ and depend on v^. 

On the other hand, for even jV, the eigenstate with the minimal energy ^l/^ -0 ^ is the 
descendent state of the (JV — 1)— particle state which has minimal energy with q^-i 7^ 0. Thus, 
E™ n (q N = 0) = E^iq^ ± 0) > 0. 

Studying more exactly this problem one can obtain (HSE!] a relation between the quasimo- 
menum 9n of the descendent state and the ancestor one 0/v-ij which takes the following form 



(-l)^ 1 . (7.7) 

<?jv=0 
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Additionally, one can define linear operator A [HEl EH] that maps the subspace Vjf^ 
of the (N — 1)— particle ancestor eigenstates with the quasimomentum O^-i = t^N into the 
iV— particle descendent states with q^ = and On = ir(N + 1) as 

A : y$ff 1=7rJV) _> V ^=AN+D) (7 8) 

It turns out that this operator is nilpotent for the eigenstates which form trajectories i-e. 
A 2 ^/ = 0. Thus, the descendent-state trajectory can not be ancestor one for (N + 1)— particle 
states. However, it is possible to built a single state [23j with g 2 = q 3 = . . . = q^ = 0, i.e. for 
only one point Uh = 0, that has En = and the eigenvalue of Baxter Q— operator defined as 

grW)-^, (7.9) 

u 

where a normalization factor was omitted. For N = 3 this state corresponds to the wave- 
function defined in (|5.51|) . 

Additional examples of the descendent states for iV = 3 and iV = 4 will be described later 
in the next Sections. 



7.1 Quantum numbers of the = 3 states 

In this section we present the spectrum of q% for three reggeized gluons. For the first time such 
solutions for Re[g 3 ] = were obtained in [48j. The authors of Ref. [48j used the method of the 
^3 eigenfunctions described in Chapter 5. Similar quantization condition was also constructed 
in [HI]. Solutions with Refgs] ^ were found in IIHIEDI- Moreover, in the latter paper the 
solutions with rih ^ are described. 

The first solution using the Baxter Q— operator method was described in Ref. [23]. Later, 
similar results were also obtained in Ref. |2Hj- R turns out that the Baxter Q— operator method 
[2B] for N = 3 is equivalent to the method of the q 3 eigenfunctions described in Chapter 5. For 
higher N > 3 only the Baxter Q— operator method was used to find the quantization values of 
gjv and the energy. 



7.1.1 Lattice structure 

Solving the quantization conditions (|6.27j) for N = 3 and for h = ^-y^, i.e. with v h = 0, we 

1/3 

reconstruct the full spectrum of (fc. It is convenient to show the spectrum in terms of q 3 rather 

1 /3 

than q 3 . Since q 3 is a multi- valued function of complex q s , each eigenstate is represented on 

1/3 

the complex q 3 —plane by N = 3 different points. Thus, the spectrum is symmetric under the 
transformation 

ql^ 3 <-> exp(27ri/c/3)g3^ 3 , where < k < 3 . (7-10) 

Additionally, symmetry (|7.4p gives a more regular structure. For the total SL(2, C) spin of 
the system h = 1/2, which means rih = and u h = 0, we present the spectrum in Figure 
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Figure 7.1: The spectrum of quantized q 3 for the system of N = 3 particles. On the left the 
total SL(2, C) spin of the system is equal to h = \, while on the right h = 1. Different symbols 
stand for different quasimomenta #3: stars 63 = boxes #3 = Air/ 3 triangles 63 = 2ti/3. 

17.11 One can easily notice that the spectrum has the structure close to the equilateral triangle 

lattice of the leading order WKB approximation ()4.2H) . Indeed, apart from a few points close 

1/3 

to the origin, the quantized values of g 3 are located almost exactly at the vertices of the WKB 
lattice. The WKB formula (23] gives 



(7.11) 



where £\ and £\ are integers, such that their sum £\ + £\ is even. Here the lattice spacing is 
denoted by 

1 r 3 (2/3) 





f 1 dx 


_ 4 l/3 7r J 


-00 \/l - x 3 _ 



2tt 



0.395175... . 



(7.12) 



The lattice of q^ 3 extends on the whole complex plane except the interior of the disk with the 
radius A 3 : 

l<7 3 1/3 | > A 3 (7.13) 

situated at g 3 = 0. 

In accordance with t)7. 1 1^ . a pair of integers £\ and £2 parameterize the quantized values of 
q^ 3 . Going further, one can calculate the quasimomentum as a function of £\ and £2. It has a 
following form 

9 3 (£i,£2) = y £i ( mod27 0- ( 7 - 14 ) 

Thus, as we can see states with the same value of Re[gg^ 3 ] have the same quasimomentum. In 
Figure 17.11 different quasimomenta are distinguished by stars, boxes and triangles. 
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Figure 7.2: The spectrum of quantized q 3 for the system of N = 3 particles. On the left the 
total SL(2, C) spin of the system is equal to h — §, while on the right h = 2. Different symbols 
stand for different quasimomenta #3: stars 63 = Torres #3 = Air/ 3 triangles 63 = 2ti/3. 
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The same lattice structure is exhibited by other spectra with different n h . However, they 

1/3 

have different corrections to the leading order WKB approximation for q 3 . These spectra 
are presented in Figures The corrections to the lattice structure depend on q 2 as seen 

in (|4.21|) . Since the WKB lattice is obtained in the leading order of the expansion for large 
conformal charges, 1 \qH 2 \ ^ IqI^Ij the corrections are bigger for lower \q^ 3 \- Later, we 
shall discuss some other features of the corrections to the WKB leading order approximation. 

As we can see in Figs. 17.1117.21 for h 6 Z we have additionally trajectories with q 3 = 0. 
They are called the descendent states because their spectra are related to the spectra for the 
N — 1 = 2 Reggeon states. We discuss this point further in next Sections. 



7.1.2 Trajectories in 

In the previous Section we considered the dependence of g 3 on n h for v h = 0. However, the 
spectrum of conformal charges also depends on the continuous parameter v h with h = ^-y^+iz/^. 
It turns out that the spectrum is built of trajectories parameterized by real parameter Each 
trajectory crosses one point (star, box, triangle) in Figs. l7.lH7.21 An example of three trajectories 
is presented in Figure l7~3l They are numbered by (£1, £ 2 ) = (0, 2), (2, 2) and (4, 2) whereas they 
quasimomentum #3(^1,^2) = 0, 47r/3 and 27r/2, respectively. 

1/3 

The trajectories cumulate at vy L = 0. When we increase Vh, q 3 tends to infinity and the 
structure of quantized charges starts to be less regular, especially for trajectories with lower 
jg^ 3 ]. We can see this in Figures [7~H where we project trajectories with h = ~ + ivh on the 
fh — plane. Here stars denote point with = 0, boxes with vy l = \ and circles = 2. Grey 
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Figure 7.3: The dependence of quantized q 3 {vh] £x, £2) on the total spin h = 1/2 + iv^. Three 
curves correspond to the trajectories with (^1,^2) = (0,2) , (2,2) and (4,2). 

lines are drawn to show the projection of the trajectories for intermediate values of u^. 

For rih 7^ we notice that the spectra start to rotate with v h . In Figures [?~4l and IT751 we 
present trajectories with positive rih = 0,1,2,3. Due to the symmetry (17. 5J) . which means 
h — ► 1 — h* or rih — > — nh with g 3 — > q 3 , or equivalently, z/^ — > —v^ but (73 — > q 3 , the spectrum 
for the negative n h is the same as for the positive ones but it rotates in the opposite direction 
with Vh- 

Some of the results presented in this Section were found in earlier works [SHIED!- Trajectories 
with quasimomentum #3 = and nu = were obtained in Refs. [SHIED!- The case for h = 2+iUh 
with quasimomentum #3 = was discussed in Ref. [SDJ- In this thesis we additionally analyse 
the spectra for h — 1 + ivn and h = 3/2 + iv h . 

7.1.3 Energy and dispersion 

For all trajectories in (g 2 , g 3 )— space we can calculate the energy of the reggeized gluons using 
Eq. 1J6.14JI . Example of the energy spectrum for trajectories from Figure 1731 with h = | + Wh 
is shown in Figure I7TH1 

The energy along the trajectories is a continuous gapless function of v^. As we can see the 
energy E 3 grows with rising \i/h\. For = it has a minimum value mm Uh E 3 (i'h; £1, £ 2 ) at 
v h = 0. In the case n h 7^ 0, due the bending of the trajectories some minima of the energy 
are moved away from v h = [50j. However, the ground state corresponds to the point (s) on 

1 /3 t— 

the plane of q 3 (see Figure ETJ closest to the origin. For N = 3 the ground state is located 
on the (0, 2)— trajectory at Uh — and nh = with quasimomentum equal 9 3 = 0. Due to the 
symmetry ()7.4|) it is doubly-degenerated and its conformal charge and energy take the following 
values: 

^ground = ±Q 2Q526 _ _ _ ^ ^ground = Q. 98868 .... (7.15) 
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In the vicinity of v h = the accumulation of the energy levels is described by the dispersion 
parameter <t 3 ()7.2ft given below in the Table 17751 

We show comparison the WKB result of Eq. (17.11^ with the exact expressions for q 3 
at h = 1/2 in Figure 17711 and Table IT7T1 One can find that the expression fj7. 1 1 j> describes 
the excited eigenstates with good accuracy. In the case where the eigenstates have smaller 
q 3 agreement becomes less accurate. Thus, for the ground state with iq 3 = 0.20526 . . . the 
accuracy of (17.11^ is ~ 20%. Obviously, in the region where the WKB expansion is valid, 
i.e. \q^ 3 \ ^> \q^ 2 \, Eq. (|7.11j) can be systematically improved by including subleading WKB 
corrections. 





(q | xact ) 1/3 


(g WK B) l/3 


-E 3 /4 


(0,2) 


0.590 i 


0.684 i 


-0.2472 


(2,2) 


0.358 + 0.621 % 


0.395 + 0.684i 


-0.6910 


(4,2) 


0.749 + 0.649 % 


0.790 + 0.684 i 


-1.7080 


(6,2) 


1.150 + 0.664 i 


1.186 + 0.684 a 


-2.5847 


(8,2) 


1.551 + 0.672 i 


1.581 + 0.684i 


-3.3073 


(10,2) 


1.951 + 0.676 % 


1.976 + 0.684 i 


-3.9071 



Table 7.1: Comparison of the exact spectrum of q 3 at h = 1/2 with the approximate WKB 
expression (I7.11j) . The last line defines the corresponding energy E 3 (0;£i,£ 2 ). 
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Figure 7.5: The trajectories of q 3 projected on Uh = 0. On the left panel h = | + ivhi while 
on the right one h = 2 + iv^. stars denotes = 0, ftcres = 1 and triangles — 



7.1 .4 Descendent states for TV = 3 

One also can notice in Figures I7~4l and 17751 that for odd we have states with g 3 = 0. For 
n h = and = it has the energy E 3 = 0, so it is lower than (|7.15jl . These states are 
descendants of the states with two Reggeons [5T1 [5~jJ [HBJ EHJ. We constructed them using 
the qz eigenfunction method. The wave-functions of these states are described by (I5.44j) . (I5.46j) . 
(I5.51B and l)5.53jl . These states have the same properties and the energy as the corresponding 
states with N — 1 = 2 particles, E 3 (q 2 , (fe = 0) = E 2 (q 2 ), with [TBI fT6]: 

'1 + In J 



E 2 (q 2 )=4Re[ 



h) + ip(h) - 2^(1)] = 8 Re 



■0 



(7.16) 



where ijj(x) = 4-lnT(x) and q 2 = —h(h — 1). Moreover, their wave-functions are built of the 
two-Reggeon states [HBJEII] and the quasimomentum 9 3 = 0. Contrary to the states with g 3 ^ 0, 
the states with q 3 = (|5.46j) couple to a point-like hadronic impact factors [5311201111], like the 
one for the 7* — > 1] C transition. 



7.1.5 Corrections to WKB 

The WKB formula for the lattice structure of the conformal charge q 3 was derived in paper 
[2T| . This formula tells us that for g 3 — > 00 

2 00 / 7 \ fc" 

1 + 



1/3 £3(2/3) 



2tt 



\Q(n) 



Q(n) 



(7.17) 
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-2-10 1 2 



Figure 7.6: The energy spectrum corresponding to three trajectories shown in Figure I731 The 
ground state is located on the (0, 2)— trajectory at v h = 0. 



where 

°'" J -2- - 1 

k=l 

and li, l 2 , n = {rii, . . . ,n N } are integers, while the coefficient 



G(n) = \{h + k) + i^-(h - k) = Y,^e m{2k - 1)/3 (7.18) 



b = ^\ (7.19) 

where star denotes complex conjugation. 

After numerical calculations we have noticed that better agreement with the exact results 
is obtained for 

3V ~ 3 <*■ - I) • (T.20, 



IF V s " 3 

In order to show this, we calculated the values of the conformal charge q 3 for h = ^-y^ for 
rih = 0, 1, . . . , 19. We evaluated numerically g 3 with Im[g 3 ] = 0, i.e. and separately 

with Refgs] = 0, i.e. a^' where the superscript (i) and (r) refers to imaginary and real parts, 
respectively. Then we fitted expansion coefficients for large \q^\ in the range . . . 5000 
with high numerical precision. In order to save space in Table l?~2l we present results only for 
n h = 0, 1, 2 and 3. 

Coefficients for k — 0, . . . , 2 agree with formula (I7.17L i.e. do = 1, a\ — 1 and a 2 = — 1, 
but as we previously mentioned in ()7.19j) and ([7.20(1 . the expansion parameter b/ \Q(n)\ 2 is 
different. This difference comes from the fact that the series formula (|7.17|) with ()7.19|) is 
derived in the limit 1 \q^ 2 \ <S \<ll^ 3 \- Since in (|7.2()jl the value of q\ is much bigger then 2/3 
the value of the parameter b from (|7.20j) in the above limit goes to (|7.19jl . One may suppose 
that therefore the factor 2/3 in Eq. ()7.20|) as subleading was omitted in the derivation presented 
in Ref. EH- 
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22 76106772 

• 1 V7 -L V7 V7 1 1 *—> 


79 325249 

1 !_/ • '.7 .7 i— J- ' 7 


268.13 




W 

"ft 


1.749245141092 


3 2779793873 

'.J * 1 1 t_7 1 t-7 '.7 ^7 1 ".7 


5.23893228 


-0 430149 

V7 • X '.7 V7 _L J- »7 


-75.81 




a, + a, 

ft 1 ft 


4.000000000000 


10.0000000000 


28.00000000 


78.895100 


192.32 



1 /3 *IT^™^T5li 

Table 7.2: The fitted coefficient to the series formula of q 3 ()7.17|) with = 0, 1, 2 and 4 



Secondly, we see that the coefficients with k > 2 start to depend on n h . Thus, to describe 

1/3 

the behaviour of q 3 properly, we have to introduce a second expansion parameter, for example 
q2- We can also notice that for k > 2 for a few first coefficients a® + ajP G Z and this sum 
does not depend on n h . 

Moreover, we can see that after the numerical fitting we obtain two different sets of the ex- 
pansion coefficient, {a^} and {a^}, defined in ()7.17|) . for real and imaginary g^ 3 , respectively. 
Thus, in order to describe full-complex values of g^ 3 in terms of the series 1)7.17)1 we have to use 

1 /3 

both sets of coefficients, one for real and one for imaginary part of q 3 . Alternatively, we can 
perform expansion with two small parameters, i.e. <?2/|2( n )| 2 an d l/|2( n )| 2 - Since the leading 
terms, i.e. with a , a\ and 02, for real and imaginary q^ 3 are equal and known analytically, 
good approximation is obtained using Eq. 1)7.17)1 with ()7.20jl and neglecting higher order terms 
with a fe > 3 . 

7.2 Quantum numbers of the N = 4 states 

In this Section we present the spectrum for four Reggeons. Earlier, some results for iV = 4 
were only presented in (2E| and some numerical results in Ref. (2Hj- Here we show much more 
data and we present more detailed analysis of this spectrum. 

For four Reggeons the spectrum of the conformal charges is much more complicated than in 
the three-Reggeon case. Indeed we have here the space of three conformal charges (<?2, ?3, ^4)- 
Thus, apart from the lattice structure in q 1 ^ we have also respective lattice structures in 
^3— space. Here we consider the case for n h = so that h = ^ + iv h . This spectrum includes the 
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ground states. For clarity we split these spectra into several parts. We perform this separation 
by considering spectra with different quasimomenta 9±(q, q) as well as a different quantum 
number £ 3l which will be defined in the solution ()7.22|) - ()7.23|) of two quantization conditions 
(14481) for iV = 4 from Ref. (23] . 

From the first quantization condition (J4.18j) one can get the WKB approximation of the 
charge g 4 as 

■ •< r 2 (3/4) r l z j 

(7.21) 



1/4 
94 



4^ 



(mod 27r) 



(7.22) 



and the quasimomentum is equal to 

2 2 

where £i, £2 and £3 are even for even £ and odd for odd £. Thus, we have two kinds of lattices: 
with 6*4 = 0, ii and with #4 = ±7r/2. They are presented in Fig. 17.71 In these pictures gray lines 
show the WKB lattice (|7.21j) with vertices at £±,£2 £ Z. To find the leading approximation for 
the charge q 3 , we apply the second relation in 1)4. 18j) with m = 2 which gives 



Im 



1/2 

94 



(*1 " h 



(7.23) 



Notice that the system (|7.21|) and (j7.23j) is underdetermined and it does not fix the charge q 3 
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Figure 7.7: The spectrum of the integrals of motion g 4 for iV = 4 and the total spin h = 1/2. 
The left and right panels correspond to the eigenstates with different quasimomenta e 104 = ±1 
and ±z, respectively. 

It turns out that after choosing one value of #4, the lattice in q 1 / 4 — space is still spuriously 
degenerated 1 and this degeneration also corresponds to different lattices in g 1 / 2 . The parameter 
£3 which is defined in ()7.23|) will be used to distinguish these different lattices. 
1 degeneration in the leading order of the WKB approximation 



69 



7.2.1 Descendent states for N = 4 



One can notice that for N = 4 and n h = we have the descendent states. They appear in 
sector with the quasimomentum # 4 = tt, which agree with (|7.7jl . The wave-functions of the 
descendent states are built of three-particle eigenstates with 6*3 = 0. These three-Reggeon 
states correspond to the lattice depicted by stars in Fig. 17.11 On the other hand, the lattice 
for the descendent state for N = 4 is presented in Fig. 17.81 on the left panel. Thus, one can 
compare the both lattices and notice that they are exactly the same. Moreover, the energy of 
these descendent state and the corresponding three-Reggeon states are also the same (|7.6j) . 




-1.5 -1 -0.5 0.5 1 1.5 -1.5 -1 -0.5 0.5 1 1.5 



Re[g 3 1/3 ] Re[g 4 1/4 ] 

Figure 7.8: The spectra of the conformal charges for N = 4 and comparison to the WKB 
expansion. On the left panel, the spectrum of g 3 with g 4 = corresponding to the descendent 
states with 6*4 = tt. On the right panel, the spectrum of g 4 for h = 1/2 and = with 64 = 0. 
The WKB lattices are denoted by the gray lines. 



7.2.2 Lattice structure for q% = 

Let us consider the spectrum with g 4 7^ and q 3 = 0, see the right panel of Fig 17.81 In this 
case the quasimomentum 64 = and the lattice structure include vertices that correspond to 
the ground state. 

Similarly to the N = 3 case we have in the g]^ 4 — space a lattice with a square-like structure 
described by (|7.21|) . In this case even numbers t\ and £ 2 satisfy £ 4 + £ 2 G 4Z. Thus, we have 
the WKB formula 

[nM y^^{^^), (7.24, 
where the vertices are placed outside a disk around the origin of the radius 

0.423606.... (7.25) 
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"43/4 


r 1 dx 


1 r 2 (3/4) 


U j 


_i Vi -x A _ 


2v^F 



As before, the leading-order WKB formula ()7.24|) is valid only for \q l J 4 \ 3> \q\ 




-4 -2 2 4 -4 -2 2 4 



Figure 7.9: The dependence of the energy, —E4/A, and the quantum number, (74/(72, with 
g 2 = 1/4 + 1/^ j on the total spin h = 1/2 + iuh along the ground state trajectory for N = 4. 

The energy is lower for points which are nearer to the origin. Similarly to the N = 3 case, 
the spectrum is also built of trajectories which extend in the [y h , q 3 , q^j — space. Namely, each 
point on the g^ 4 — lattice belongs to one specific trajectory parameterized by the set of integers 

{£i,e 2 ,...}. 

The ground state for N = 4 is situated on a trajectory with (£1,^2) = (4,0). We find that 
for this trajectory g 3 = Img 4 = 0, whereas Re[g 4 ] and £4 vary with v h as we show in Figure I7~9l 
An accumulation of the energy levels in the vicinity of v h = is described by Eq. (|7.2|) with 
the dispersion parameter 04 given below in Table 1731 

On the q^ 4 — plane the ground state is represented by four points with the coordinates 
^2) = (±4, 0) and (0, ±4). Due to a residual symmetry g 1 / 4 <-> exp(ikTv/2)q\' 4: , they describe 
a single eigenstate with 

g f ound = , gf ound = 0.153589 . . . , Ef ound = -2.696640 .... (7.26) 

with h = 1/2. It has the quasimomentum 6*4 = and, in contrast to the N = 3 case, it is 
unique. 

The comparison of (|7.24j) with the exact results for 54 at h = 1/2 is shown in Figure l7~8l 
and Table 17731 One can see that the WKB formula (I7.24J) describes the spectrum with a good 
accuracy. 

7.2.3 Resemblant lattices with ^3 = 

In the previous Section we introduced the parameter £ 3 which helps us to distinguish different 
lattices. This parameter takes even values for 6*4 = 0, 7r and odd ones for # 4 = ±n/2. Let us 

1 /2 

take £ 3 = 0. It turns out that in this case the spectrum lattice in the q 3 —space is similar 
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(£ l /2,e 2 /2) 


/ exact) 1/4 


( g WKB)V4 


-EJA 

^ I 


(2,0) 


0.626 


0.599 


0.6742 


(2,2) 


0.520 + 0.520 i 


O.599 + 0.599 i 


-1.3783 


(3,1) 


0.847 + 0.268 % 


0.899 + 0.299 i 


-1.7919 


(4,0) 


1.158 


1.198 


-2.8356 


(3,3) 


0.860 + 0.860i 


0.899 + 0.899 i 


-3.1410 


(4,2) 


1.159 + 0.574 i 


1.198 + 0.599 i 


-3.3487 



Table 7.3: Comparison of the exact spectrum of g 4 at g 3 = and h = 1/2 with the approximate 
WKB expression 17.241 The last column shows the exact energy E4. 

to the corresponding lattice in the q l J A — space, i.e. considering only the leading order of the 
WKB approximation the g^ 2 — lattice in comparison to the g]^ 4 — lattice is rescaled by some real 
number. An example of such a lattice for 6*4 = is shown in Figure 17.101 One can notice 
that the non-leading corrections to the WKB approximation cause the bending of the lattice 
structure in Fig. 17.101 for g]/ 4 concave whereas for g 1 / 2 convex. Moreover, we can see the 
g]^ 4 — lattice as well as g^ 2 — one have a similar structure to the lattice with g 3 = presented in 
Fig. 17.81 These lattices also do not have the vertices inside the disk at the origin g 3 = g 4 = 0. 
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Figure 7.10: The spectra of the conformal charges for iV = 4 with 6*4 = 0, £ 3 = and £ 4 = 1. 
On the left panel the spectrum of g^ 4 , while on the right panel the spectrum of g 1 / 2 



Substituting 

g i/ 2 = r 3 e^ 3 and g] /4 = r 4 e^ 4 (7.27) 
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Figure 7.11: The spectra of the conformal charges for N = 4 with #4 = 0, £3 = and £4 = 2. 
On the left panel the spectrum of q l J A while on the right panel the spectrum of q\ 



into (|7.23|) we obtain a condition for the leading order of the WKB approximation 

2 



sin (2(03 - 4 )) . 



(7.28) 



Thus, for £3 = and for a scale A = r 3 /r 4 > we have 3 = 4 . This means that the vertices 
on the qH" 2 — lattice have the same angular coordinates as those from the q l J 4 — lattice. Looking 
at the numerical results in Figure 17.101 we notice that the missing quantization condition for 
£3 = in the leading WKB order should have a form similar to 

?3 



Re 



1/2 
1a 



A/ 



(7.29) 



where A^ 4 = r 3 /r 4 G 1 is a constant scale for a given lattice £ 4 . 

It turns out that for a specified quasimomentum we have an infinite number of such lattices. 
They differ from each other by a given scale Xe 4 . For example for 6*4 = we have another 
lattice, shown in Fig 17.111 Its scale A2 differs from the scale Ai of the lattice from Fig. 17.101 
The resemblant ql^ 4 — lattices for #4 = are described by (|7.21|) with the integer parameters t\ 
and £ 2 satisfying t\ + £ 2 £ 4Z. 

The similar lattices also exist in the sector with the quasimomentum #4 = n. Some of 
them are presented in Figures 17.121 and 17.131 For 6*4 = tt the resemblant g] 7 ' 4 — lattices are also 
described by (|7.2H) but the integer parameters t\ and £ 2 satisfy t\ + £ 2 G 4Z + 2. 



7.2.4 Winding lattices with £ 3 ^0 

In the case with £3 7^ we have much more complicated situation than for £3 = 0. According 
to (|7.28j) the angles 3 and 4 defined in (|7.27|) are no more equal. Moreover, they start to 
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Figure 7.12: The spectra of the conformal charges for A^ = 4 with # 4 = ir, £3 = and £4 = 1. 
On the left panel the spectrum of q l J A while on the right panel the spectrum of q^ 2 
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Figure 7.13: The spectra of the conformal charges for iV = 4 with 6* 4 = ir, £ 3 = and £4 = 2. 
On the left panel the spectrum of q l J A \ while on the right panel the spectrum of q l J 2 



depend on the scale A = r 3 /r 4 . 

An example of such a lattice, with £ 3 = 1 and #4 = — n/2 we show in Figure EHH For 
this case, the q l J 4 — lattice is defined by (|7.21|) with £\ and £ 2 odd integer numbers satisfying 

^ ^ 1 /2 

£\ + £2 £ 4Z. In Figure I7.14| in order to present the correspondence between the q 3 — and 
q l J 4 — lattice, we depict only some vertices of the lattice which extends in the whole plane of 
the conformal charges except the place nearby the origin g 4 = q 3 = 0. As we can see the 
ql^ 2 — lattice is still a square-like one but it winds around the origin qlj 2 = 0. Looking at Figure 
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Figure 7.14: The winding spectrum of the conformal charges for N = 4 and h = 1/2 with 
6*4 = — 7i / 2 and £ 3 = 1. On the left panel the spectrum of q l J A , while on the right panel the 
spectrum of q]J 2 

I7.14[ let us start from 3 = and 4 = — 7r/4 where 3 and 4 are defined in (|7.27jl . Thus, 
the difference 3 — 04 = 7r/4 so that our scale at the beginning is A = \/2. In this region 

1/2 

the vertices of the g 3 —lattice are in the nearest place to the origin. When we go clockwise 
around the origin of the lattices by decreasing 4 , we notice that 3 also decreases but much 
slower. Thus, according to l)7.28j) the difference 3 — 4 changes. Moreover, due to (|7.28p . the 

1/2 

scale A also continuously rises. Therefore, the g 3 —lattice winds in a different way than the 
g]^ 4 — lattice. After one revolution the vertices of the g]^ 4 — lattices are at similar places as those 
with 4 decreased by 2ir. This provides additional spurious degeneration in g 1 / 4 . However, 
after revolution by the angle 27r, the vertices in g 1 / 2 — space have completely different conformal 
charges g 3 . The spurious degeneration in the g^ 2 — lattice does not appear. 

We have to add that we obtain the second symmetric structure when we go in the opposite 
direction, i.e. anti-clockwise. Moreover, the winding lattice, like all other lattices, extends 
to infinity on the g^ 2 — and g]^ 4 — plane and does not have vertices in vicinity of the origin, 
1z = 1i = 0. However, the radii of these empty spaces grow with 3 4 . 

Additionally, due to symmetry of the spectrum (|7.5|) we have a twin lattice with g& — > gj*. 
Furthermore, the second symmetry (17.41) produces another lattice rotated in the g 3 ^ 2 — space by 
an angle tt/2. Notice that this symmetry (17.4ft exchanges quasimomentum 84 <-> 2tt — 64 (mod 2tt) 
Thus, the spectrum for #4 = n/2 is congruent with the spectrum of # 4 = — 7r/2 but it is rotated 
in the g 3 ^ 2 — space by tt/2. 

As we said before for 6*4 ± 7i/2 we have winding spectra with odd £3. Similarly, for even 
4 ^ we have winding spectra with # 4 = 0, tt, thus, the lattice with the lowest non-zero £3 
corresponds to \£ 3 \ = 2. We present some points of this spectrum in Fig. \7.15\ for which £\, £ 2 
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Re[g] /4 ] Re^ j 

Figure 7.15: The winding spectrum of the conformal charges for N = 4 with h = 1/2, 84 = 
and £3 = 2. On the left panel the spectrum of q l J 4 while on the right panel the spectrum of qj 2 
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l/2l 
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15 



in (I7.21B are even and £\ + £2 G 4Z + 2. In this case we start with 3 = 7r/4 and 4 = which 
implies the beginning scale A = \J2. The spectra wind as in the previous case. 

Similarly, for 64 = 11 we have also spectrum with the lowest \£ 3 \ = 2. It is defined by 
(|7.21j) with £1, £ 2 even and £\ + £ 2 G 4Z. In this case the angles 3 = and 4 = — 7r/4 so 
we also have the beginning scale A = \/2, as depicted in Figure l7~T6l In order to describe the 
winding spectra better we may introduce an integer parameter £4 which helps us to number 
the overlapping winding planes of the spectra and name spuriously-degenerated vertices in the 
q]f 4 — plane. 

To sum up, even for a given quasimomentum we have many lattices which overlap, so 

1/2 

that vertices of the lattices, especially in q 3 —space, make an impression of being randomly 
distributed. However, as we have shown above, those spectra may be distinguished and finally 
described by (I7.21J) and (I7.23J) . Still, there is a lack of one nontrivial WKB condition which 
would uniquely explain the structure of the resemblant and winding lattices. 



7.2.5 Corrections to WKB 

Let us consider the spectrum of the conformal charge g 4 for N = 4 with q 3 = and h = ^-y^. 
It turns out that for n h 7^ it has similar square-like lattice structure like that with n h = 0, see 
Fig. l7.8l on the right panel. Similarly to the case with three reggeized gluons we have evaluated 
the conformal charges q^ with even rih with high precision. We have done it separately for g 4 
with Im[g 4 ] = and Re[g 4 ] = 0. Next, we have fitted coefficient expansion of the WKB series, 
aj^ and respectively. 
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Figure 7.16: The winding spectrum of the conformal charges for N = 4 with h = 1/2, 64 = 71 
and £3 = 2. On the left panel the spectrum of q l J 4 while on the right panel the spectrum of q^ 2 



In [21] the series formula for q)/^ looks as follows 

^3/2 



1/4 
Q4 



2r 2 (i/4) 



Q(n) 



b - f b 



(7.30) 



where 



Q(n) = ^ n 



k=l 



Mr(2k-l)/4 = ( Z l ! „• ^2 

\/2 



(7.31) 



and £1/2, ^ = {fti, ■ • • , ^tv} are integer. 

Here for iV = 4, similarly to the N = 3 case (|7.20ft . we have a different expansion parameter 
bj |Q(n)| 2 and in this case the parameter 



h 4 ( * 5 

7T \ 4 



(7.32) 



is decreased by 5/4. The expansion coefficients of the series (|7.3()j) are shown in Table [7~H 
Similarly to [21] the coefficients do = 1 and Oi = 1. The remaining coefficients depend on rih- 
Moreover, the coefficients with k > 1 are different for real g^ 4 and for imaginary q\^ A but 
one may notice that for k — 1, 2 the sum + G Z and does not depend on n h . Thus, to 
describe the quantized values of q]/ 4 more generally one has to use both sets of the coefficients, 



(r) 



and a k * , or perform the expansion with two small independent parameters, i.e. q2/\Q( n ) 



and l/\Q(n)\ 2 . 

Using the series (I7.30J) with (|7.32jl and coefficients from Table EU gives good approximation 
of the conformal charges g 4 with g 3 = 0. However, if someone wants to have a better precision 
one has to introduce an additional expansion parameter. 
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coef. 


k — 2 


fc = 3 


jfe = 4 


jfe = 5 
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(r) 

(7,7 


2.9910566246 


-24.021689 


91.591 


645.5 
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fir,. 


-4 9910566246 


28.021689 


-148 830 

_L X v7 • v 7 ' J V 7 
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-2.0000000000 


4.000000 
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2302.2 
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-1 3991056625 


4 674008 
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-0.6008943375 
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20.388 


-200.8 
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a, + a. 


-2.0000000000 


4.000000 
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-296.5 
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(r) 

"fe 


-1.351212983 
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-11.322 


43.164 
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-0.648787017 
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-0.096 
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a, + a. 


-2.000000000 


4.000000 


-11.418 


43.775 
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fit. 
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1 883461 
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2 116539 


-4.8117 
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^k 1 "fc 


-2.0000000000 


4.000000 


-10.2365 


30.172 
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(r) 
fir,. 


-0 6326570719 


1 197694 


-2 94579 

■ «.7 J- .7 1 ' 7 


8.3372 




w 

a. 


-1.3673429281 


2.802305 


-5.96129 


12.1238 




(»•) , (») 
a, + a, 

k 1 re 


-2.0000000000 


4.000000 


-8.90708 


20.4610 



Table 7.4: The fitted coefficient to the series formula of q l J A ()7.3()jl with nn = 0, 2, 4, 6 and 8 

7.3 Quantum numbers of the states with higher TV 

In the previous Sections we presented the spectra for N = 3 and N = 4 particles. One can notice 
that the latter spectrum is much more complicated than the one for N = 3. This complexity 
grows for larger N. Thus, in this Section we discuss only the ground states for AT = 2, . . . , 8 
which we have obtained after solving numerically the quantization conditions (|6.27j) [221 ESI- 

For more than N = 8 Reggeons numerical problems appear. Firstly, the {g^}— space be- 
comes more dimensional, i.e. it has (2iV — 2) real dimensions. This means that the method of 
solving non-linear equations [HHl for higher N needs far more time. Secondly, in the case N > 8 
the series (|6.7|) and 1)6. 17j) become less convergent and recurrence relations for their coefficients, 
which include differences of large numbers, start to require bigger precision. Thus, in order to 
overcome this problem one have to use multi-precision libraries. 

7.3.1 Properties of the ground states 

Performing the numerical calculations we found for each sector from N = 2 to 8 states which 
have the minimum energy. It turns out that these states have common properties for even N. 
However, the properties of states with even N are different than the ones with odd N. 

For even N the ground states are in sector where = so that h = \ + iv h . They are 
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situated on the trajectories where the odd conformal charges vanish, 

?3 = ?5 = • • • = Qn-i = 0, (7.33) 

while the even conformal charges are purely real, 

Im[g 2 ] = Im[g 4 ] = . . . = Im[q N ] = 0. (7.34) 

The minimal energies for these trajectories are at v h = 0. The states are symmetric under (|7.4jl 
and (17.5ft . However, due to (|7.33ft and (|7.34jl the ground states for even A are not degenerated. 
Moreover, the ground states with even A have negative energies En < 0. 

For odd A we have two types of the ground states: the descendent ones and the others that 
have qN ^ 0. 

The descendent ground states are in the sector with = so that h = 1 + iv^, and they 
energies En^n = 0) = vanish. Along all trajectories qN = so these states are compositions 
of (N — 1)— particle states. 

On the other hand, the ground state trajectories for odd A with qN ^ are in the sector 
with rih = so h = \ + iv^. For these trajectories even conformal charges are real, 

Im[g 2 ] = Im[g 4 ] = • • • = Imfeiv-i] = (7.35) 

while the odd conformal charges are imaginary, 

Re[g 3 ] = Re[g 5 ] = . . . = Re[q N ] = 0. (7.36) 

Similarly to the case with even N, these ground states are situated at = 0. However, they 
have positive energies En^n 7^ 0) > 0. Moreover, due to (J7.4I) and (17.5ft the ground state is 
double-degenerated. We show exact values of the conformal charges and the energies for the 
ground states with qN 7^ in Table 17.51 Because of the degeneration for odd iV we have two 
equivalent sets of the ground states with q^ <-> (— l) k qk- 

7.3.2 Energy for different iV 

In order to analyse the dependence of the energy as a function of the number of particles iV we 
plot this dependence in Figure I7T71 As we can see we have two cases: one for odd and one 
for even N. 

For even N the lowest energy is for the ground state with N = 2 and it grows with N as 

E » N- 1.3143- (7 - 37) 
Contrary to the even N case, for odd N we have the highest energy for the ground state with 
N = 3 and it falls down with N as 

-2.0594 

N A - 1.0877 V ; 

In the both formulae for large A we obtain the limit of the ground-state energy as 

\E N \ ~ ^ ^ 0. (7.39) 
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Table 7.5: The exact quantum numbers, g^, and the energy, E N , of the ground state of N 
reggeized gluons in multi-colour QCD. The dispersion parameter, a^, defines the energy of the 
lowest excited states, Eq. (|7.2j) . 
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Figure 7.17: The dependence of the ground state energy, —En/4, on the number of particles 
N. The exact values of the energy are denoted by crosses. The upper and the lower dashed 
curves stand for the functions 1.8402/ (JV - 1.3143) and -2.0594/(iV - 1.0877), respectively. 



7.3.3 Energy dependence on Uh 

The ground states are situated on the ground-state trajectories which are parameterized by a 
continuous parameter v h . Let us consider the dependence of the energy E N {vh) along these 
trajectories, which are plotted in Figure l7~T8l 

We noticed in the previous Section that we have minima of the energy at = and these 
energies satisfy inequality 

E 2 < E 4 < Eq < E 8 < < E 7 < E 5 < E 3 . (7.40) 
These values of E N {y h = 0) coincide with those depicted by crosses in Figure lTTTl Thus for u h , 



80 



12345 12345 

Figure 7.18: The dependence of the energy —E^iy^) / '4 and the "highest" integral of motion 
\qn\/Q2 with q 2 = (1/4 + uf) on the total spin h = 1/2 + ivh along the ground state trajectory 
for different number of particles 2 < N < 8. At large z^, — E 8 > . . . > — _E 3 > — E 2 on the left 
panel and | (fa Ate I < • • • < 193/92 1 on the right panel. Dashed lines correspond to odd N whereas 
solid lines to even N. 

the order of the energy levels tell us that we have anti- ferromagnetic system 2 . In the vicinity 
of i>h = the energies behave as (|7.2jl where measure accumulation of the energies. The 
values of <Tat are shown in Table IT751 

Going to the larger we notice that the energies grow. For ~ 1 we have a quantum 
phase transition. The energies of our system reorder and for the large v-^ the ground-state 
energies of our system satisfy 

< E 2 < E 3 < E 4 < E 5 < E 6 < E 7 < E 8 . (7.41) 

The reason for this behaviour is that for the total SL(2, C) spin h — 1/2 + Wh with the large 
v h the system approaches a quasi-classical regime |H]l E21 IH| , in which the energy E N (u h ) and 
the quantum numbers have a universal scaling behaviour 

E N (u h ) ~ 4Log|g#| , \q N \~C N ul, (7.42) 

with Cn decreasing with N. As can be seen from the right panels in Figures l?~9l and 17. 18[ this 
regime starts already at Uh ~ 2 [23J. 

2 In anti-ferromagnetic systems we have two different sectors: with odd and even N numbers of particles. 
For example: in the simplest anti-ferromagnetic periodic spin-chain with two possible spin orientations: up and 
down, for odd N the lowest energy is when the neighbouring spins have opposite orientations. Since for odd 
N we have always one pair of neighbouring spins with the same orientation the ground states for odd N has 
higher energies than the the ground states for even N. 
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Chapter 8 



Other quantization conditions 



In this Chapter we show for completeness another solution to the Baxter Equation (|3.3|) that 
has been presented in Refs. (23 [2B| by de Vega and Lipatov. In their method two additional 
conditions were assumed: the conformal charges should satisfy Im[z fc gfc] = and the holomorphic 
energy of each from N linearly independent solutions of the pertinent Baxter equation should 
be equal. As was explained in Ref. (2HJ these conditions are too strong. 



8.1 De Vega and Lipatov's solution 

In the papers [23 EH] the authors consider the Baxter equations ()3.3jl and 1)3.4)1 with complex 
spins (s — l,s — 1). The equal values of the complex spins cause that the Baxter equations in 
the holomorphic and anti-holomorphic sectors are identical. However, is this case the complex 
spins (s — l,s = 1) do not satisfy 1)2.58)1 so we have to use the different form of the scalar 
product given by p. 42)1 . 

It has been shown in Ref. |22j that the simplest N— Reggeon solution to the Baxter equation 
()3.3)) may be conveniently written in a form of a sum over poles of orders 1 up to iV — 1 situated 
in the upper semi-plane of complex variable u: 



oo p(N-2), \ oo 
r=0 v ' r=0 



a r (q) K(q) z r (q) 



(u — ir) 1 ^^ 1 (u — ir) N ~ 2 (u — ir) 



(8.1) 

where q = {g 2) 13, ■ ■ ■ , In} an d P^h~q'\ u ) are polynomials in u of degree N — 2, whereas 
a r (q),b r (q), . . . are some residual coefficients. Substituting 1)8.1)1 to 1)3.3)1 we obtain the re- 
currence relations between the polynomials -P r .^„ (u) which allows us to calculate them suc- 



(JV-2), 



u 



cessively starting from PQ, hq ' 

The solution 1)8.1)1 is normalized by the constraint 

1™<,; 2) W = ^(?) = 1 (8.2) 

and the remaining coefficients of the polynomial Po^^iu), or equivalently the coefficients b r (q), 
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c r (q), • • • , are calculated from the condition 

]jmQ( N - 1 \u',h,q)~u h - N (8.3) 



which ensures (2jj that the solution (j8.1|) satisfies the Baxter equation (|3.3j) for u — > oo. The 
condition ()8.3|) leads to 

(u — ir) 



lim u N ~ \ r ' h ' g ^ = (8.4) 



which gives N — 2 equations and with 1)8. 2j) fixes all the polynomial coefficients 2 \u), or 
equivalently a (q), b {q),.... 

The second independent solution may be written in terms of the first one as 

oo p(N-2) ( , 

Q<°>(«; h, q) = Q^i-u; h, -q) = £ , (8.5) 

r=0 v 7 

where -q = (q 2 , -q 3 , • • • , (-) n Qn)- 

It turns out, that there is a set of Baxter functions Q^(u) with £ = 0, 1, . . . , N — 1 which 
have poles both in the upper and lower half— u planes: 



qM(«;M = E 



r=0 



fit — ir)* f— ix — zr)^ 1 * 



(8.6) 



where the polynomials P^ h q (u) and P^^ q {—u) are fixed by the recurrence relations fol- 
lowing from the Baxter equation and from the condition that the solution should decrease at 
infinity more rapidly than u~ N+2 . 

Using these functions in both holomorphic and anti-holomorphic sectors it is possible to 
construct the Baxter function 

Q qrq (u,u) = J2Ct,iQ {t \u ] h,q)Q^ l \u ] h,q) (8.7) 

t,l 

choosing the coefficients C t) i, appropriately. 

In this method the coefficients C t ,i are fixed by the normalization condition [27J. It leads to 
the requirement that Q q ^(u,u) do not have poles at 

u = im and u = —im for \m\ > . (8-8) 

The pole at it = u = is removed by the measure of the (s = l,s = 1) scalar product 



(12.421) . Because in the bilinear solution we have products of the poles (u — ir)~ s and (u — ir')~ s 
we expand Q q , q (u,u) around the spurious poles (18.81) and we equate the residual expansion 
coefficients to zero. This results in the quantization conditions for the conformal charges q and 

q- 

The authors of Refs. [23 EHJ claim that the same conditions are given by imposing equality 
of holomorphic energies 

* = i ^^[ uNp S£M] (8 - 9) 
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for all solutions Q^(u). However, our system is two-dimensional and the holomorphic energy 
is not a physical observable and it may depend on a calculation method. The real observable 
is the total energy defined as 

E = 2i lim ^^ln \(u-i) N -\u-i) N -\ N u N Qaa(u,u)] . (8.10) 
8.1.1 Three Reggeon states 

For simplicity, let us consider the N = 3 Reggeon case. The Baxter equation 1)3.3)1 with s — 1 
reduces to 

B 3 (u;h,q 3 ) = [2u 3 - h(h - l)u + q 3 ] Q(u; h, q 3 ) 

-(u + ifQiu + z; h, q 3 ) - (u - ifQiu - z; h, q 3 ) = . (8.11) 

Here the authors of Refs. [23 [28] claim that iq 3 should be real in order to obtain single- valued 
wave-functions. However, we have already shown ()5. 23)1 - 1)7.11)1 that it is possible to construct 
a single- valued function without this assumption. 

To find polynomial coefficients in Q^'(u,u) we define auxiliary functions f r , which look as 
follows 

oo 

h(u;h,q 3 ) = 

1=0 

h(u;h, q ) = ± a -^Mi. (8.12) 

^ — ' — IU — I 

1=0 

They help us to impose the conditions coming from behaviour of Q(u; h,q 3 ) at u — > oo. Sub- 
stituting ()8.12j) to ()8.11)1 we obtain the recurrence relations for the residua of ()8.12)) : 

(r + l) 3 a r+ i(h,q 3 ) = [2r 3 + h(h - l)r + iq 3 ] a r (h,q 3 ) - (r - l) 3 a r _i(/i, q 3 ) 

(r + l) 3 b r+1 (h,q 3 ) = [2r 3 + h(h-l)r + iq 3 ] b r (h, q 3 ) - (r - l) 3 b r ^{h, q 3 ) 

+ [ 6 r 2 + h{h - 1)] a r {h, q 3 ) - 3(r + l) 2 a r+1 (h, q 3 ) - 3(r - l) 2 a r _n(/i, q 3 ) . (8.13) 

We can choose, 

a (h,q 3 ) = 1 and 60(^,^3) = (8.14) 

what gives as 

ai(h, q 3 ) = iq 3 , h(h, q 3 ) = h(h - 1) - 3z'g 3 . (8.15) 

Thus, we can calculate all the coefficients a r (h,q 3 ) and b r (h,q 3 ). 

Now, we can construct Q®(u,u) as a linear combinations of the auxiliary functions 1)8.12)1 
in a form 

Q {2) (tt; h, q 3 ) = f 2 (u;h,q) + B(h,q 3 )fx(u;h,q) , 

Q^\u;h,q 3 ) = f 1 (u;h,q) + C(h,q 3 )f l (-u;h,-q 3 ), (8.16) 



a-i{h,q 3 ) bi(h,q 3 ) 
(—iu — I) 2 —iu — I 
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where B(h,q 3 ) and C(h,q 3 ) are coefficients we will fix using the condition (I8.4jl by demanding 
that the Baxter equation is satisfied at infinity. Taking the limit u — > oo of (jS.llj) with (|8.16|) 
we obtain 

B 3 (oo;h,q 3 ) = [h(h - 1) - 2] lim [u Q(u; m, fi)] . (8.17) 

u^oo 

Applying (|8.17|) to ()8.16|) we get the coefficients in a form 

B(h,q 3 ) = - ^° br [ h h ,q3 \ and C(h,q 3 ) = ^ a f> q3 \ . (8.18) 

Therefore, the solutions Q( 2 \u; h,q 3 ) and Q^'(u; h,q 3 ) are completely determined. 
Due to the symmetry (I7.4J) we can construct the third solution 

Q^(u;h,q 3 ) = Q^\-u;h,-q 3 ). (8.19) 

Moreover, if we multiply a solution of the Baxter equation by a periodic function of %u with 
period 1 the result is also a solution of the Baxter equation (|8.11jl . For instance 

(u; h, q 3 )n cot(ivra) (8.20) 

is also the solution to (jS.lljl and thanks to cot(i7rw) it has second order poles at u = im for 

1 771 1 > 0. 

All these solutions are related by the equation [23 [2H] 



[X{h, q 3 ) - coth(vra)] (u; fc, g 3 ) = Q (2) («; K fe) - C(/i; g 3 )Q (0) (u; fc, g 3 ) , (8-21) 

where g 3 ) can be obtained by performing the Laurent expansion of (I8.21JI around the poles 
(|8.8[) . In Refs. [23I2E1 it is shown that it is possible to construct the solution (|8.7j) with h = h 
and g 3 = — g 3 for which poles at (|8.8|) disappear. It has the following form 

Q q ,q(u,u) = Q {2 \u- h,q 3 )Q^(u;h,q 3 ) - Q (0) (m; h, g 3 )Q (0) (u;h,q 3 ) . (8.22) 

Moreover, in Refs. [23 EH| the authors assume that the condition X(h, q 3 ) = gives 
quantization of q 3 . Accidentally, for solutions with X(h,q 3 ) = the holomorphic energies 
(I8.9JI of the solutions (u; h, q 3 ) and it cot(iwr)Q^ (u; h, q 3 ) are equal what gives equivalent 
condition to X(h, q 3 ) = which looks like 

B(h, , 3 ) - h(h, ,3) - A- + -E + £ = . (8.23) 

^3 l Qsj^ 2 r-l %q 3 r + 1 



However, the way of cancelling poles presented in [SUES] is not unique. Indeed, the condition 
(|8.23|) is too strong and it is only satisfied for some solutions with Reg 3 = 0. There are other 
solutions, for example with Img 3 = 0, for which the holomorphic energies are not equal. 

In Refs. [23 [2H1 for N = 4 Reggeons the quantization conditions are also introduced by 
the equality of the holomorphic energies. Similarly to the N=3 case, the numerical values of 
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the conformal charges which result from these conditions agree with Refs. [23J only for h = h, 
Reg 3 = and Img 4 = 0. 

Contrary to the solutions described in (2HJ the solutions presented in this Chapter have a 
simple pole structure, so they look much simpler than those defined by (I6.1j) - (l6.16j) . However, 
the series in (|8.12j) are much slower convergent than the series defined in (|6.7j) and (|6.17j) . More- 
over, in [2H1 the quantisation conditions ()6.27|l come from single-valuedness and normalization 
of the Reggeon wave-functions explicitly, whereas in |27| l28j except for the normalization con- 
dition, the quantization conditions (|8.23J) follow from the quality of the holomorphic energies 

(JEU). 

Thus, if one wants to find a spectrum of the conformal charges numerically the method 
presented in |23| is much more convenient. 
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Chapter 9 



Anomalous dimensions 



In the preceding Chapters we used the reggeized gluon states to describe the elastic scattering 
amplitude of strongly interacting hadrons. It turns out that the Reggeon states might be also 
used to describe deep inelastic scattering of a virtual photon 7*((? At ) off a (polarized) hadron 
with momentum p^ |26[ 128] . In this case we also perform calculations in the limit of the low 
Bjorken x but in the region where 

M 2 « Q 2 « s 2 = (p„ + q,f = g2(1 ~ X) , (9.1) 

with 

Q 2 = -q 2 M 2 =p 2 . (9.2) 

Notice that in the previous case (|2.2|l there was only one hard scale present, namely the mass 
of the scattering particles which were assumed to be similar. 

In the limit (|9.1|) . similarly to (|2.2j) . the moments of the structure function F 2 (x,Q 2 ) = 
F(x, Q 2 ) may be rewritten as a power series in the strong coupling constant a s : 

F(j, Q 2 ) = / dx xi- 2 F(x, Q 2 ) = £ ^~ 2 Fn{j, Q 2 ) , (9.3) 

J ° N=2 

where a s = a s N c /ir. 

The leading contribution comes from the Reggeon states with the intercept j — ► 1 and in 
this limit 

where the conformal charges are denoted by q = (q2,Q2i ■ ■ ■ 1 Qn } Qn) an d P3(Q), ft${M) 

are the impact factors of virtual photon and scattered hadron, respectively, while the energies 
E N are the eigenvalues of the Reggeon Hamiltonian (|2.31j) . 

Performing the operator product expansion (OPE) one may expand the moments of the 
structure function F^(x, Q 2 ) in inverse powers of hard scale Q 

hj,Q 2 )= E ^Y, c r^^s{Q 2 )){p\oi J mp), (9.5) 

n=2,3,... ^ a 
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where the expansion coefficients are the forward matrix elements of Wilson operators 0*^(0) 
of increasing twist n > 2 and they satisfy renormalization group equations 

Q 2 ^M°IMp) = <U) (p\o a n,Mp) ( 9 - 6 ) 

with the anomalous dimensions of QCD 7^(j). The parameter a enumerates operators with the 
same twist n. It also appears in (j9.5j) where C%(j, a s ) denote the coefficient functions. 
The anomalous dimensions 7^(j) may be written as a power series 

oo 

J a nti)=Y.^nU){*s(Q 2 )n k - (9-7) 
k=l 

For small Bjorken x (j — > 1), the moments of the structure function Fiy(j,Q 2 ) take a form 



1 fM^ n ~ 2 ~ 2 ^ 



F(hQ') = —, U 2_ C '»(J>'-W)) I-] • (9-8) 

n=2,3,... a 



It turns out that solving the Schrodinger equation ()2.34|) we are able to find poles which 
give the main contribution to (|9.4|h Moreover, combining ()9.8|) and (|9.4|) we are able to find 
the dependence of 7^ n (j) on the strong coupling constant near j — > 1. In this Chapter we show 
the way one can find this dependence using techniques developed in the previous Chapters. 

9.1 Expansion for the large scale Q 

The moments of the structure function are defined in (19. 8 j) . In the limit j — > 1 they can be 
rewritten as (|9.4|) where the impact factors are given by 

$(Q) = / d 2 z (%*\* q (z )) , /#(M) = J d 2 z (* q (2o)\%) (9.9) 

and the functions ^q(z ) are orthonormal with respect to the scalar product ()2.35|) : 

<* q (zo)|%(4)> = / nd 2 z k * q ({z}-z )^ q i{z}-z* )y = 5V\z -z' )5 qq/ , (9.10) 
^ fe=i 

where 5 qq > = 5{vh — v' h )5 nhn > 8g£/ with £ = {£2, £3, ■ ■ ■ ,£2(n~2)} defined in (|7.1j) . Due to the 
scaling symmetry of the reggeized gluon states 

*q(\zi, Xz 2 , ...,Xz N )= \ 2 - h -H q {z x , z 2 , . . . , z N ) (9.11) 

we are able to calculate dimensions of the impact factors (|9.9j) 

$(Q) = cfyQ- 1 -* 1 * , Pp{M) = M- 1+2i » h , (9.12) 

where and Cp are dimensionless. Substituting (I9.12JI into (19 .4jl and expanding in powers 
of the ratio (M/Q) 

1 /-oo / 7t/f \ _1 + 2i "h 



£ n h >0 
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where q = q(uh',h n ,£). We shall calculate the integral in ()9.13)) by performing analytical 
continuation in v^. Next, we shall close the integration contour in the z/^— complex plane and 
integrate by summing residua closed by the contour. The corresponding poles come from the 
denominator of (19.131) so they are placed at Uh where 

j-l = a a E N {q(u h ;n h ,£))/A (9.14) 

is satisfied. Here the parameters nu and I = {^1,^2, • • • ,&2(N-2)} are integer numbers and 
they enumerate the quantizated energy levels. The formula for (19.13)) is valid for j > = 
1 - a s min q E N (q)/ 4. 

In this way a contribution to the moments of F(x, Q 2 ) coming from a given pole in takes 
the following form 

F N (J,Q 2 )-^ . (9-15) 



Q 2 \Q 

where we show the dependence on j explicitly. Comparing exponents in (|9.8j) and (J9.15I) we 
obtain relation 

7n(j) = {n- l)/2 - iv h {j) = [n- (h(j) + h(j))}/2 , (9.16) 

where h(j) and h(j) are the SL(2, C) spins defined in (12.601) and j n (j) is the anomalous dimen- 
sion of the twist— n operator 1)9.6)1 . 

Now we can calculate j n (j) which with a s (j) — > behave as j n (j) — > 

ln(j)=^-^ + 0(a s 2 ). (9.17) 

Substituting ()9.14)) into ()9.17j) one finds that j n (j) — > corresponds to E N {v h ) — > 00. Thus, 
we can combine the expansion of the energy E^{yh) around its poles in the complex v^— plane: 



E N (q) = -4 



C-l 

h c + ci e + . . . 

e 



(9.18) 



with the small anomalous dimensions j n (j) = ~~ e w h ere according to (|9.17j) iv h = iz/^ ole + e. 
Inverting the series 1)9.18)) and using (|9.14)1 we obtain 



7n(j) = "C-l 



+ . . . 



(9.19) 



where coefficients Ck = Ck(n,rih,£) in Eqs. ()9.19)) and (19.181) are identical. Thus, evaluating 
these coefficient for the energy ()9.18)) we may calculate coefficients for the anomalous dimension 
expansion in the coupling constant a s . 

Moreover, we can notice that the positions of the energy poles 

E N (q) ~ n 1W9 (9.20) 

iv h - (n - l)/2 

determines the twist n of 

iv h = (n-l)/2 with n > iV + n ft (9.21) 
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with the number of Reggeons N and the Lorentz spin n^. 

For N = 2 formula for the energy (THJ 0] is given by an analytical expression 

E 2 (v h ,n h ) = 4 (Jl^^ + i v A + Jl±pL- iv ^\- 2^{l) \ , (9.22) 



where rih are nonnegative integers. In this case the anomalous dimensions for the leading twist 
n = 2 have been calculated in Ei 



7 2 (J) = + 2C(3) 4 + 2C(5) 6 + O(of) , (9.23) 

where is the Riemann zeta-function. The general formula for N = 2 and arbitrary n > 2 
and nfe has been obtained in [2~5] : 

7.0') = + (2^(1) - V(l + Kl + n) - V(l + n)) f ' + O(o») . (9.24) 

9.2 Analytical continuation 

The reggeized gluon functions ^f(z*k) defined in 1)2. 54j) are normalized with respect to the scalar 
product 1)9.101) . Performing the analytical continuation is Uh it turns out that the normalization 
condition 1)2 .63)1 is no longer satisfied. The quantization conditions (I2.63|) and 1)6.271) become 
relaxed. Thus, conditions q k = are not necessarily satisfied so that h ^ 1 — h*. However, 
we still have the conditions coming from the single-valuedness of the Reggeon wave- functions. 
The latter condition (16.271) ensures that the two-dimensional integrals (I9.9jl are well defined. 
In (I3.18j) the energy E^iyh) is a smooth function of real v^. After the analytical continuation 
into the complex v^— plane the energy spectrum exposes pole structure at Re[^] = 0. 

In Quantum Mechanics the problem of analytical continuation of energy E(g) as a function 
of a coupling constant g has been studied in various models [63 | 165]. In our case Vh plays a 
role of such a coupling constant g. 

Generally, the energy E(g) is a multi-valued function of g. It turns out that the number of 
energy levels for real g is equal to the number of branches in the complex g— plane. In order to 
analyse the global spectrum one may glue together different sheets corresponding to its branches 
and study E(g) as a single- valued function on the resulting Riemann surface. Thus, the models 
may have a complicated structure of spectrum in the complex g— plane. Additionally, these 
spectra may consist of a few disconnected parts due to some additional symmetry [64J. 

9.2.1 Two level model 

In order to understand the analytical continuation more clearly let us consider the simple model 
with two energy levels e± and e^- 

H=( ei 9 ) , (9.25) 
V 9 £2 
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where we add interaction by introducing a small coupling constant g. 

The energy eigenvalues of the interacting Hamiltonian (|9.25|) are given by 

E ± to ) = fi±« ± y / (^) 2 + 9 2 ■ (9.26) 

Here we have two branches E + (g) and E_(g) which are shown in Fig l9.ll where for illustration 
we have chosen with e\ = and e 2 = 1. For the real g the energy levels E±(g) do not cross. 




Figure 9.1: Spectral surfaces for two energy level model. On the left panel energy levels as a 
function of a real coupling constant g. On the right panel real energy surfaces as function of 
complex g 

However, we have two branching points at g± = ±|(ei — e 2 ) for which E + (g±) = E-(g±) = ei 1 2 e2 . 
The energy is a smooth two- valued function of complex g. The spectral surfaces cross each other 
for Re[g] = where —ig > —ig + or — ig < —ig-. Going around these branching points we can 
pass from one surface to the other one. In this case, the spectral curve, which is a function 
defined on one complicated Riemann surface, consists of two branches E + (g) and E^(g). 

One has to add that in this model the energy does not have any poles which occur in our 
SL(2, C) Heisenberg model. 

9.2.2 SL(2, C) Heisenberg model for two Reggeons 

The energy for two Reggeons with nonnegative is defined by (I9.22J) . We notice that analytical 
continuation of (|9.22|) exhibits the poles at 

n — 1 

iv h = ± — - — with the twist n > 2 + n h . (9.27) 

The residua of these poles are related to the anomalous dimensions given by ()9.24|) . The leading 
twist n = 2 corresponds to the pole at v h = —i/2 with n h = 0. 
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The energy (|9.22j) does not have any branching points. Thus, the spectral surfaces numbered 
by the Lorentz spin do not mix with each other in the complex z/^— plane. The spectral 
surfaces are functions defined on trivial separated Riemann surfaces. 

9.2.3 SL(2, C) Heisenberg model for more Reggeons 

For the case with more than two Reggeons the energy spectrum is much more complicated. 
The energy E^iy^ nt,£) depends not only on Vh and rih but also on 2(iV — 2) integer numbers 
I = {£1,^2, ■ ■ ■ ,h(N-2)}- Similarly to the N = 2 case the energy surfaces for N > 2 possess 
poles for imaginary v^. Additionally, they also have branching points. Thus, the energy spec- 
trum consists of infinite number of spectral surfaces, which are functions of complex with 
complicated Riemann structure. 

The energy in (|3.18|) is a smooth function of real u^. In this region we have quantization 
conditions coming from the single- valuedness of the Reggeon wave-function as well as the nor- 
malization condition which ensure that Hamiltonian (|2.32j) is hermitian. The latter condition 
gives conjugation relations between the conformal charges 

q k = q k * with k = 2, . . . , N. (9.28) 

Performing the analytical continuation into the complex z/^— plane we have to relax the quanti- 
zation conditions ()2.63j) . The single- valuedness condition survives. However, the normalization 
condition is not valid anymore. Thus, the relations (|9.28|) in complex v h are not valid, so that, 
q k and q k may be independent. 

We notice that the N— Reggeon spectrum contains the states with given rih and £ for which 

q k {vh\ rih, £) = MkWh, n h , I). (9.29) 

These states satisfy (|9.29j) on the whole complex z/^— plane. They are of utmost interest because 
they give the higher order contribution to the structure function (|9.3|) . Other interesting states 
are the descendent states. 

For iV > 3 the energy E^{y} l \nh- I () is a multi-valued function of parameterized by 
and t which number the spectral surfaces. In the complex v^— plane the spectrum possesses 
poles at 

Tl — 1 

ifh = i — - — with the twist n > N + n h . (9.30) 

It turns out that it also possesses the infinite set of branching points which are square-like type 
similarly to (|9.26|h In order to construct the complex curves En^v; rih,£) we apply Eqs. (|3.18|) 
developed in Refs. [2B, 23j. After analytical continuation we get 

E N = [e(h, q) + e(h, -q) + (e(l - h*, q*))* + (e(l - h\ -g*))*] , (9.31) 

where the SL(2, C) spins h and h are given by ()2.6flj) and the conformal charges are denoted 
as q = {q k }, q = {q k }, -q = {{-l) k q k } and -q = {(-!)*%}. As one can see E N (-q, -q) = 
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Ejsr(q,q). Since for real Vh we have h = 1 — h and q = q* the energy E^{q,q) takes real values 
on the real u^— axis. 

The function e(h, q) in (|9.31jl is defined for arbitrary complex h and q by 



e(h, q) = i-^ In [e N Q(i + e; /i, g)] 



(9.32) 

e=0 



where Q(u; h, q) has a form 



i 



Q(u;h,q)= / dzz m - l Qx{z). (9.33) 



Going further, Qi(z) satisfies the differential equation (16 .2|) which comes from the Baxter 
equation (|3.3j) . Around z = 1 the function Qi(z) corresponds to one of the solutions in 1)6.16)) 
with asymptotics Qi(z) ~ (1 — zY~ h ~^ and around z = from (|6.6|) with Qi{z) ~ ln^ -1 2. It 
is convenient to normalize Q(u; h, q) as 

Q<? + e,h,q) = ±- ^ £ -^4 + ° (^) ■ (^) 

where the function e(h,q) appears in (|9.31|l . 

Thus, solving the quantization conditions (|6.27|) and using (I3.18jl and (19.31^ we are able 
to evaluate the energy E N {y h ) on the whole complex u h — plane. Then, we find the poles 
which appear only for Re[^] = 0, and we evaluate expansion coefficients of E N (v h ) around its 
poles. Next, using (I9.18|) and (I9.19|) we the calculate expansion coefficients of the anomalous 
dimensions (|9.19|) in the strong coupling constant a s . These numerical results will be presented 
in the next sections. 



9.3 Spectral surfaces 

The two-Reggeon energy spectrum is described by the analytical formula (|9.22)l . It is a mero- 
morphic function of Uh and depends on the integer conformal Lorentz spin nh- It possesses 
poles at JH2H). 

For higher iV > 3 we do not have an analytical formula, so we have to calculate the spectrum 
numerically using Eqs. (|9.31|) - ()9.34|) . The energy E N {y h ^n h ,t) is a multi-valued function and 
its branches are numerated by Lorentz spin Uh and integer parameters I = {£1, £2, ■ ■ ■ , ^2(^-2)}- 
These spectra, comparing to the N = 2 case, except the poles at (|9.3()L have additionally an 
infinite set of branching points at z/^ r,fc = v h which allows us to glue some branches along their 
cuts. This causes that the spectrum consists of a set of complicated spectral curves. Thus, we 
consider the energy E^iyh.n^t) as a multi- valued function on complicated Riemann surfaces. 
The states from a given Riemann surface have the same quantum numbers, i.e. Lorentz spin 
Uh, C— parity, quasimomentum, Bose symmetry. These numbers impose additional rules for 
parameters £ which describe a single spectral curve on one Riemann surface. 
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Figure 9.2: The energy spectrum of the A/" = 3 Reggeon states E 3 (i/h;nh,£) for = and 
£ = (0,£ 2 ), with £ 2 = 2, 4, . . . , 14 from the bottom to the top (on the left). Analytical 
continuation of the energy along the imaginary v^— axis (on the right). The branching points 
are indicated by open circles. The lines connecting the branching points represent KeE 3 . 

Let us consider the A^ = 3 Reggeon state along the spectral curves labelled with £i = 
and £2 =even positive numbers. This gives condition q k + q k = 0. An example of these states 
is presented in Figure l9~2l On the left panel we plot the energy levels along the real Uh— axis 
whereas on the right panel we continue them analytically. Here we show the states for imaginary 
u h only. 

For real Vh the energy is real and the energy levels are smooth functions of real v^. Even 
though that we have only one slice of the spectral curve we are able to observe that, firstly, we 
have poles at it / l° le = 3/2, 5/2, 7/2, 9/2, 11/2 and, secondly, different energy levels collide at 
= 2.70, 3.29, 4.73, 5.34, 5.74, 5.94. which are denoted by circles in Fig. 19.21 Moreover, we 
have also branching points at complex iv^ = 1.723+i .248 where the ground state £2 = 2 collide 
with the next state £2 = 4. Thus, all these energy levels are glued in the complex z/^— plane 
and they form one multi-valued spectral surface. 

The complex branching point values imply that the contribution of the corresponding 
square-root cuts to the integral 1)9.4)1 scales as l/[Q 1+2lu h T hi 3 / 2 Q] and, therefore, it breaks the 
OPE expansion 1)9.5)1 . Although such corrections are present in 1)9.13)1 for given n h and £, they 
cancel against each other in the sum over all states [26] . 

Another interesting example of spectral surfaces are N = 3 Reggeon states with the Lorentz 
spin Uh — 1 and q 3 = q 3 = 0. They are descendent states and their energy levels coincide 
with the levels for A^ = 2 Reggeon states. Performing analytical continuation we notice that 
qai^h) = lz{vh) = and E 34 (u h ) = E 2 {v h , 1) on the whole complex i^-plane, where E 2 {v h \ 1) 
is defined by (EH2l . 

Therefore, the energy of this state E 3 ^{vh) is a single-valued meromorphic function on the 
complex v h — plane. At — the energy E 3 ^{^h) vanishes and at this point we have the 
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"physical" ground state for the system of N = 3 Reggeons. Similarly to the N = 2 case, the 
descendent states have poles at (|9.27|) . 

In the above example of the N = 3 Reggeon states we have presented two special states 
E%(vh — 0, Vh — 1) = and E 3 (uh — 0, Uh = 0)/4 = 0.24717, which are considered as the 
ground states for the odderon in QCD with the intercept jo = 1 — ol s E^ O033J. We notice 
that the minimal twist corresponds to the position of the energy pole in the complex v^— plane 
which is the nearest to the origin, q k = q k = 0. 

Thus, for rih = 1 we have the minimal twist n min = 3 whereas for rih = we have the 
minimal twist n m i n = 4. Therefore, the leading contribution to the structure function (|9.3j) 
with the intercepts jo = 1 and jo < 1 scales at large Q 2 as ~ 1/Q 3 and ~ 1/Q 4 , respectively. 



9.4 Energy poles and the anomalous dimensions 

In order to evaluate the anomalous dimensions of N Reggeon states around 7 a (a s = 0) = 
we work out the Laurent expansion of the energy function around its poles (|9.18|) . As we said 
before the energy possesses the poles at iu h — (n — l)/2 with the twist n > N + rih and the 
Lorentz spin rih. 

Firstly, let us consider the descendent states for iV = 3 Reggeons with n h = 1. They are 
related to N = 2 Reggeon states, so they have equal energies E^ d {v h ) = E 2 (vh, 1) defined in 
(I9.22L The pole which is the closest to the origin is located at ivh = 1. Thus, {h, h) = (2, 1) 
and the twist n = h + h = 3. Expanding (I9.22J) around this pole we get 

£ 3 ,d(2 + e) = 4 + 1 - e - (2C(3) - 1) e 2 + . . j . (9.35) 

Next, using (j9.35|) and (|9.18j) with (|9.19jl we find the twist-3 anomalous dimension corresponding 
to the odderon state with the intercept jo = 1 as 

7 f-'0') = - (-^) 2 + (2C(3) + 1) (-^J + . . . , (9.36) 

where the subscript and the superscript denote the twist and the number of Reggeons in the 
given state, respectively. 

For the N = 3 case with g 3 ^ we do not have analytical formula for the energy. However, 
applying (|9.31|) - (|9.34j) we are able to fit numerically the expansion coefficients of the energy 
around its poles. 

Let us consider the spectral surfaces with rih = 0, which we show in Fig. 19.21 Applying 
(|9.31jl - ()9.34|) . we expand the energy in the vicinity of the poles at h — 1/2 + iv h = 2,3,4,5 
with E 3 = E 3 (h + e): 

£ 3 (2 + e) = 4 (V 1 + ~ - ~ e+ 1.7021 e 2 + . . . 
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E 3 {3 + e) = 4 ( 2e- 1 + — - 1.6172 e + 0.719 e 2 + . . . 
^ a) (4 + e) = 4 fe _1 + ^ - 0.6806 e- 1.966 e 2 + . . . 
^ b) (4 + e) = 4 ^2 + ^ -3.2187 e + 3.430 e 2 + . . . 



^ a) (5 + e) = 4(2e- 1 + — - 2.0687 e+ 1.047 e 2 + ... 



E ( 3 h) (5 + e) = 4 ^2 e" 1 + ^ - 2.4225 e + 0.247 e 2 + . . . J . (9.37) 

Here e — > and ellipses denote C(e 3 ) terms. Since the energy is defined on complicated Riemann 
surfaces, therefore for one Riemann surface at the same given we have many branches and 
at some given we may have many poles. Thus, we introduce the superscripts (a) and (b) to 
distinguish different branches of spectral curves as well as different poles at the same Vh- 

In this case the nearest pole is at Uh = 3/2. Therefore, the minimal twist equals n = h+h = 
4. Now applying (|9.18jl and (|9.19|) we find from the first relation in (|9.37|) the anomalous 
dimension corresponding to the odderon state with the intercept jo < 1 as 

^U) = ^ - 5 (^r) 2 - j (^t) * - wm C^i) 4 + " " 

Similarly, we may calculate the other anomalous dimensions form (|9.37|) . However, to save the 
space we do not present them here. 

All expansion series formulae for the energies around its poles may be written as 

E 3 (h + e) = 4 + 7 (/0 + 0(e)) , (9.39) 

where the parameter R is integer and it is equal to 1 or 2. Moreover, the function j(h) is a 
rational number and is related to the energy of the Heisenberg SL(2, R) magnet model. More 
details concerning the properties of j(h) can be found in Refs. 66 67L Because it is not the 
topic of this thesis we omit this interesting subject. 

Now, let us consider cases with N > 4 Reggeon states. For even N the main contributions 
comes from the sector with rih = 0. The closest poles to the origin are located at iv h = (N—l)/2 
so the lowest twist n = N. These states build Pomeron states. For instance, the expansion 
series of the energy for N = 4 and N = 6 states which give the highest contribution look as 
follows 

/2 1 

£ 4 (2 + e) = 4f - + l--e- 1.2021e 2 + ... 

£ 6 (3 + e) = 4^ + ^-^e- 0.238e 2 + ... ) . (9.40) 
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Applying Eqs. (|9.18j) and (|9.19L we obtain the twist— A anomalous dimensions of the A— Reggeon 
states 

2 / _ \ 4 



^—(J) = + - 13.6168)-^ ) 

— / — \ 2 / — \ 3 
(JV=6),- '>« / \ . / a s \ / n„ 



ii'(S) = + + " 33 2 HAJ + " <9 ' 41) 

The higher corrections are large and grow with A so the expansion series (19.40)1 has finite radius 
of convergence and its value decreases with A. This agrees with the fact that the intercept of 
A— Reggeon state scales at large A as [jjv — 1| ~ 1/A [23] . 

In the case with odd A, we consider two interesting sectors: with n h = and rih = 1. In 
the first sector the pole closest to the origin is located at ivh = A/2. These gives the minimal 
twist n = N + 1. For A = 3 we have written expression for the energies and the anomalous 
dimensions in ()9.37|) and (|9.38ft . For A = 5 they look like 



A3 7 101 n \ . , 

£ 5 (3 + e) = 4 (- + -- — e- 0.1136 e 2 + ... .J , (9.42) 

™ -j% + KA) 2 -K^) 3 - i3 - 032 (A) 4+ -- 

These states possess the odderon quantum numbers and have intercepts jo = 1 — a s E^ lower 
than 1. 

In the second sector, with odd A and n h = 1 the minimal twists are smaller and they are 
equal to n — N. It corresponds to the poles of the energy at iv h = (A — l)/2. At A = 5 we 
have the Laurent expansions as follows 

£ 5d (3 + e) = 4 ( 3+ 2 ^ + 1.36180 - 0.4349 e - 0.315 e 2 + . . . ) , (9.43) 



3 - 1 + 3,56524 G^i) 2 + 1,8743 (A) 3 " 11,219 (A) 4 



7 f =5) (j) = , Q % +3.565241 -7——— ) +1.8743(^-1 -11.219(^-1 +.... 



For A = 3 the results are shown in (|9.35l) - (|9.36l) . The poles belong to the same spectral 
surfaces that contain the ground states for given odd number of particles A. Because, these 
ground states are descendent ones, they have minimum energy = 0, the intercept jo — 1 
and finally they have the odderon quantum numbers. 
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Chapter 10 

Summary and Conclusions 



In this work we have considered the scattering processes in the Regge limit 1)2.2)1 where the 
compound reggeized gluon states, i.e. Reggeons, propagate in the t— channel and interact with 
each other. We have performed calculations in the generalized leading logarithm approximation 
(GLLA) [3 H3E1, in which a number of Reggeons in the t— channel is constant. In this approach 
the reggeized gluon states appear in two different kinematical regions. Thus, we have considered 
both of them, i.e. the elastic scattering processes 1)2. lj) as well as the deep inelastic scattering 
processes ()9.1)1 . We attempted to find a scattering amplitude of the above processes with multi- 
Reggeon exchange. However, a structure of reggeized gluon states as well as their properties 
have turned out to be so reach, complicated and interesting that in this work we have focused 
on description of the Reggeon state properties as well as on analysing the spectra of the energy 
and integrals of motion. 

In order to simplify the problem we have applied the multi-colour limit [10], which makes 
the N— Reggeon system SL(2, C) 1)2.26)) symmetric 1)2.271) and completely integrable. In this 
limit the equation for the N— Reggeon wave-function takes a form of Schrodinger equation 
()2.34)) for the non-compact XXX Heisenberg magnet model of SL(2, C) spins s [68, 44, 69|. Its 
Hamiltonian describes the nearest neighbour interaction of the Reggeons [13 [T4j propagating 
in the two-dimensional transverse-coordinates space (12.91) . The system has a hidden cyclic and 
mirror permutation symmetry 1)2.68)) . It also possesses the set of the (N — 1) of integrals of 
motion, which are eigenvalues of conformal charges [38j, qu and q k , ()2. 66)1 - 1)2. 67)1 . Therefore, 
the operators of conformal charges commute with each other and with the Hamiltonian and 
they possess a common set of the eigenstates. 

Eigenvalues of the lowest conformal charge, q2, may be parameterized (15.21) by the complex 
spins s and the conformal weight h, where h can be expressed by the integer Lorentz spin 
and the real parameter v h related to the scaling dimension 1)2. 60)) . Solving the eigenequation 
for q 2 we have derived ansatzes for N— Reggeons states with an arbitrary number of Reggeons 
N as well as arbitrary complex spins s (I5.ll)) . Since the N = 3 Reggeon ansatz separates 
variables, the q%— eigenequation can be rewritten as a differential equation of a Fuchsian type 
with three singular points ()5.16)1 [18]. We have solved this equation by a series method. Gluing 
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solutions for different singular points, (|5.24jl . (|5.29jl and (jf>.33j) . and taking care for normalization 
and single-valuedness of the Reggeon wave-function 1)5. 23j) we have obtained the quantization 
conditions ()5.41j) - l|5.42j) for integrals of motion (q 2 , q 2 , Qz, Q3) which we have solved numerically 
[23 lEH] • For q 3 = q 3 = the series solutions have a simple form and we are able to resum them. 
Thus, we have obtained analytical expressions for the three- Reggeon wave-functions with g 3 = 
()5.44)) - l)5.53p . For higher N we obtained a set of N — 2 non-separable differential equations 
which are hard to solve even numerically. We have presented explicitly such equations for 

iv = 4 

In order to solve the Reggeon problem for more than three particles one can make use of 
the more sophisticated technique, i.e. the Baxter Q— operator method [HHJ - It relies on the 
existence of the operator Q(w, u) depending on the pair of complex spectral parameters u and 
u. The Baxter Q— operator has to commute with itself (13.1)1 and with the conformal charges 
(13.2)1 . It also has to satisfy the Baxter equations (I3.3)1 - (I3.4I) . Furthermore, the Q— operator has 
well defined analytical properties, i.e. known pole structure (|3.12j) and asymptotic behaviour 
at infinity ()3.14|) . The above conditions fix the Q— operator uniquely and allow us to quantize 
the integrals g&. In turns out [38J that the Reggeon Hamiltonian can be rewritten in terms of 
Baxter Q— operator (I3.15j) . Therefore, combining together the solutions of the Baxter equations 
and the conditions for qk with the Schrodinger equation we can calculate the energy spectrum 
(J3.24)) . Moreover, we are able to determine the quasimomentum of the eigenstates (|3.20j) . i.e. 
the observable which defines the properties of the state with respect to the cyclic permutation 
symmetry (12.68)1 . 

One can solve the Baxter equations (|3.3Jl - (|3.4j) using the quasi-classical approach (|4.3|) . 
Following Refs. [21] we have presented the WKB approximation for the quantized values of qk 
(14.181) and explained the structure of their spectrum (|4.21l) . Applying this approximation, one 
introduces a small parameter 77 and performs expansion of the Baxter function for large values 
of qk- The single-valuedness requirement of the Baxter function gives conditions for q^ in a form 
of Bohr-Sommerfeld integrals (|4.15ft . Adding the condition coming from the possible values of 
the quasimomentum (|4.16l) one obtains the relations between the quantized values of qk (14.18)1 . 
However, not all quantization conditions can be obtained from (I4.18jl . In turns out that in the 
leading order of the WKB approximation for a given value of v h the quantized values of 
correspond to vertices of equilateral lattices l)4.21)h For lower conformal charges the conditions 
are underdetermined, i.e. in the N = 4 case one has only one of two required conditions (17.23ft . 

In order to obtain the exact values of qk we have used the method [2HEH] which makes it pos- 
sible to calculate numerically the spectrum of the conformal charges as well as other observables, 
i.e. the energy ()6.14|) and the quasimomentum (|6.24jl . The exact method consists in rewriting 
the Baxter equations and the other conditions for Q— operator eigenvalues as the A^— order dif- 
ferential equation l|fi.2j) . This equation can be solved similarly to the ^3— eigenequation (I5.1fi|) . 
Thus, one obtains the conditions for quantized qk (|6.27j) and the formulae for the energy (|6.14j) 
and the quasimomentum (|fi.24j) . 
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We have applied this exact method for states with N — 2, . . . , 8 particles. The quantum 
values of form continuous trajectories, Fig. 17.31 This fact is related to the dependence of 
q2 on the real parameter 1)5.2)1 . The energy is a continuous function of so the energy 
spectrum is gapless. This energy continuity leads to additional logarithmic dependence of the 
scattering amplitude 1)7.3)1 on the total energy s. The energy changes continuously along the 
trajectory 1)7.2)1 and reaches the lowest value in the vicinity of u h = 0. Moreover, since the qk 
spectrum possesses symmetries, 1)7.4)1 and 1)7.5)1 . the energy for some cases is degenerated. 

For N = 3 we have calculated the behaviour of the q 3 — spectrum for the conformal Lorentz 
spins rih = 0, 1, 2, 3 and the scaling dimension 1 + liv^. Some results for rih > were presented 
before in Ref. [22j for — 1 and Refs. [HD] for = 3. The quantized values of g^ 3 for 
given n h and fixed v h exhibit the WKB lattice structure ()4.21)) . which for N=3 takes a form 
an equilateral-triangle lattice (I7.ll)) . Figs. 17. 1117.51 The non-leading WKB corrections move the 
quantized value of (73 away from the lattice and cause that that the quantized values of q 3 lie 
outside a disk located around the origin of the lattice, i.e. near q 3 = 0. However, for odd n h 
there are states with q 3 = 0. They are called descendent states because their wave- function is 
effectively built of N = 2 Reggeon states. The non-descendent state with the lowest energy 

belongs to nn — sector and its energy is positive 1)7.15)1 . This state is double-degenerated and 

1/3 

it appears to be the nearest one to the origin of the q 3 lattice. However, the ground state for 
iV = 3 is the descendent one with n h = 1 and energy E 3 = 1)7.16)1 . Having found the exact 
values of q 3 we are able to calculate corrections to the WKB approximation (17.20)1 . Table l7T_l 
These corrections differ from the corrections obtained earlier in Ref. [21]. The difference seems 
to be caused by using only one expansion parameter 77 for two various conformal charges, qi 
and q 3 . The obtained corrections are subleading to the WKB approximation [21] which is an 
expansion for large values of conformal charges, i.e. 1 <C \q^ 2 \ ^ IqI^I- 

In the case with N = 4 particles we have constructed the spectrum for rih = depicting 
complicated interplay between the conformal charges: q 3 and q^. Such a complete analysis has 
been performed here for the first time. Earlier, for N = A full spectrum of q^ was shown in 
Ref. [2S] (however the corresponding q 3 spectrum was not discussed) and some values of q^ were 
found in Ref. [2Jj. The spectrum of q^ 4 has a structure of square-like lattice (17.24)1 . Fig. 17.71 In 
this case the spectrum of q 3 is more complicated. It turns out that it has a few possible forms. 
Firstly, there are simple states with q 3 = 0, Fig. 17.81 Moreover, we have found the q^ 2 — lattices 
whose distribution of vertices is similar to the distribution of vertices of g]^ 4 — lattice. We have 

__________ 1/2 

called this structures as resemblant lattices, Figs 17.1(1117.131 The next q 3 —lattices are called 
winding lattices. They wind around the origin and in course of this winding the distance 
between vertices increases and also the lattice goes away from the origin, Figs. I7.14ll7.16l The 
WKB approximation does not describe these structures (17.23)1 exactly because one quantization 
condition is missing. Finally, we have also considered the descendent states with g 4 = 0, Fig. 

____ -j^ /Q 

17.81 They have quasimomentum 84 = n and the structure of their q 3 —lattices is the same 

1/3 (fW™^H 

as for q 3 —lattices ()7.11)1 in the three-Reggeon case with 6* 4 = 0. For the N = 4 Reggeon 
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states with q 3 = we have calculated corrections to the WKB approximation 1)7.30)1 . Table [Ol 
Similarly to the N = 3 case they come from the fact of using only one expansion parameter 
i] for two different conformal charges, g 4 and q 2 . One has to notice that these corrections are 
subleading in comparison to the WKB limit |2 1 1 . 

For N > 4 we have more conformal charges so their spectrum is more complicated. There- 
fore, in these cases we have been concentrating only on the ground states, Table 17.51 Thus, 
we have shown the dependence of the ground state energy as a function of particle number 
N = 2, . . . ,8 (T7~38l) - (r7~37l) [22J, Fig. I7T71 When we enlarge number of particles N the energy 
seems to go to zero as 1/JV. For even N the energy of the ground state is negative so the 
contribution of the even ground states to the scattering amplitude increases with the total 
energy s. For odd N we have to consider separately the descendent states with q^ = and 
non-descendent ones. For the latter states the energy is positive so the contribution to the 
scattering amplitude decreases with s. For the states with q^ = the lowest energy En = 0. 
Therefore, the descendent states include the ground states for odd N. Additionally, we have 
calculated the energy along the ground state trajectories for N = 2, . . . , 8 as a function of u^, 
Fig. 17.181 These functions are symmetric in u^. The minimum energy is in = 0. The order 
of these energy levels tell us that the system for Vh = is anti-ferromagnetic one. Increasing 
v h we have observed reordering of the energy levels which makes the system for a large v h 
ferromagnetic. 

Next we have made comments on the approach of the Baxter Q— operator method presented 
in Ref. [2B1 EZ|- The authors of this approach solve the Baxter equation by using the pole 
ansatz (|8.1jl whereas the quantization conditions come from the normalization requirement and 
condition of equality of holomorphic energies 1)8.9)1 for different holomorphic solutions to the 
Baxter equation. We claim the the latter condition is too strong because the Reggeon system 
is two dimensional so the holomorphic energy is not a physical observable. 

At the end, we have passed to the deep inelastic scattering region where Reggeon states allow 
us to investigate properties of the structure function 1)9.3)1 of the scattered hadron [21]. Per- 
forming analytical continuation of the energy into complex we have calculated the anomalous 
dimensions as functions of the strong coupling constant 1)9.7)1 and the respective twists 1)9.21)) 
[26] which have been determined for each N separately. The case with N = 2 was discussed in 
Refs. [231 [2H]. It turns out that after the analytical continuation the energy for N > 2 is multi- 
valued function defined on the complicated Riemann surfaces and possesses branching points 
and poles [S3 EH EZ21 E6J, Fig. 19.21 The poles 1)9.30)1 occur on imaginary axis of Uh and their 
positions correspond to the twist n of the structure functions while expansion coefficients of the 
energy around the poles 1)9.18)1 are related to expansion coefficients of anomalous dimensions 
■y(a s ) in the strong coupling constant a s ()9.19)1 . Thus, we have found [2E] that for N Reggeon 
system the minimal twist n — N. However, for even N it is formed by the states belonging to 
the rih = sector while for odd N it corresponds to the — 1 sector. 

Thus, we have become convinced that although the Q— Baxter method is complicated it 
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is a powerful tool. It allows us to solve the reggeized gluon state problem for more than 
N = 3 particles. The above calculations are of interest not only for perturbative QCD but 
also to statistical physics as the SL(2, C) non-compact XXX Heisenberg spin magnet model 

(HUES EH). 

One has to mention that in this work we have not calculated all pieces of an elastic scattering 
amplitude but only the part that is responsible for the dependence on the total energy s (|7.3|) . 
which is determined by the intercept j = l — agE^/A. Similarly, in the deep inelastic scattering 
processes we have only identified exponents of the operator product expansion (OPE) series 
which are related to the twists and anomalous dimensions of the structure function F 2 (x, Q 2 ) 
()9.3|) . In order to find the full scattering amplitude one has to calculate the impact factors 
()9.9|) . i.e. overlaps between the Reggeon wave function and the wave functions of the scattered 
particles. These impact factors strongly depend on the scattered system and they are hard to 
calculate due to a large number of integrations. Therefore, with more than N = 2 Reggeons 
only overlaps of the simplest Reggeon states were calculated, i.e. in the r\ c — > J*(q) diffractive 
process [Hll EH EI], the double-diffractive production of J/ifj photon-photon processes 
[74 | l75l l76l 177] as well as the proton-proton scattering process (ZH] and other similar processes. 

One has to remember that apart from the GLLA method presented in this work there 
are other techniques which also make use of the reggeized gluons. Firstly, one may use the 
extended leading logarithm approximation (EGLLA) [H3 EHl EH] where the number of gluons 
in the t— channel may fluctuate. Moreover, Balitsky and Kovchegov [8QJ ED E2j derived a non- 
linear differential (BK) equation that for the two-Reggeon Pomeron contains the linear BFKL 
equation and additionally the non-linear terms which for higher energy s describe the saturation 
effect |82 [ 1831 IMl 170) . Furthermore, one can use the dipole model [H3 EEl E3 EH EH , the colour 
glass model [HH1 EDI EH E2] or the Wilson loop operators model (HBl El] in order to derive 
or explain diagrammaticly the BK equation. All these methods seem to be equivalent in the 
leading order to the BFKL. However, they include different non-leading corrections which cause 
that the different approaches better describe the scattering processes in different kinematical 
regions [34|. 

This work opens the way for further studies related to the multi-Reggeon states. Using the 
method presented in Refs. (2TJ E2 EH] it will be interesting to calculate the spectrum of the 
conformal charges for = 4 and 7^ and compare them to results obtained in Refs. [27, 28j. 
Additionally, increasing precision one my try to find the energy of the ground states for > 8 
and confirm the dependence En ~ 1/N. It is also of interest to perform analytical continuation 
of the energy for the states which have not been considered here and calculate their twist and 
anomalous dimensions. Moreover, we have mentioned that there are attempts to calculate the 
scattering processes with Reggeon exchange where a number of exchanging reggeized gluons 
varies in the t— channel, i.e. Extended GLLA method [HSlEEllZn]- These diagrams start to be 
important at higher energy s [34J and describe saturation which might allow to unitarize the 
scattering amplitude. To complete this picture one should also include the diagrams with N > 2 
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Reggeon states and the diagrams containing vertices of N Reggeons going into M Reggeons 
for any N and M. These diagrams would allow to calculate the DIS processes as well as the 
elastic scattering processes, i.e. onium-onium or proton-proton scattering. Furthermore, we 
can go beyond the multi-colour approximation and try to calculated NLO BKP diagrams. The 
next-to-leading order corrections to the BFKL kernel (NLO BFKL) was calculated by Fadin 
and Lipatov in Ref. [HE] and seems to play an important role for really high energy s. This 
calculation needs further studies. One has to notice the analysis of the multi-colour limit seems 
to be interesting for the anomalous dimensions of the structure functions in DIS processes and 
compare to Ref. (HOI - At the end we have to say that there is a connection between the GLLA 
approach and Af = 2 SUSY gauge theories presented in Ref. [97J. It is very impressive and this 
topic also needs more research. 
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Appendix A 

5X(2, C) invariants and other variables 



Let us consider a difference of coordinates (z\ — z 2 ). It changes under the SX(2,C) transfor- 
mation (|2.27j) as 

(4-4) = (ff^-fS) ={cz l + d)-\cz 2 + d)-\z 1 -z 2 ) . (A.l) 

One can see that during the transformation (jA.ljl the factors {cZ{ + d)" 1 appear in front of the 
difference. These factors also exist in the transformation law of the wave-function (I2.55j) . In 
order to cancel these factors one can construct a fraction where the same additional factors 
appear in the denominator and numerator of the constructed fraction variable. 
The fraction variable 

(Zi — Z 2 )(z 3 — Zq) 

x = r = (z 1 z 2 z 3 z ) (A.2) 

(Zl - Z )(z 3 - z 2 ) 

is invariant under the SL(2, C) transformations l|2.27)l . It is only one independent invariant for 
four coordinates (HHl [HH| - One can see that the fractions coming from the SL(2, C) transforma- 
tion simplified because: 



• the variable is a function of the coordinate differences (JAJJ), 

• the coordinates in the numerator and denominator of the variable x are the same. 

Therefore, simplifying the partial fractions we can obtain expression like ((azi + b)(czj + d) — 
(azj + b)(czi + d)) which with making use of ad — cb = 1 goes to (zi — Zj). As we can see, we 
have to build our invariants from differences of the coordinates. 

It is easy to see that performing permutations of coordinates we can construct six different 
dependent invariants 

Oiz 2 ^o) = x, {z z z 2 z x zq) = 1/x, 

(z^z^q) = 1/(1 -x), (z 1 z 3 z 2 z ) = 1-x, (A. 3) 

(z 3 z 1 z 2 z ) = (x-l)/x, (z^z^zo) = x/(x-l). 

Let us take another product of z%j 

W ' = iz< ( 4> z> ) = {CZ0 + d) \z ~z )(f z ) = (CZ ° + dfw ■ (A " 4) 
\. z \- z o)K z 2- z o) [zi - z ){z 2 - z ) 
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One can see that since in the denominator we have two additional z variables after the trans- 
formation we obtain a multiplying factor (czq + d) 2 . Similarly for w h we obtain 

w ' h = (cz + d) 2h w h . (A.5) 

Now one can compare this transformation to the transformation of the SL(2, C) wave-function 

(HMD . 

For N = 3 the variables w = -, — ( Z3 ~ Z2 ^ and x = ! Zl ~ Z2 w 23 ~ z °| transform under the cycling 
permutation (|2.68j) as 



(* - 1) M -> riij -> (* - ^ 

(— x) — > — yj — > > (— x) , 

u> — > u>(x — 1) — > it?(— x) — > u> 



(A.6) 



while under the mirror permutation: 

(-*) -> (ij -> (-^) 
u; — > w x — >■ u> 

For higher A we have more invariants. All of them can be constructed [HT| from 

N N N 



(A.7) 



^ (^r-l — z r){z r +l — Zo) 



Y[x r = (-l) N ; £(-l) p H ^ = °- ( A - 8 ) 



(Zr-l — Zo)(z r+ i — Z r ) . 
\ / \ / r= \ r= i j,=r+l 

Variables subject to the two conditions: 

f[x r = (-lf (A.9) 

r=l 

and 

AT N 

^(-l) r J] x fe = 0. (A.10) 

r=l fc=r+l 

From (IA.9j) we have ari = (-1)*/ IX=2 x r- From (jSHSj) we can calculate 



T-rAT 

llfc=3 Xfc 



X 2 = — ^ . (A.ll) 



Interchangeably, one can derive from (|A.9jl xjv = (—^) N /Yl^=i x r an d next from (|A.10jl 



(-1) 



N 



Xn ~ 1 ~ wv-i, .wTTiv-2 : • ( A - 12 ) 



E^r(-i) r nitv 

Thus, we can see that we have for N particles N — 2 independent invariants built of the particle 
coordinates zs. 
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Appendix B 

Solution of the eigenproblem 



B.l Case for N = 2 

For two particles, i.e. N = 2, we do not have any invariant variables Xj, (|A.8jl . Substituting 
(15. lj) into the equation 1)5. fl) we obtain 

fe^io)* 10 ^)^^)* 13 = (~h(h - 1) + Sl ( Sl - 1) + s 2 (s 2 - 1)) (zio)* 10 (^)^(si 2 )* ia . (B.l) 
Next, solving (|B.ljl for fc^- we get the eigenstates: 

2. ^10,^20) = (^io)- m - Sl+S2 (^o)" 1+/l+sl ^ 2 (^i2) 1 -' l - Sl - S2 . 

The first solution transforms with the proper factor (cz + d) 2h (czi + d) 2si (cz 2 + d) 2s2 . The 
second one has h — > 1 — h. 

A full eigenstate of q 2 with N = 2 can be presented as a linear combination of the difference- 
coordinate products. It comes from the fact that during the SL(2, C) transformation the factor 
(cZi + d) 2s have to be obtained 1)2. 5 5|) . Thus, the eigenstate, up to some normalization factor, 
has a form [15)13)116] 

h / \ si / \ S2 



*(Z 10 , ^ 20 ) = (^ 1 o)-' l - Sl+S2 (^o)- /l+S1 - S2 (^12) ?i - Sl - S2 =~ ' ZU 



Z 20 Z W 



\ ( z 20 V 1 / gig V 



(B.2) 



For the homogeneous model, i.e. s = s. 



^(z 10 , Z 20 ) = (Z 1 0)- ,1 (^20)- /I (^12) ,l - 2S = f-^-) • (B.3) 

\ ^10^20/ 

In the special case for s = it leads to 

h 



^(zio, ^20) = M-^-'W* = ( — ) (B.4) 



-Y 



while for s = 1 it gives 



^12 \*10*20/ Z 1Q Z 2D \*10*20/ 
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B.2 Case for N = 3 



In this case, the solutions with z can be reduced to one independent function 

^(*10,Z20,*30) = (z 10 )- /l + fc 3 2 -«i+^+^3( 2l2 )fci 2 ( Z20 )-fc3 2 -^-2s 2(z3o) -/ l+fcl2+Sl+S2 - S3 
X (z 3 l) h ~ k32 ~ kl2 ~ ai ~ S2 ~ S3 {z 32 ) k32 = 
= (^ 10 )- /l - Sl+S2+S3 (^0)- 2S2 (^0)- /l+Sl+S2 - S3 (^3l) h - S1 - S2 - S3 {^) kl2 (^) fcM 

(B.6) 

and to the other one with h — ► 1 — h. It is invariant under the cyclic permutation and contains 
two free parameters, k\ 2 and k 32 . 

For N = 3, there is one independent invariant variable (|A.2j) . e.g. x = x 2 . Next, we 
multiply our wave-function by an arbitrary function of x, F(x). There is only one pattern of 
the eigenfunction which satisfy the eigenequation for q 2 and the SL(2, C) transformation 1)2. 5 5j) 
law. This eigenfunction looks like 

^(z 10 ,z 20 ,z 3 o) = (^io)-' 1 - Sl+S2+S3 (^o)- 2s2 (^3o)" ?i+Sl+S2 - S3 (^3i)^ Sl " S2 - S3 i 71 (a:) = 

ft— Sl— S2— S3 (B.7) 

Fix). 



(21O) 2S 1(02O) 2S2 (Z3O) 2S 3 V 210230 

For the homogeneous model (s = Sj) we have an ansatz 
*(z 10 ,z 20 ,z 3Q ) = (z 10 )- h+s (z 20 r 2s (z 30 r h+s (z 31 ) h ~ 3s F(x) = (-^-) %(*) 



^iqZ 2 qZ 3 q) 2s V ^10^30/ 



(B. 



Substituting s = we have 



*(*io, Z20, z 3Q ) = (z 10 )- h (z 30 )' h (z 31 ) h F(x) = ( ) F(x) , (B.9) 



whereas for s = 1, it is 



- h+1 (z, n )- 2 (z,n)~ h+1 (z^ h - S F(x) = — ( " 



*(zio, ^o, z 30 ) = (z 10 )- h +\z 20 )-\z 30 )- h+ \z 31 )^F(x) = ^ -^L_ F(x) 

{ZioZ 2 oZ 3 q) \Zi Z 30 I 



(B.10) 



B.3 Case for TV = 4 

Here, like in the previous cases, we have also four 2x2 solutions. When we take solutions with 
Zo they reduce into two functions: 

*(Z 10 , £20,-^30, *40) = {z w )- k ^- k ^- k ^~ 2s ^Z 12 ) k ^{z 2Q )-^^ 

X ( K Zl4 : ) kli {z 3 l) kjl (Z 32 )~ h ~ fc 24-fcl2-fe3l-fc34-fcl4-Sl-S2-S3-«4 

(B.ll) 

and the other one is with h — > 1 — /t. The function fjB.l 1|) is invariant under cyclic permutation 
and contains two free parameters. 
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In the N = A case we have N — 2 = 2 independent invariant variables. Let us take 
X\ = r^zTTT^zf^ an d £2 = ^ 1 ~ Z2 ^ 23 ~ Z °| . We multiply our wave- function by an arbitrary 

\Z4. Zq )\Z2 Z\ ) yZ\ Zq ) \Z% Z2 ) 

function of these variables, F(x 1 ,x 2 ). Considering f)2.55() . we obtain a solution for the q 2 - 
eigenequation 



* (ZIO, Z20, Z 30 , ^40) = (^ 10 )- fc «-fc3i-fei4-2 Sl ^ i2 )fc 12 ( /220 )-/ t +fe3i+fe 3 4+fci4+.i- S2 + S 3+ S 4( 224 > 



mo; 

x(2 1 4) fcl4 (2 31 ) fc31 (23 2 ) _/t_fc24_fcl2 " fc31_fc34_A:i4_Sl_S2_S3_S4 /(a;i,a;2) 

h— Si— S2— S3 — S4 

(Zl()) 2s l (^20) 2S2 (^3C)) 2s 3 (Z4()) 2s 4 I Z 10 V 30 



/ \ n— SI— S2 — S3— S4 

( ) F( Xl , x 2 ) 

\ Z10230 / \ 11 



(B.12) 



For the homogeneous model (s = Sj) we have an ansatz 

/ \ h— 4s 

(^10^20^30^40) 2 

For s = 



*(zio, ^20, ^30, ^4o) = 7 ^7 ( ^ 31 ) F(xx, x 2 ) . (B.13) 

\ £10^30 / 



^(^10,^20,^30,^40) = (^10) h {z3o) h (z 3 i) h F(x 1 , x 2 ) = ( 2:31 ) F(xi,x 2 ) 

\^10^30/ 

and for s = 1 

*(*!(>, *20, *30, *4o) = (2l0)^ +2 (220) _2 (^30)^ +2 (^40) _2 (23l)^ 4 i 71 (xi,X 2 ) = 



(B.14) 



(■210^20 ^30^40; 



(—) 

\^10^30/ 



F(xi,x 2 ). (B.15) 
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Appendix C 



Solutions for N = 3 



C.l Solutions for 5 = around x = + 
C.l.l Solutions for q% ^ and h Z 

The eigenequation for q 3 is a differential equation of the third order so around each 

singular point, x = 0,1, oo, it has three independent solutions. Around x = + we have an 
indicial equation 

(h-n-r)(r + n — l)(n + r) = 0, (C.l) 

so its solutions are given by T\ = h, r% = 1 and r 3 = 0. As we can see we have two cases when 
h £ Z (one solution with Log(a?)) and h G Z (one solution with Log(x) and one solution with 
Log 2 (x)). As we will see below we have to consider separately solutions with q 3 = 0. In the 
first case 

m 2 (x) = x r2 ££L a nir2 x n , (C.2) 
113(2;) = x r3 Xir=o & n,r 3 ^ n + a; r2 E^=o a„, r2 x n Log(x) , 
where ao, r is arbitrary (e.g. equal to 1) 

(iq 3 - (h - 2r)(h - r)r) 

a>i,r — — 77 — 1 rrr, — \ — a o,r K^- 6 ) 

[h — 1 — r)r(l + r) 

and m = n + r 

(h-m+ l)(h -m+ 2)(m - 2; 



(/i — m){m — l)m 



whereas &o,r 3 = • Qo,r 2 an d &i,r 3 is arbitrary. One can notice that coefficient &o,r 3 is well 



&n-2,r + 

(^3-(l + fe- m )( m -l)(fe-2(m-l))) 

H 77 w 77; a n-l,r , 

(n — m){m — ljm 



defined only for g 3 7^ 0. Moreover, 

, _ 2+(h-3)h-iq 3 , 3h-8 _ , 6+fe(fe-6) /p r\ 

"2,r 3 — 2(2-/i) U1 > r 3 2(2-h) U1 > r 3 "T 2(2-h) °> r 3 ^ -t V 
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and 



7 _ (h+l-m)(h-m+2)(m-2) i . iq- A -{l+h-m)(m-l){h-2(m-l)) 7 . 

°n,r 3 — (h-m)(m-l)m °n-2,r s + (7i-m)(m-l)m °n-l,r 3 + 

+ (h-m){m-\)m a n-3,r 3 (h-m)(m-l)m a n-2,r 3 + 

2m-3m 2 +h(2m-l) 

(ft— m)(m— l)m re— l,r3 > 



where m = r 3 + n. 



C.1.2 Solution with q 3 = 



In this case a , r defined in (|C.2jl is arbitrary (e.g. equal to 1). For the first solution with r± — h 
we have a^ ri = whereas for the third one a 1>rs is arbitrary (let us take 0). It turns out that 
we do not need Log— solutions. The second solution is more complicated. One can derive an 
exact formula for 

nk-h_ T(l-h + n) 
=1 FT! ~ r(i-fc)r(n + 2) ao - ra (C - 7) 

and performing summations 

1l\ {x) X n Q"n,,r\X X Oo,ri j 

m 2 (x) = x r2 ££1 a n , r2 x n = a , r2 x J2n=o r(i-fc)r(n+2) xH = - a o,r 2 i((l - ^ - 1) , ( c - 8 ) 

^3(-^) X 3 /~] T! _q d nr;K X ^0^3 • 

Gathering together this all solutions we have 

u(x) = A + B{-x) h + C{x - l) h (C.9) 

where A, B, C are arbitrary. The above solution was presented by Lipatov and Vacca in Refs 
. |5IIESj. 

C.1.3 Solution with q% ^ 0, q<i = and /i = 1 

For h = 0, i.e. g 2 = 0, and g 3 ^ we have a different set of solutions. Here we have three 
solutions of the indicial equation (|C.1J) which are integer r± = 1, r 2 = 1, r 3 = so the solutions 
are given by 

^l{x) X *y^ Jn _ n Q-n.ri *E : 

U2 (x) = X r2 ^r=0 ^n.ra^" + ^ E^O an,n^™L0g(x) , (CIO) 
U 3 (x) = Of 3 ^r=0 Cn.ra^™ + ^ J2n=0 b n ,r 2 X n Log(x) + X ri Yln=0 a n ^ ri X n hog 2 (x) , 



where ao, r is arbitrary (e.g. equal 1) 



r + r 2 (2r-3) 

Ql.r — ^77^ — ; — \ : — a 0,r l < ^- li J 

r^(l + r) — iqs 



and m = n + r 



(m-3)(m-2) 2 (m - 2)(m - l)(2m - 3) - sg 3 ^ 10 , 

On.r — TTo 0"n-2,r H 7 7TT5 dn-l.r > l ( ^--'-^J 

[m — l) z m [m — l) z m 
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whereas 



h,r 2 = 7^(-mb ,r 2 + «o,n - 5ai jn ) (C.13) 



, (n-l) 2 (n-2) , . (2m 3 -3m 2 +m-iq s ) h . l+6n(n-l) , 

°«.' - 2 n 2 (n+l) °n-2,r 2 ~T n 2 (n+l) w »i-l,r 2 T n 2 (n+l) a n-l,)-i"t" ,~ , \ 

(n-l)(5-3n) (2+3n) 



+ n 2 (n+l) a "-2,ri „(„+l) a n,»*i 

while Co jr . 3 = ^~Qo,ris c i,r 3 is arbitrary and 



c 2,r 3 — — n*93 c i,7-3 + bo, r2 ~ 5fci,r 2 + 3a 0jri — 4ai )ri , (C.15) 



_ (n-2) 2 (re-3) . -ig 3 -6+n(13+n(2n-9)) __ 2(n-2)(3n-8) . , 

C n ,r 3 — n(n-l) 2 C ™-2,r 3 i" n(n-l) 2 C n _i jr3 „(„_l)2 »n-2,r 2 + 

I 26-36n+12n 2 i 6n 2 -8n+2 i , 2(7-3n) 6(3-2n) (T 1 fil 

* 1 " n(n-l) 2 °n-l,r 2 n(n-l) 2 U »V2 „(„-l)2 «n-3,n „(„_ 1 )2 «n-2,n "I" 

I 2(2-3n) 
T n(n-l) 2 " , '- 1 . r i ' 

C.1.4 Solution with q2 = qz = and /i = 1 

In this case a 0) r is arbitrary (e.g. equal to 1). For the first solution with r± = (h = 1) we have 
a l ri = and for the third one ai jr3 with r 3 = is arbitrary (let us take 0). The second solution 
is more complicated. We need Log— solutions 



u 2 {x) = x 7 " 2 yj b n , r2 x n + x ri a„ in a; n Log(a;) . (CI 7) 

n=0 n=0 

Using recurrence relations with b± jr2 = |ao,n (and a ntn = for n > 0) 

b 2 ,r 2 = 3&l,r 2 , , , 

(n-l) 2 (n-2) , (n-l)n(2n-l) , 
"n,r 2 n 2 (n+l) w "-2,r 2 T n 2( n+1 ) "n-l,r 2 

one can derive an exact formula for b n r2 = 6 l r2 and performing summations with arbitrary 

& ,r 2 (= 0) we have x J^^Lo b n ,r 2 x n = 26 l r2 (x + Log(l — x) — xLog(l — x)), so that 

m 2 (x) = a 0) n (sc + Log(l — x) — xLog(l — x) + xLog(x)) , (C.19) 

Gathering together all this solutions we have 

u{x) = A + B{-x) + C((x - l)Log(x - 1) + (-x)Log(-x)) , (C.20) 
where A, B, C are arbitrary. See also Ref. [2j 
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C.1.5 Solution with q% ^ 0, = and h = 

For /i = 0, i.e. q 2 = 0, and arbitrary g 3 ^ we have three solutions 

M2 ( x ) = a^X^o 6n,raa: n + a: ri Z)^=o a n,na; n Log(x) (C.21) 

U 3 (x) = X rs J2n=0 C n,r 3 X n + 2x"' 2 J2n=0 b n , n X n Log(x) + X Tl J2n=0 On.n^Log 2 (x) 

where r 1 — 1, r 2 = and r 3 = 0. Here a ,r is arbitrary (e.g. 1) 

_ -iq 3 + 2r 2 

ai ' r _ r (l + r) 2 a °' r [ ' 

and m = n + r 

_ _ {m-2f -iq 3 + 2{m-lf 

u n,r — n u n—2,r ' , 1 \ t u n—l,r ; l^.z,OI 

rrr (m — l)m z 

whereas &o,r 2 — ~j^ a Q,r\ and b\ >r2 is arbitrary while 

&2,r 2 = ^((2 - iq 3 )bi,r 2 + 6a ,n - 8ai, ri ) , (C.24) 



, (n-2) 2 , -ig 3 +2(»-l) 3 , _ (n-2)(3n-4) , 

"n,r 2 — n 2 (%-2,r 2 "T „2( n _!) "n-l,r 2 n 2 (n-l) u "-3,ri-|- 



I 6(»-l) , 2-3n 

n 2 u "-2,ri ~r n 2( n _!) "n-l,n 



The coefficient Co ir3 = -^-(2ao, n + &i, r2 ) and arbitrary and 



(C.25) 



l 5 

C2,r 3 = ^(2 - iqz)ci,r 3 + 3&i, r2 - 46 2 ,r 2 + 3a , ri - -ai,n , (C.26) 

(n-2) 2 , -ig 3 +2(n-l) 3 _ 2(n-2)(3n-4) , , 

t-n,r 3 — n 2 (-n-2,r 3 1~ n 2( n _i) ^n-l,r 3 n 2 (n-l) u n-2,r 2 T 

I 12(n-l) , | 2(2-3n) , 2(5-3») , 12 , (T 97^1 

H ^ 0n-l,r 2 + n (n-l) <V 2 + „2( n _i) «n-3,n + ^ fl n-2,n+ l^- z 'J 

I 2(l-3n) 
r n 2 (n-l) U "- 1 .'-l • 

C.1.6 Solution with h = and ^2 = 93 = 

In this case a 0) r is arbitrary. For the first solution with r± — 1 we have recurrence relations 
with ai jn = |a , ri 



(C28) 



(22,ri 3^1, ri 

— (n-1) 2 _ 2ni 

a n,ri — n+1 Qn-2,n + n+1 a n-l,n • 



Mi (a;) = a; ri ^a njri x n = z z n a 0|ri j = -Log(l - x) . (C.29) 



Summing series up we derive an exact formula for 

n=0 n=0 

In the second solution, with r 2 = 0, we have arbitrary a^ r2 (let us take 0) and solution is 

U2 (x) = ao,r 2 - The third one is with Log— solutions 

00 00 
u 3 (x) = x ra ^ h,r 3 x n + x r ' 2 ^2 a n>r2 x n Log(x) . (C.30) 

n=0 n=0 
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Using recurrence relations with arbitrary 6 ,r 3 

^2,r 3 — \^l,rz 



h ("-2)^ | 2(n-l) 2 , ( C - 31 ) 

"n,r 3 n 2 "n— 2,r 3 T n 2 "n— l,r 3 



one can derive an exact formula for b nt r 3 = &i,r 3 ^ and performing summations with arbitrary 
ko,r 3 (= 0) we have Yin 

,r 3 Log(l — x). All these solutions look like 

Ui(x) = X ri Yln=0 a n,r 1 X n = -ai )f . 3 Log(l - x) , 

m 2 (x) = x r2 Xl^Lo a n,r 2 x n = «0,r 2 , (C.32) 

U 3 (x) = X r3 ^7=0 b n,rs xn + a ,r 2 Log(x) = -& ,r 3 Log(l - x) + a ,r 2 Log(x) . 

Gathering together all this solutions we obtain 

u(x) = A + BLog(-x) + CLog(x - 1) , (C.33) 
where A, B, C are arbitrary. See also Ref. [27]. 



C.1.7 Solution with q$ ^ and h = 2 

For h — 2 and ^3 7^ we have three solutions 

u 2 (x) = x r2 ^~ =0 & n,r 2 ^ n + a^ 1 Er=o an,n;£ n Log(£) , (C.34) 

U 3 (x) = X rz Y,n=0 C n,r 3 X n + 2x T2 Y,n=0 &n,r 2 ^ n Log(x) + X Tl Y,n=0 a^r^Log 2 (x) 

with ri = 2, r 2 = 1, r 3 = 0. Here ao,n is arbitrary (e.g. 1) 

-iq 3 + 2r(r-l)(r-2) 

a l,r 1 — T^) 7\ a 0,ri (^.OOJ 

r{r z — 1) 



and m = n + r 1 



(m-3)(m-4) -ig 3 + 2(m - l)(m - 2)(m - 3) 

-«n-2,ri H / ^77 77 <2n-l,ri , (O.ODJ 



jl m(m — 1) (m — 2)(m — l)m 

whereas &o,r 2 = — if~ a o,n and b\^ 2 is arbitrary while 

&2,r 2 = ^(-«<?3&l,r2 + 4a 0,n - Hai,n) > (C.37) 




(C.38) 



, _ ( n - 2 )(n-3) 7 , -iq 3 +2(n-2)(n-l)n , ll+3(n-4)n 

Wn,r 2 — n(n+l) u n-2,r 2 ' n(n 2 -l) u n-l,r 2 n (n 2 -l) "n-3,ri-r 

6(n-2)n+4 l-3n 2 
" r n(n 2 -l) u n-2,ri "T „(„2_i) u n-l,ri • 

Moreover, Co >r3 = j|-&o,r 2 and Ci jr3 = _ 1^2£g^2ig_ _ 6 °°^ ri and C2 >r3 is arbitrary while 

c 3 ,r 3 = g(-ig3C 2 ,r 3 + 2& o,r 2 + 8&i, r2 - 226 2 ,r 2 + 12a , n - 12a hri ) , (C.39) 



_ (n-3)(n-4) -«, 3 +2(n-3)(n-2)(n-l) 2(3n(6-n)-26) , , 

t-n.rs — n(n-l) c "-2,r 3 T n(n-l)(n-2) ^n-l,r 3 i* „(„_i)( n _ 2 ) u n-3,r 2 -^ 

, 4(ll-12m+3m 2 ) i __ 2(2-6to+3to 2 ) , _ 6(n-3) . 12 , 

"T" n(n-l)(n-2) °™-2,r 2 n(n-l)(n-2) °n-l,r 2 n(n-l)(n-2) a "-4,ri + n ( n _i) a n-3,ri T 

6 - 
n(n-2)" n -2^i ■ 



(C.40) 
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C.1.8 Solution with q$ = and h = 2 

In this case a ,r is arbitrary (e.g. equal to 1). For the first solution with r 1 = 2 we have 
remaining coefficients a n> o^ ri = 0. The second solution, with r 2 = 1, has an arbitrary a^ r2 (we 
can take it as 0), ai )T2 — an d third one with arbitrary ai >rs , <i2,r 3 (we also set them to 0). All 
these solutions look like 

U\{x) X ^En=0 Q"n,ri% % ^0,r 3 j 

U 2 (x) = X T2 J2n=0 a n,r 2 X n = Xa ,r 2 , ( C -41) 
u 3( x ) x 3 En=o" B .' , 3''' ^0,r 3 • 

Gathering together this all solutions we have 

u(x) = A + B(-x) 2 + C{x-1) 2 , (C.42) 

where A, B, C are arbitrary. As we can see this solution corresponds to the solution with q 3 = 
and arbitrary h {0, 1}. For other integer h # {0, 1} we can observe the same correspondence. 



C.2 Solutions for 5 = around x = 1 
C.2.1 Solutions for 93 7^ and h Z 

Similarly to the case x = for x = 1~ we have independent solutions where an indicial equation 
has a form 

(h-n-r)(r + n- l)(n + r) = 0. (C.43) 

Its solutions have following values: r± — h, r 2 — 1 and r 3 = 0. As we can see we have two cases 
when h £ Z (one solution with Log(x)) and h G Z (one solution with Log(x) and one solution 
with Log 2 (x)). We will see below that we have to consider solutions with q 3 = 0, separately. 
In this case, for g 3 7^ 0: 

ui(x) = (l-^E^oVil 1 - 1 )"- 

u 2 (x) = (l-xrE^^fl-^- (C44) 

« s (x) = (i-xrEZ bn,r 3 a-xr + (i-xrj:n=o^,r 2 (i-xrLog(i-x), 

where a , r is arbitrary (e.g. 1) and 

_ (ig 3 + (h-2r)(fr-r)r) 

Ol,r — 77 ; \ /, . \ a 0,r 

(h — 1 — r)r(l + r) 

and m = n + r 

(/i - m + l){h - m + 2)(m - 2) (?g 3 + (1 + /i - m)(m - l)(h - 2(m - 1))) 

0"n,r — 77 77 i\ a n-2,r 77 77 77 Q>n-l,r j 

[h — m)[m — l)m (h — m)(m — l)m 

(C.46) 

whereas 6 ,r 3 = Qo,r 2 ( one can notice that coefficient 6 ,r 3 is valid only for g 3 7^ 0) and &i jr3 
is arbitrary 

_ »g 3 +(fe-2)(fe-l) r _ 8-3fe - , 6+fc(fe-6) 4? n 
"2,r 3 — 2(2-/i) W1 > r 3 2(2-h) U1 > r 3 "r" 2(2-/i) "O-^ V^-^V 
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and 

7 _ (h+l-m)(h-m+2)(m-2) i __ iq 3 +(l+h— m)(ro-l)(ft-2(m-l)) 7 . 

°".'"3 — (h-m)(m-l)m °n-2,r 3 (h-m)(m-l)m °n-l,r 3 + 

, fc 2 +fe(7-4m)+(m-2)(3m-4) __ /i 2 ~6fe(m-l)+6(m-l) 2 . /p aq\ 

+ (ft-m)(m-l)m "n-3,r 3 (h-m)(m-l)m a n-2,r- A + \yi.10) 

2m-3m 2 +h(2m-l) 

(/i— m)(m— l)m n ~li r '3 

where m = r 3 + n. 

We can easy see that above coefficients correspond to coefficients for solution around x = 0~ 
with g 3 -> -53. 



C.2.2 Solution with g 3 = 

Here ao, r is arbitrary (e.g. equal to 1). For the first solution with r% = h we have ai )n = and 
for the third one a^ r3 is arbitrary (let us take 0) . We do not need Log solutions. The second 
solution is more complicated. One can derive an exact formula for 

nk — h T(l — h + n) 
, _J+i = r(i-/or(n + 2) a °- ra ( } 

and performing summations 

Ui(x) = (1 - x) n Y^Lo a n,n(l ~x) n = (1 -x) h a 0jri , 



U 2 {X) = (1 - X) r2 J2n=0 a n,r 2 (l ~ x) n = fl ,r 2 (l ~ x) J2n=0 r(l-h)r(n+2) ^ ~ ^ 

= — ao,r 2 ^( a; — 1) ) 
u 3 (a;) = (1 - a;) r3 X^=o a n,r 3 (l - ^) n = a ,r 



(C.50) 



Gathering together all this solutions we have 

u{x) = A + B(-x) h + C(x - 1)\ (C.51) 
where A, B, C are arbitrary. The above solution was presented by Lipatov and Vacca in Refs. 

|531 EH- 

C.2.3 Solution with q% ^ 0, q<i = and h = 1 

For h = 1, i.e. g2 = 0, and q 3 ^ we have three solutions to the indicial equation (|C.43jl which 
are integer n = 1, r 2 = 1, r 3 = so solutions are 



(C.52) 



Mil 




= (1 


— X 




En=0 a n,r 1 (1 


— X 


r 


j 


U 2 i 




= (1 


— X 






— x] 


,n 


+ ^ E^O a n,n(l - x) n L0g(l - x) , 






= (1 


— X 


Y 3 


X/n=0 C n,r 3 (1 


— £ 


,n 


+ 2(1 - XY 2 J2n=0 - x)"L0g(l - X 






+ ( 




x) 


Z^n=0 "n,ri I 


1 - 




)"Log 2 (l-x), 



where ao, r is arbitrary (e.g. equal to 1] 



_ iq 3 + r + r 2 (2r -3) 
r z (l + r) 
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and m = n + r 



(m-3)(m-2) 2 (iq 3 + (m - 2)(m - l)(2m - 3)) 

0"n,r — / 7T5 On-2,r H 7 TT^ Ora-l,r j (Lv.54j 

(m — lj z m (m — lj z m 



whereas 



6l,ra = ^(m b 0,r 2 + Oo.ri - 5ai,n) , (C.55) 



, _ (n-l) 2 (n-2) , (iq 3 +2n 3 -3n 2 +n) , l+6n(n-l) 

"n,r 2 — n 2 (n+l) u n-2,r 2 < n 2 {n+l) u n-l,r 2 ' n 2 (n+l) u n-l,»*iT 

(n-l)(5-3n) _ (2+3n) 

~r n 2 (n+l) ""-2,ri „(„+l)Vi ■ 



(C.56) 



Moreover, Co, r3 = — Oo,n c i,r 3 is arbitrary and 
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c 2,r 3 — 2*?3 c i,r 3 + ^o,r 2 — 5&i jr2 + 3a 0i ri — 4ai jn , (C.57) 

_ (n-2) 2 (n-3) . ig 3 -6+n(13+n(2n-9)) _ 2(n-2)(3n-8) i . 

Cn,r 3 — n(n-l) 2 C n-2,r 3 + n(n-l) 2 C„-l,r 3 n ( n -l) 2 °n-2,r 2 + 

2(13+6(n-3)n) , _ 2(l+n(3w-4) , 2(7-3n) _ 6(3-2n) ,m 

^ n(n-l) 2 u "-l,»*2 n(n-l) 2 U ™> r 2 ^ ra(ri-l) 2 u ™-3,ri „( n _l)2 "n-2,n T ^.uoj 
I 2(2-3») 
^nCn-l) 2 "™- 1 . 1 -! • 

C.2.4 Solution with ^2 = 93 = and /i = 1 

In this case a , r is arbitrary (e.g. 1). For the first solution with r\ = (h = 1) we have a 1>ri = 
and for the third one a^ rs with r 3 = is arbitrary (let us take 0). The second solution is more 
complicated. We need Log— solutions 

oo oo 

u 2 (x) = (1 - x) r \J2 (1 ~x) n + x r \J2 a «.n ( X " ^"Logl 1 - • (C59) 

n=0 n=0 

Using recursive relations with b\ >r2 = |a , n (and a„ in = for n > 0) 

^2,r 2 = 3^1, r 2 5 

, . (n-l) 3 (n-2) , (n-l)n(2n-l) , ( C - 60 ) 

"n,r 2 n 2 (n+l) u ™-2,r 2 T n 2( n+1 ) "n-l,r 2 

one can derive an exact formula for 6 ra r2 = bi >r2 ^ n ^ n and performing summations with arbitrary 

&o,r 2 (= 0) we have 

oo 

xJ2 bn,r 2 (l - x) n = 26i >r2 ((l - x) + Log(:r) - (1 - x)Log(x)) (CM) 

n=0 

what gives us 

ui(x) = (1 - x) ri Y^7=o a n,n (1 - a;) n = ao,n = (1 - a?)ao,n , 

u 2 (x) = a , ri ((l - x) + Log(x) - (1 - x)Log(x) + (1 - x)Log(l - x)) , (C.62) 
u 3 (x) = (1 - x) r:i Er=o a ",r 3 (l - x) n = a , r3 . 
Gathering together all this solutions we have 

u(x) =A + B(-x) + C{(x - l)Log(x - 1) + (-x)Log(-x)) , (C.63) 

where A, B, C are arbitrary. Here we have also obtained the consistent solution. 



122 



C.2.5 Solution with q% ^ 0, q<i = and h = 

For /i = 0, i.e. q2 = 0, and q^ ^ we have three solutions 



«2(x) = (i-^E^^ii-xr + fi-irE^vfi-wi-x), 

U 3 (X) = (1 - x)^ £~ Q Cnir3 (l - X y + 2(1 - X)** En=0 Kr 2 (± ~ x) ra Log(l - X ) + 

+ (1 - x)^ E"=o V(! - ^) n Log 2 (l - x) , 
where coefficients ri = 1, r 2 = 0, r 3 = are integer while ao, r is arbitrary (e.g. 1) 

m + 2r 2 

ai ' r ~~ r (i + r )2 a °- r ^• bb ' ) 

and m = n + r 

(™ ~ 2) 2 ig 3 + 2(m - l) 3 

(m — l)m 2 

whereas &o,r 2 = 7~ a o,ri arid 6i r2 is arbitrary while 



h,r 2 = J((2 + i93)6i,r2 + 6a o,n - 8ai, ri ) , (C.67) 



, _ (n-2) 2 , iq 3 +2(n-l) 3 , _ (»-2)(3n-4) 

"n,r 2 — n 2 "n-2,r 2 T n 2( n _x) "n-l,r 2 n 2 (n-l) u n-3,riT 

1 6(n-l) 2-3n - 

n 2 "™-2,ri ~r „(„_!) "n-l,n • 

Moreover, Co jr3 = ^-(2ao, n + &i,r 2 ), Ci,r 3 is arbitrary and 



(C.68) 



1 5 

C2,r 3 = 7(2 + iqz)ci,r 3 + 3&i, r2 - 4& 2 ,r 2 + 3a , ri - -ai, n , (C.69) 



(n-2) 2 j g3+2 ( n -l) 3 2(n-2)(3n-4) , 

<-ri,r 3 — n2 (-n-2,r 3 T „2( n _ 1 ) ^n-l,r 3 n 2 (n-l) u n-2,r 2 ^~ 

12(n-l), , 2(2-371), 1 2(5-3w)„ ,12, 



i^n-JJ, ^2-3wj , , 2(b-dn) 12 , 7f) \ 

"I" n2 "n-1,7-2 "T n (n-l) W ™,r 2 n 2 (n-l) u ™-3,ri T „ 2 «n-2,ri"T ^.iuj 



I 2(l-3n) 
r n 2 (n-l) U "- 1 .''l • 



C.2.6 Solution with 92 = 93 = and h = 

In this case a ,r is arbitrary (e.g. 1). For the first solution with r\ — 1 we have recurrence 
relations with ai ;n = §ao,n 



° 2 ' ri ~ 5 a Vi' (C.71) 

(n-1) 2 , 2nf_ v ' 

u n,ri n+ i u n-2,r\ "T n+1 "n-l,7-i • 



Summing series up we derive an exact formula for 

00 00 

U1 (X) = (1 - x)* 1 ^,(1 " X) n = (1~X) ]T(1 - x)^ — — = -Log(x) . (C.72) 



rt=0 n=0 
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In the second solution, with r 2 = 0, we have an arbitrary a^ r2 (let us take 0) so the solution 
is u 2 (x) = a , r2 . The third one is with Log solutions 

oo oo 

u 3 (x) = (1 - x) r * " *) n + (1 - X T Yl a "^(l - ^) n Log(l - x) . (C.73) 

n=0 n=0 

Using recurrence relations with arbitrary &o,r 3 

b 2 ,r 3 = 2^3, 



, (n-2) 3 , 2(n-l) 2 , 

u n,rs n 2 u n—2,r 3 ~T~ „2 "re— l,r3 



one can derive an exact formula for 6 n r . 3 = 6i, r3 r and performing summations with arbitrary 
6 ,r 3 (= 0)we have Y^=o^n,r 3 (^ ~ x ) n — ~b 0rs Log(x). All these solutions look like 

ui{x) = (1 - x) ri Er=o a «,n( 1 - X T = -ai,r 3 Log(x) , 

M 2 (x) = (1 - a;) r2 $^^=0°n,ra(l ~ X T = a 0,r 2 , 

U 3 (x) = (1 - X) r:i J2n=0 b n.r 3 (l ~ x) n + a ,r 2 Log(l - x) = -fe ,r 3 Log(x) + a ,r 2 Log(l - x) . 

(C.75) 

Gathering together all this solutions we have 

u{x) = A + BLog(-x) + CLog(x - 1) , (C.76) 
where A, B, C are arbitrary. 



C.2.7 Solution with <? 3 ^ and h = 2 

For integer ft, = 2 and 7^ we have three solutions to the indicial equation (IC.43j) which are 
integer rx = 2, r 2 = 1, r 3 — so solutions are 

ui(x) = (1 -x) ri Enlo^nC 1 -a;)", 

m 2 (x) = (l-x) r2 ^ =0 6 nir2 (l-x)" + (l-x) ri ^ =0 a ri>r . 1 (l-x) n Log(l-x), 
u 3 (z) = (l-xrE"=o c ^ 3 (l-x)" + 2(l-x) r2 Er=oVr- 2 (l-^) ri Log(l-x) + 
+(1 - x) n £~ a»,n(l - a;)™Log 2 (l - x) , 

where ao,n is arbitrary (e.g. 1) 

ig 3 + 2r(r- l)(r-2) 



and m = n + ri 



(m-3)(m-4) iq 3 + 2(m - l)(m - 2)(m - 3) 

-0>n-2,n H / ^77 7T a n-l,n j (Vj./yJ 



m(m — 1) " 1 (m — 2)(m — l)m 



whereas &o,r 2 = ^ a o,ri and 6i r2 is arbitrary while 



h,r 2 = ^(mh,r 2 + 4a , n - llai,n) > (C.80) 
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, _ (n-2)(n-3) , ig 3 +2(n-2)(n-l)n , ll+3(n-4)n 

Wn,r 2 — n(n+l) °™-2,r 2 ~T~ n (n 2 -l) °n-l,r 2 n(n 2 -l) u ra-3,riT 
6(n-2)n+4 , j-3n _ 
^ n(n 2 -l) u n-2,n T „( n 2_i) "n-l,n • 

Therefore, Co, r3 = — i^^o,r 2 an d Ci i7 . 3 = 4 ( b °- r 2+ bl - r 2) + 6 °° 3 T ' 1 and C2 >r3 is arbitrary while 
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n(n-2) a ™-2. r l " 



(C.81) 



c 3 ,r 3 = g(ig3C2,rs + 2 &o,r 2 + 8&i, r2 - 22b 2 , r2 - 12a , n + 12ai >ri ) , (C.82) 

_ (n-3)(n-4) ig 3 +2(n-3)(n-2)(n-l) 2(3n(6-n)-26) , 

^n,r 3 — n(n-l) c n~2,r 3 "I n(n-l)(n-2) ^n-l,r 3 "T „( n _i)( n _ 2 ) w n-3,r 2 -^ 

, 4(ll-12m+3m 2 ) 7 _ 2(2-6m+3m 2 ) , 6(ra-3) . 12 , 

+ n(n-l)(ra-2) °«-2,r 2 n (ra-l)(n-2) n-l,r2 n(n-l)(n-2) a «-4,ri + „(„_i) a n-3,n T 



(C.83) 



C.2.8 Solution with qs = and h = 2 

In this case a , r is arbitrary (e.g. 1). For the first solution with r\ = 2 we have other coefficients 
On>o,ri = 0. The second solution, with r 2 = 1, has arbitrary a^ T2 (we can take it as 0) and 
a>2,r 2 — 0. The third one is with arbitrary ai )r3 ,a2, r3 (we also set them to 0). All these solutions 
look like 



Ui(x) 


= (1 


- x) 1 En=0 (X n,r 1 (1 


-x) n 


= (1- 


•^0 ^0,r 3 j 




u 2 (x) 


= (1 


— X) 2 En=0 a n,r 2 (1 




= (1- 


-^) ( ^0,r 2 5 


(C.84) 


u 3 (x) 


= (1 


— x) r3 En=0 a «>»* 3 (l 




= a ,r 3 







Gathering together all this solutions we have 

u(x) = A + B(-xf + C{x-l f , (C.85) 

where A, B, C are arbitrary. As we can see this solution corresponds to the solution with ^3 = 
and arbitrary h {0, 1}. For other integer h £ {0, 1} we can observe the same correspondence. 



C.3 Solutions for s = around x = oo 

C.3.1 Solutions for (/3 7^ and h Z 

When x = oo + we have an indicial equation 

(h + n + r)(h - 1 + r + n)(n + r) = (C.86) 

so its solutions are ri = 0, r 2 = 1 — h and r 3 = h. As we can see we have to cases when 
h £ 7L (one solution with Log(x)) and /i G Z (one solution with Log(x) and one solution with 
Log 2 (a;)). As we will see below we have to consider solutions with g 3 = 0, separately. In the 
first case 

U1 (X) = (1/^^=0^(1/*)", 

u 2 (x) = (l/xrE:=oV 2 (V^, (C87) 
«s(x) = (l/x)^Er=o^, 3 (l/^) n + (l/^) r2 Er=o«n,, 2 (l/a:)"Log(l/x), 
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where ao, r is arbitrary (e.g. 1) 



(iq 3 + r(h + r)(h + 2r)) 
(1 + r)(h + r)(h + 1 + r) 



Q-0,r 



(C.88) 



and m = n + r 

(m — 2)(m — 1)(2 — /i — m) 



(/i + m)(/i + m — l)m 



ig 3 + (h - 1 + m)(m -l)(h + 2(m - 1)) 

Q"n—2,r H 77 I \ / I ~ 1 \ ^n— l,r ) 



(/i + m)(/i + m — l)m 



(C.89) 



whereas &o,r 3 
is arbitrary 



Qo,r 2 ( one can notice that coefficient 5 ,r 3 is valid only for q 3 ^ 0) and bi :T3 



, iq 3 +2+(h-3)h h 

v 2,r 3 2(2-ft) li^a 2(2-ft) L " i : r 3 1 2(2-ft) " u ) r 3 



-3ft n , 6+fe(fe-6) 
r"l,r 3 "T ofo_M "0,-. 



(C.90) 



and 



, __ i( ?3 +(l+/t-m)(m-l)(/i-2(m-l)) , . 



(ft+l-rn)(ft-ro+2)(ro-2) 
r 3 (ft— m)(m— l)m "— 2,r 3 (ft— m)(m— l)m 

, ft 2 +ft(7-4m)+(m-2)(3m-4) __ fe 2 -6fe(m-l)+6(m-l) 2 

a n-3,r- 3 (ft-m)(m-l)m a n-2,r 3 + 



(/i— m)(m— l)m 
2m-3m 2 +fe(2m-l) 

(ft— m)(m— l)m n— l,r3 > 



(C.91) 



where m = r 3 + n. 



C.3.2 Solution with q 3 = 

In this case ao, r is arbitrary (e.g. 1). For the first solution with r\ = we have ai ri = and 
for the third one ai i7 . 3 is arbitrary (let us take 0) . Here it also turns out that we do not need 
Log solutions. The second solution is more complicated. One can derive exact formula for 

nk — h r(l — h + n) , 
; __ l ~k+l = T{l-h)T{n + 2) a ^ {Cm) 



and performing summations 



U X {x) = (l/x) ri J2n=0 a n,n( 1 / X ) n = a 0,r 1 , 

=0 r(i-ft)r(n+2) ' 



U 2 (x) = (1/xr EZo*n,r 2 (l/xr = «0,, 2 (1/*) ^ E~0 FW) W 



= -a ,, 2 i((a;-l) h -x ft ), 

Ua(x) = (l/x) r3 E"=0«n,, 3 (l/xr = ao ,r 3 (l/a;)- h - 



(C.93) 



Gathering together all this solutions we have 

u{x) = A + B{-x) h + C(x - 1)\ (C.94) 
where A, B, C are arbitrary. The above solution was presented by Lipatov and Vacca in Refs. 

[nsiini]- 
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C.3.3 Solution with q% ^ 0, q^ = and h = 1 



For integer h — 1, i.e. g 2 = 0, and g 3 ^ we have three solutions to the indicial equation (|C.86J) 
which are integer r 1 = 0, r 2 = 0, r 3 = — 1 so solutions are 



Ul (x) = (i/xy^Z *n, ri (i/xT 

u 2 (x) = (l/xy>j:Z b n , r2 (l/xr + (l/xyij:™ =0 a nin (l/x) n Log(l/x) 
«s(x) = (l/x) r3 Er=oCn,,3(l/^) n + 2(l/x)^Er=o^ 2 (l/x) m Log(l/x) + 
+ (V^) ri E"=o«n,n(l/^) n Log 2 (l/x) , 



where ao, r is arbitrary (e.g. 1) 



__ iq 3 + r(l + r)(l + 2r) 
ai ' r ~ (l + r )2(2 + r) G °' r 



and m = n + r 



(m - 2)(m - 1) ; 
(m + l)m 2 



On-2,r + 



2^3 + (m — l)m(2m — 1) 
(m + l)m 2 



whereas 



(C.95) 



(C.96) 



(C.97) 



(C.98) 



y n,r 2 



(n-l) 2 (n-2) , 

, in "n— 2,r 2 



n 2 (n+l) 
(n-l)(5-3n) 
^ n 2 (n+l) u "-2,ri 



, fa+2m 3 -3m 2 +m) i . l+6n(n-l) . 

n 2 (n+l) u n-l,r 2 ' n 2( n+1 ) «n-l,n~r 



(2+3n) 



n(n+l) a ™' r i 



and Co )r 



a o,ri C\r 3 is arbitrary while 



(C.99) 



c 2 



'■3 



«?3Ci,r 3 + &o,r 2 - 5&i, r2 + 3a , ri - 4a i 



,^1 > 



(C.100) 



-n,r 3 



(n-2) 2 (n-3) 

n(n-l) 2 c "-2,r 3 I - 
;i3+6(i 
n(n- 
2(3n-2) 



-6+n(13+n(2n-9))-iq 3 
n(n— l) 2 



Cn-1 



■'.:-! 



2(13+6(n-3)n) , _ 2(3n-l) , , 2(7-3n) _ ^ , 

n(n-l) 2 u n-l,r 2 n (n-l) l V'"2 n(n-l) 2 u n-3,n „( n -l) 2 u ri-2,n T 



2(n-2)(3ra-8) 7 . 
n(n-l) 2 °n-2,r 2 + 
6(3-2n) 



(C.101) 



n(n-l) u "-l>n ' 

C.3.4 Solution with q^ = q% = and /i = 1 

In this case a ,r is arbitrary (e.g. 1). For the first solution with r\ = we have a^ n = and 
for the third one ai >ra with r 3 = is arbitrary (let us take 0). The second solution is more 
complicated. We need Log— solutions 



1*2(2;) = (1/xr ^6 n , r ,(l/x) n + (l/x) ri ^a n , ri (l/x) n Log(l/x) . 



n=0 



n=0 



Using recurrence relations with b\ yT2 = ^ao, ri (and a n ri = for n > 0) 



&2.r. 



2,r 2 



y n,r 2 



3 U 1|)"2 5 

(n-l) 2 (n-2) 



, . («-l)«(2n-l) , 



J (n+1) u n-l,r 2 



(C.102) 



(C.103) 
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one can derive an exact formula for 6 n r2 = b\^ 2 ^ n ^ n and performing summations with arbitrary 
fro,r 2 (= 0) we have ^^L &„ ir2 (l/x) n = 2& lr2 (1 — Log(l — 1/x) + xLog(l — 1/x)), so that 

Ui{x) = (l/x) ri Yln=Q a n,r 1 ( 1 / x ) n = °0,ri , 

u 2 (x) = a , ri (l-Log(l-l/x) + xLog(l-l/x) + Log(l/x)), (C.104) 

Gathering together all this solutions we have 

= A + B(-x) + C((x - l)Log(x - 1) + (-x)Log(-x)) , (C.105) 
where A, B, C are arbitrary. 



C.3.5 Solution with 7^ 0, q^ = and h = 



For /i = 0, i.e. g 2 = 0, and ^3 7^ we have three solutions of the indicial equation (IC.86J) which 
are integer r\ = 1, r 2 = 0, r 3 = so solutions are 



ui{x) = (l/x) ri ^ =0 a„, ri (l/x) n , 

u 2 (x) = (l/xY^ZoK,r 2 (l/x) n + (l/x)^j:^ =0 a njri (l/x) n Log(l/x), 
u 3 (x) = (l/x)^^ =0 c n , r3 (l/x)™ + 2(l/x)^Er=o^ 2 (l/x)™Log(l/a;) + 
+ (V^) ri E" = o«n,n(l/x)"Log 2 (l/x) , 



where ao, r is arbitrary (e.g. 1) 



and m = n + r 



iq 3 + 2r 2 
r(l + r) 



r a o,- 



(m - 2) 2 ig 3 + 2(m - l) 3 

Q'n—2,r H 7 7\ o *^ra— l,r j 



(m — l)m 2 



whereas &o,r 2 = -^Oon an d b± r2 is arbitrary while 



(C.106) 



(C.107) 
(C.108) 



h,r 2 = t((2 + iq^)bi,r 2 + 6a , ri - 8a 1)ri ) , 



(C.109) 



Furthermore, Cq, 



» - 3 



(n-2)f_, ig 3 +2(n-l) 3 , 

n 2 W n _ 2 ,r 2 T n 2f n _i1 fn-1,7-2 



ri 2 (n— 1) 
I 6(n-l) 2-3n 
"I" n 2 «n-2,n "I" n (ri.-l) "n-l,n 



(n-2)(3n-4) 

n 2 (n-l) a «-3,ri 



+ 



— (2ao, n + &i,ra) and arbitrary and 



C 2 . 



1 5 

7(2 + iq3)ci,r 3 + 3& i,r 2 - 46 2 ,r 2 + 3a , ri - -a lin , 



(C.110) 



(C.lll) 



n 2 (n— 1) 



■^3+2(ra-l) 3 

n 2 (n-l) u n-2,r 2 - 

, _ 2(5-3») 12 

n(n-l) U «' r 2 n 2 (n-l) Un - 3 . r i ^ n^n-i , y 



(n-2) 2 ^ 3 +2(n-l) 3 _ 2 (n-2)(3n-4) , , 

„2 u n—2,r3 ~T n 2 (n—l) u n—l,rs ,.->,., > "n—2.r-> 



■ 12(n— 12 7 2(2-371) ; u: 

1 n2 "n— l,r 2 T "n,r 2 "T „2f„_i-\ 3,ri "T „2"n- 2,ri"r 

1 2(l-3n) 



(C.112) 
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C.3.6 Solution with q^ = qs = and h = 

In this case a r is arbitrary (e.g. 1). For the first solution with r\ = 1 we have recurrence 
relations with ax,n — | a o,ri 



(C.113) 

„ hU—U—n _|_ £Ll_ n 

u n,ri n+1 n— 2,fi ~r n _(_i u n- l,ri • 



(n-1) 2 „ i 2n : 



Summing series up we derive an exact formula for 

OO OO 

Ul (x) = (l/xrJ2 a n,rM^) n = 0-/x) ^(l/x) n a ,n— Y = ~Log(l - l/x) . (C.114) 



n=0 n=0 



In the second solution, with r 2 = 0, we have arbitrary a lr2 (let us take 0) so the second 
solution is v,2(x) = ao, r2 - The third one is with Log solutions 

OO OO 

u 3 (x) = (l/x) rs ^6 n , r3 (l/x) n + (l/x) r * ^a n , r2 (l/x)"Log(l/a;) . (C.115) 

n=0 n=0 

Using recurrence relations with arbitrary &o,r 3 

7, ("-2) 3 ft , 2(»-l)» , 

u n,r 3 n 2 "n— 2,r3 T n 2 "n— l,ra 

one can derive an exact formula for b n>rs = bi >r3 - and performing summations with arbitrary 
&o,r 3 (= 0) we have ^2^ =0 b nt r 3 (l/ x) n = — &o,r 3 Log(l — (l/x)). All these solutions look like 

111 (Z) = (l/x) ri Y^n=0 a n,n{l/x) n = -«l,r 3 Log(l - l/x) , 
M 2 (x) = (l/x)*" 2 Y2n=0 a n,r 2 xn = a 0,r 2 > 

u 3 (x) = (l/x) r3 Xir=o & «^ 3 (l/x) n + a ,r 2 Log(l/x) = -6 ,r 3 Log(l - l/x) + a ,r 2 Log(l/x) . 

(C.117) 

Gathering together this all solutions we have 

u(x) = A + BLog{-x) + CLog(x - 1) , (C.118) 
where A, B, C are arbitrary. 

C.3.7 Solution with q 3 ^ and h = 2 

For ft. = 2 and g 3 ^ we have three solutions to the indicial equation IJC.86J) which are integer 
7*1 = 0, r 2 = —1, r 3 = —2 so solutions are 

Ul (x) = (l/^'^oV.iV^, 

u 2 (x) = (l/xrEr=o^ 2 (V^ n + (l/^) ri E~oVn(l/^) n Log(l/x), 

«S(X) = (l/xPEZo^,r 3 (l/xr + 2(l/xrEZoKr 2 (y^Og(l/x) + 

+ (V^) ri E~o«n,n^Log 2 (l/^), 
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where ao, n is arbitrary (e.g. 1) 



and m = n + r\ 



" ( r + l)(r + 2)(r + 3) (CU20) 



(m-l)(m-2) ig 3 + 2m(m 2 - 1) rnon 

On-2,ri + 7 rTTT7 , 7~\ a n-l,n , {Kj.LZL) 



(m + l)(m + 2) (m + 2)(m+l)m 

whereas &o,r 2 — ijr^o.ri and 6i >r2 is arbitrary while 



&2,r 2 = ^{mh tr2 + 4a 0jri - llai, ri ) , (C.122) 



, (n-2)(n-3) , , zg 3 +2(f 1 -2)(n-l)ra , ll+3(ra-4)n , 

u n,r 2 n(n+l) u "-2,r 2 1" n (n 2 -l) w n-l,r 2 n(n 2 -l) u n-3,riT 

I 6(n-2)n+4 , l-3n - 

" r n(n 2 -l) u «-2,n -r n ^ n2 _^ u n _i in . 
Moreover, Co, r3 = — 7^r^o,r 2 and ci jr . 3 = 4< ~ b °''"^ fel ' r2 ^ + 6 "° 3 n and C2, r3 is arbitrary while 



(C.123) 



c 3 ,r 3 = -^{iq^c 2 ,r 3 + 26 ,r 2 + 8b hr2 - 22& 2jr2 - 12a ,n + 12ai,n) > (C.124) 



_ (n-3)(n-4) , ig 3 +2(n-3)(n-2)(n-l) , 2(3n(6-n)-26) , , 

^n,r 3 — n(n-l) c «-2,r 3 1 n(n-l)(n-2) L n-l,r 3 I - „(„_i)( n _ 2 ) w n-3,r 2 "^ 

, 4(ll-12m+3m 2 ) 7 _ 2(2-6m+3m 2 ) i 6(ra-3) . 12 , 

^ n(n-l)(n-2) °"-2,r 2 n (n-l)(n-2) °n-l,r 2 n{n-l){n-2) a «-4,ri + n ( n _i) a n-3,ri + 
6 



n(n-2) °™-2,- 



C.3.8 Solution with ^3 = and h = 2 



(C.125) 



In this case ao, r is arbitrary (e.g. 1). For the first solution with r\ = we have other coefficients 
an>o,ri = 0. The second solution, with r 2 = —1, has arbitrary ai ;r2 (we can take it as 0) and 
«2,r 2 = 0. and third one with arbitrary ai r3 , a 2jrs (we also set them to 0). All these solutions 
look like 

Ul(x) = (l/x) ri Ysn=O a n,rA l l X ) n = «0,r 3 , 

u 2 (x) = (l/x) r *j:Z an,r 2 (l/xT = xao,r 2 , (C.126) 

Gathering together all this solutions we have 

u(x) = A + B{-xf + C{x-l f , (C.127) 

where A, B, C are arbitrary. As we can see this solution corresponds to the solution with q 3 = 
and arbitrary h {0, 1}. For other integer h £ {0, 1} we can observe the same correspondence. 
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Appendix D 



Coefficients in the eigenequations for 

7V = 4 



Coefficients for the eigenequations of q 3 and g 4 for N = 4 presented in (|5.55j) and (|5.56|) look 
like 

*o,o = -^liK-h + si + s 2 + s 3 + s 3 )((s 2 - s 4 )(l-/i + si + s 2 - S 3 + S 4 ) + 

+ (s 2 2 - s 2 (-l + /i + si - s 3 + 2s 4 - 2s 4 xi)+ (D.l) 
+s 4 (l + si - s 3 + s 4 - 2(1 + s 4 )xi + h(—l + 2xi)))x 2 )) - q-i , 

ti t o = —i((si + s 2 )(l - h + si + s 2 - s 3 + s 4 ) - (2 + h - h 2 + si(3 + si) - s 2 + 

+4/is 2 - 3s 2 2 - 3s 3 - 2s 2 s 3 + S3 2 + s 4 + s 4 2 - 2si(s 2 + s 3 + s 4 (xi - 1))+ (D.2) 
-2s 2 s 4 xi)x 2 - (s 2 + s 3 - 2)(1 + h + si - s 2 - s 3 + s 4 - 2s 4 a;i)x 2 2 ) , 

*Q,1 = (rra-l)^ (Kih ~ Si ~ S 2 + S3 - S 4 )(si + S 4 + (s 1 + S 2 - 2)xi) + 

+ ((si + s 4 )(l - 2h + 2si - 2s 3 + 2s 4 ) - 2(1 + h 2 + s 2 + s 2 2 - s 3 + s 2 s 3 + 
-S3 2 + (4 + s 2 )s 4 + 2s 4 + si(2 + s 2 + s 3 + 2s 4 ) - h(2 + s : + 2s 2 + 3s 4 ))xi + 
-(si + s 2 - 2)(1 + 2s 4 )xi 2 )x 2 + (-(si + s 4 )(l - h + si - s 2 - s 3 + s 4 ) + 
+ (2 + /i 2 + si 2 + s 2 + s 3 - 2(s 2 + S3) 2 + h(-3 - 2si + s 2 + s 3 - 6s 4 ) + 
+7s 4 + (s 3 - s 2 )s 4 + 5s 4 2 + si(3 - s 2 + s 3 + 6s 4 ))a;i + 
+ (-2 - h 2 - 2s 2 + s 2 2 - 3s 3 + s 2 s 3 + si(-2 + s 2 - 3s 4 ) + 
-(7 + 2s 2 + 5s 3 )s 4 - 5s 4 2 + h{3 + si + s 3 + 6s 4 ))xi 2 )x 2 2 )) , 

*2,o = - X2 (x 2 -i) (i( x i - l)xi(h - si - s 2 + s 3 - S4 + (2 - 2h + 3si - 2s 3 + 3s 4 + 

+ (-4 + si + s 2 - s 4 )xi)a; 2 + (-2 + h- 2 Sl + s 2 + s 3 - 2s 4 - 2hxt+ (DA) 
+ (4 + si + s 2 + 2s 3 + 4s 4 )xi)x 2 2 )) , 



(D.3) 



to,2 = -i(x 2 - l)x 2 (h - 2 - 2si - 2s 2 + s 3 - s 4 + 
+ (4 + h + si - 2s 2 - 2s 3 + s 4 - 2s 4 xi)x 2 ) , 



(D.5) 



h,i = -i(x 1 (2/i-3s 2 + (2/i-4(s 2 + s 3 ) + (s 2 -2)xi)x 2 + 

+ (s 2 + s 3 - 2)(xi - l)x 2 2 + 2(s 3 + x 2 )) + si(x 2 - 1 + X!{x!X 2 - 3))+ (D.6) 
+s 4 (x 2 - 1 - 2xi(l + 2x l x 2 ))) , 
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^3,0 = - l) 2 xx 2 , 



(D.7) 



*o,3 = -i(x 2 - l) 2 x 2 2 , 



(D.8) 



t 21 = i(xi - l)x l (x 2 - 1 + 2x 1 x 2 ) , 



(D.9) 



ixi(x 2 - l)x 2 ((xi - l)x 2 - 2) 



(D.10) 



and where 



/i,o = ( X2 -\y2 X2 {{-h + si + s 2 - s 3 + s 4 )((s 2 - l)(si + s 4 ) - (-2 + si + 
+s 2 )(l + s 4 )a:i) + (-(si + s 4 )(-h(-2 + s 2 ) + si(s 2 - 2) + 3s 3 + 
-2s 4 + s 2 (-2 + s 2 - 3s 3 + s 4 )) + (-5s 2 + 3si 2 s 2 + s 2 2 + s 3 3 + 5s 3 + 
-2s 2 s 3 + s 3 2 - s 2 s 3 2 + h 2 {-l + s 2 - s 4 ) + (-5 + 4s 2 2 + s 2 (5 + s 3 ) + 
-s 2 (7 + 3s 3 ))s 4 + 2(s 2 - 4)s 4 2 - s 4 3 + h{5 + s 1 - 4sxs 2 - 2s 2 2 + 
+ (9 + s 1 - 3s 2 )s 4 + 2s 4 2 ) + si(-5 + s 2 + 4s 2 2 + 6s 2 s 4 - 3s 3 (l + s 4 ) + 
-s 4 (8 + s 4 )))xi - (-2 + si + s 2 )(-/i + si + s 2 + s 4 )(l + 2s 4 )a;i 2 )a: 2 + 
+(— (si + s 4 )(si + s 2 - /i(l + s 2 ) - 3s 3 + s 4 + s 2 (si + s 2 + 3s 3 + s 4 )) + 
+ (3s 2 - 5s 3 + 2(s 2 2 + s 2 (l - s 2 )s 3 + (-1 + s 2 )s 3 2 ) + 3s 4 + 2(s 2 2 + 
+4s 2 (l + s 3 ) - s 3 (4 + s 3 ))s 4 + (6 + 3s 2 - s 3 )s 4 2 + s 4 3 + /i 2 (l + s 4 ) + 
+ Sl 2 (l - 2s 2 + s 4 ) + si(3 - 2s 2 2 + s 3 + (7 + s 3 )s 4 + 2s 4 2 + s 2 (3+ 
+4s 3 + s 4 )) + h(-3 + s 3 + 2si(-l + s 2 - s 4 ) + (-7 + s 3 )s 4 - 2s 4 2 + 
-s 2 (3 + 2s 3 + 3s 4 )))xi + (2 Sl 2 s 2 + s 3 + s 2 (-3 + s 2 2 + (-2 + s 2 )s 3 ) + 
+/i 2 (-l + s 2 - 2s 4 ) - 3s 4 + (s 3 - s 2 (7 + 5s 3 ))s 4 - (7 + 3s 2 + 2s 3 )s 4 2 + 
-2s 4 3 + h(3 + si + s 2 - 3sis 2 - 2s 2 2 + s 3 - s 2 s 3 + 2(4 + si + s 2 + s 3 )s 4 + 
+4s 4 2 ) + si(3s 2 2 + s 2 (2s 3 + s 4 ) - (3 + 2s 3 + s 4 )(l + 2s 4 )))a; 1 2 )a; 2 2 + 
+ ((si + s 4 )(-/i + si + s 2 + s 3 + s 4 )) - (/i 2 s 2 + 2sis 2 + si 2 s 2 + 2s 2 2 + 
+s 2 3 - 2s 3 - sis 3 + s 2 s 3 + 4sis 2 s 3 + 2s 2 2 s 3 - s 3 2 + s 2 s 3 2 + 
+ (5s 2 2 - s 3 (5 + si + s 3 ) + s 2 (2 + 5si + 6s 3 ))s 4 + (4s 2 - s 3 )s 4 2 + 

+ /i(s 3 - 2s 2 (l + Si + S 2 + S 3 ) + (— 5S2 + S3)S4))^1 + 

+ (/i 2 s 2 + 2s 2 2 + s 2 3 - 2s 3 + 2s 2 s 3 + 

+3s 2 2 s 3 - s 3 2 + 2s 2 s 3 2 + (2s 2 (l + 2s 2 ) + 5(-l + s 2 )s 3 - 2s 3 2 )s 4 + 

+ (3s 2 - 2s 3 )s 4 2 + sis 2 (2 + s 2 + 2s 3 + 3s 4 ) + 

-h(-s 3 + s 2 (2 + si + 2s 2 + 3s 3 ) + 4s 2 s 4 - 2s 3 s 4 ))a;i 2 )x2 3 ) , 
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/o,o = ( X2 -i)2 (s 2 s 4 {h 2 - h + si- Ahsi + 3si 2 + s 2 - 2hs 2 + 4sis 2 + s 2 2 + s 3 + 

— S3 2 + S4 — 2/iS 4 + 4siS4 + 2s 2 S4 + S4 2 + 2(/lSi — Si 2 — SiS 2 — /iS 3 + 

+2sis 3 + s 2 s 3 + S3 2 - s 1 s 4 + s 3 s 4 + (1 - h + 2si + s 2 + s 4 ) 

(-/i + si + s 2 + s 3 + s 4 )xi)x 2 + x(-si 2 - si(l + 2s 2 + 4s 3 + 2s 4 ) + (D.12) 

+2si(l + s 2 + 2s 3 + s 4 )xi + h 2 (2xi - 1) + 

+ («2 + S3 + s 4)(-l - s 2 - s 3 - s 4 + 2(1 + s 2 + 2s 3 + s 4 )xi) + 

+ 2si + 2s 2 + 2s 3 + 2s 4 - 2(1 + Sl + 2s 2 + 3s 3 + 2s 4 )a;i))a; 2 2 )) - g 4 , 

/o,i = -^(s 4 (si 2 (l + 2x ± x 2 - x 2 2 ) + si{l-h + 2s 2 - s 3 + s 4 + (-2-/1+ 
+4s 2 + 3s 3 + s 4 + 2(1 - h + 2s 2 + s 4 )a;i)a; 2 - (1 + 5s 2 + s 3 + 2s 4 + 
+2/i(xi - 1) - 2(-2 + 4s 2 + 2s 3 + s 4 )a;i)a; 2 2 - (-2 + s 2 + s 3 )x 2 3 ) + 
+s 2 2 (l + (2xi - l)x 2 )(l + x 2 (3 + x 2 )) + (1 + h - s 3 - s 4 )x 2 
x(h + h(2xi - l)x 2 - (2 + s 4 - 2x 2 ){l - x 2 + 2x ± x 2 )+ (D.13) 
+s 3 (l + x 2 (-2 + x 2 - 2x1X2)))+ 

+s 2 ((l + (-1 + 2xi)x 2 )(l + s 4 + x 2 + As 4 x 2 + (-3 + s 4 )x 2 2 ) + 
-h{\ + (2xi - l)x 2 (l + x 2 (4 + x 2 )) + 
+s 3 (-l + x 2 (3 - 2x 2 2 + 4^(2 + x 2 ))))))) , 

/2,o = (^ZTp^ - si - s 2 + s 3 - s 4 )(2 + si - s 2 + s 4 + 

+ (-4 + si + s 2 - s 4 )xi) + (-Ah + 4si - 2/isi + 2si 2 + 4s 2 + Sis 2 - s 2 2 + 
— 6s 3 — 3siS 3 + 3s 2 s 3 + 4s 4 — 2hs 4 + 4siS 4 + s 2 s 4 — 3s 3 s 4 + 2s 4 2 + 
-(h 2 - sis 2 - s 2 2 + 3sis 3 + 3s 2 s 3 - S3 2 + 5sis 4 + 4s 2 s 4 - 3s 3 s 4 + 5s 4 2 + 
+ ll(si + s 2 - s 3 + s 4 ) - h(ll + Si + 6s 4 ))xi + 

-(-h + si + s 2 + s 4 )(s 1 + s 2 - 2(2 + s 4 ))xi 2 )x 2 - (s x 2 + 2siS 2 + s 2 2 + 

-3sis 3 + 3s 2 s 3 + 2s x s 4 + 2s 2 s 4 - 3s 3 s 4 + s 2 - h(2 + s 1 + s 2 + s 4 )+ 

+2(si + s 2 - 3s 3 + s 4 ) + h(7 + 2si + 3s 2 - s 3 + 5s 4 )xi + 

-(si 2 + 2s 2 2 - s 3 (13 + 2s 3 ) + 7s 4 - 6s 3 s 4 + 4s 4 2 + 

+si(7 + 3s 2 + s 3 + 5s 4 ) + s 2 (7 + 4s 2 + 6s 4 ))xi + h 2 {x x - l)x : + 

-h(5 + si + 2s 2 + s 3 + 4s 4 )xi 2 + (s 2 2 - 3s 3 + 5s 4 - s 3 s 4 + 3s 4 2 + 

+si(5 + s 2 + 2s 3 + 3s 4 ) + s 2 (5 + 3s 3 + 4s 4 ))xi 2 )x 2 2 + 

—(xi - l)(—hs 2 + sis 2 + s 2 2 - 2s 3 - S1S3 + s 2 s 3 + s 2 s 4 - s 3 s 4 + 

-(sis 2 + s 2 2 - 4s 3 + s 2 s 3 - s 3 2 + h(s 3 - s 2 ) + s 2 s 4 - 3s 3 s 4 )xi)x 2 3 )) , 

/o,2 = s 4 x 2 (2 - h + 2si + 2s 2 - s 3 + s 4 + (-6 - s x + 3s 3 + 

+2(-/i + 2(1 + si + s 2 ) + s 4 )xi)x 2 + (D.15) 
-(si + /i(2xi - 1) - (2(-2 + s 2 + s 3 ) + s 4 )(-l + 2x0)x 2 2 ) , 
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fi,i = -^ziihsi - si 2 - 2s\s 2 + S1S3 + hs4 - 2sis 4 - 2s 2 s 4 + s 3 s 4 - s 4 2 - 2hxx + 
+2s 4 Xi + hsix 1 — s\ 2 Xi + 2s 2 Xi + hs 2 xi — 2siS 2 ^i — s 2 2 x 1 — 2s 3 x 4 + 
+s 1 s 3 x 1 + 2s 4 X! — 2/is 4 xi + 2s 4 s 4 a; 1 + 2s 2 S4^i — 2s 3 s 4 x 4 + 2s 4 2 x 1 + 
-(-(si + s 4 )(2 -2h + 2si + 3s 2 + 3s 3 + 2s 4 ) + (h 2 + s 1 + 3s 2 + 
+As 1 s 2 + 2s 2 2 - 3s 3 + 3sis 3 + 3s 2 s 3 - s 3 2 + (5 + 6s 4 + 10s 2 - 6s 3 )s 4 + 
+5s 4 2 - h(3 + si + 3s 2 + 6s 4 ))^i + ((-2 + s 4 + s 2 )(-/i + s 1 + s 2 ) + 
+ (-2 + Ah - 5«i - 5s 2 )s 4 - 4s 4 2 )xi 2 ):r 2 + (-(si + s 4 )(4 - ft + s x + 
-3s 3 + s 4 ) + (ft 2 + si + si 2 + 3s 2 + 5sis 2 + 2s 2 2 + s 3 + s 4 s 3 - 2s 3 2 + 
+/i(-3 - 2si - 3s 2 + s 3 - 4s 4 ) + lls 4 + 4(si + s 2 - 2s 3 )s 4 + 3s 4 2 )xi+ 
+ (-h 2 + si - s 2 - 4sis 2 - 3s 2 2 + s 3 - 2s 4 s 3 - 4s 2 s 3 - 5s 4 + 
-3sis 4 - 4s 2 s 4 + 5s 3 s 4 - s 2 + h(l + si + 4s 2 + s 3 + 2s 4 ))xi 2 )a; 2 2 + 
+ (-2 + s 2 + s 3 )(xi - l)(si + s 4 - (-/i + s 2 + s 3 + 3s 4 )xi)x 2 3 ) , 

/3,0 = X2 (i- X2 ) x l 2 ( x l ~ l ){~ h + S l + S 2 ~ S3 + S 4 - (Si + S 2 - 2s 3 + S 4 + 

- 1) - (si + s 2 + s 4 )x l )x 2 + s 3 (x 1 - l)x 2 2 ) , 



(D.16) 



(D.17) 



^0,3 



-s 4 (x 2 - l)x 2 2 (l + (2xi - l)x 2 ) , (D.18) 



/2,i = - 2s 2 + s 3 - 2xi + (3si + 3s 2 - 2s 3 + s 4 )xi + 

+ (s 3 (6:ri — 3) + Xi(3 — As 1 + s 4 — 2xi + 3«iXi)+ 

+s 2 (3 + xi(-5 + 3xi)))x 2 + ( Xl - 1)(3 + 3s 3 (x! - 1) - Sl x x + 2s 4 xi)x 2 2 + 
-(-2 + s 2 + s 3 )(xi - 1)W - ft(x 2 + 2xi - 1)(1 + (xi - l)x 2 ))) , 



(D.19) 



/i,2 = x 2 (-(h - 2s 2 + s 3 )xi + si(l + 2x x - x 2 )(l + (xi - l)x 2 )+ 

+ Xl x 2 (-3 + 3s 3 - hx x + 2s 2 xi - (3 + h - 2s 2 - 2s 3 )(x 1 - l)x 2 )+ (D.20) 
+54(0;! - l)(x 2 - 1)(1 + (-1 + 3xi)x 2 )) , 

/ 3 ,i = (xi - l)xi 2 (l + (xi - l)x 2 ) , (D.21) 

A,3 = -xi(x 2 -l)x 2 2 (l + (xi-l)x 2 ), (D.22) 

h,i = a:i(l + Xi(x 2 -2)-x 2 )x 2 (l + (xi-l)x 2 ), (D.23) 

/ 4 ,o = /o,4 = 0. (D.24) 

If we take homogeneous spins s = S\ = s 2 = . . . our equations simplify and the coefficients 
for q 3 look like 

*o,o = l)(h-As)s(x 1 - l)x 2 -q 3 , 
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*i,o = ^zr ) (i(2(h-2s)(s+(s-l)x 1 )-2(-s(l-2h + 2s) + 
+ (-1 +h- 4s)(-l + h- 2s)x! + (s - 1)(1 + 2s)x 1 2 )x 2 + 
-((h -2)(h- l)(xi - 1)^ + 4s 2 a;i(3xi - 1) - 2s(h - 1+ 
-3(h - 2)x 1 + (Ah - 7)x 1 2 ))x 2 2 )) , 

t ,i = -i((l + h)x 2 (-2 + h + 2x 2 )+4:S 2 (l + x 2 )(l+x 1 x 2 )-2s(-l + h+ 
+2hx 2 + (1 + h + 2x l )x 2 2 )) , 



and for g 4 



t2,o = ^=iy (gi ~ + 2s + 2(-l + - 2s + 2x ± )x 2 + 

+ (2 - h + 2s + 2(h - 2 - 4s):r 1 ):r 2 2 ) , 



= —2i(xi(l + x 2 )(h + x 2 — xix 2 ) + s(—l + x 2 + 
+xi 2 (x 2 - l)x 2 — x±(l + x 2 )(3 + x 2 ))) , 



(D.26) 



(D.27) 



(D.28) 

fo,2 = -i(x 2 - l)x 2 (-2 + h - 4s + (h - 2(s - 2 + sxi))x 2 ) , (D.29) 



(D.30) 

3,0 = ~i(x! - 1)V 2 , (D.31) 



to,s = -i{x 2 -l) 2 x 2 2 ) (D.32) 
h,\ = - l)x 1 (2xix 2 + x 2 - 1) , (D.33) 
ti, 2 = ixi(x 2 - l)x 2 (x 2 ( Xl - 1) - 2) (D.34) 



/o,o = (^(s 2 (/i 2 (l + z 2 )(l + (2zi-l):r 2 )+ 

-h{\ + 8s) (1 + x 2 ){\ + (2X! - l)x 2 ) + 2s(2(l + x 2 )(l - x 2 + 2 Xl x 2 )+ (D.35) 
+s(7 + x 2 (2 - 9x 2 + 16xi(a:2 + 1)))))) - g 4 , 

/i,o = ( X2 -i)2 X2 {h 2 xix 2 {-\ + (ari - 1)(-1 + s(x 2 - l))x 2 ) + h(x 1 (3x 2 - 2) 
x(l + (xi - l)ar 2 ) - 2s 2 (xx - l)(ar 2 - 1)(1 + 2:r 1 :r 2 + (4xi - l)a; 2 2 ) + 
+s(2 + x 2 (2(x 2 - 2) + Xl (x 2 - I0)(x 2 - 1) - x x 2 {2 + (x 2 - ll)x 2 )))) + 
+2s(2s 2 (x 1 - l)(x 2 - 1)(1 + x 2 + 3x lX2 + (5xi - 2)x 2 2 ) + 
+s(-2 + x 2 (3 + 10 Xl (x 2 - 1) - x 2 2 + m 2 (3 + {x 2 - 12)x 2 )))+ 
+xi(2 + x 2 (-5 + x 2 (2 + x 2 ) - xx(-3 + x 2 (4 + x 2 )))))) , 



(D.36) 



/o,i = ^(s(-(l + h)x 2 (h-2 + 2x 2 )(l + (2x 1 -l)x 2 )-As 2 (l + x 2 (A + x 2 

x (3 + 2x 2 ) + xi(3 + x 2 ){\ + 3x 2 ))) + 2s(-l + h(l + (2 Xl - l)x 2 )(l+ (D.37) 
+x 2 (4 + x 2 )) + ar 2 (-l + (7 - 5x 2 ):r 2 + 2x^-1 + x 2 (Ax 2 - 1)))))) , 
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h,o = { X2 -i)*x 2 (^i(-4g - h 2 xix 2 {l + (x! - l)x 2 ) + 2s(3x 2 - x 2 3 + 
-2x 1 2 x 2 (x 2 (3 + x 2 )) + xi(x 2 - l))(-4 + x 2 (7 + 3x 2 )) + 
-s{-{x 2 - l) 3 + Sa^x, 2 + Xl {l + x 2 {7 + (x 2 - 9)x 2 ))))+ (D.38) 

+ (xi - l)x 2 )(2 - 2x 2 + xi(5x2 - 4)) + s(-(x 2 - 1) 3 + 
+8x! 2 x 2 2 + x 1 (l+ x 2 (7 + (x 2 - 9)x 2 )))))) , 

= sx 2 (-h(l + x 2 )(l + (2x l - l)x 2 ) + 2(1 - 3x 2 + 2x 2 (x 1 + x 2 - 2x 1 x 2 ) + 

( D.39 ) 

+s(2 + x 2 (l - 3x 2 + 5x^1 + x 2 ))))) , 

= ^ L T j(-/i 2 x 1 x 2 (l + (xi - l)x 2 ) + 2s(-3s + 2xi + sxi + (2 + 3(xi - l)xi+ 
+s(4 + xi(5xi - 13)))x 2 - (2(2 + (a* - 4)xi) + s(-l + xi(8xi - 5)))x 2 2 + 
-(s - l)(xi - l)(5xi - 2)x 2 3 ) + h(- Xl (l + (xi - l)x 2 )(2 + x 2 (2x 2 - 1)) + 
+2s(l + x 2 (-2 + x 2 + xx(5 - xi + 4(xi - l)x 2 + (x : - l)x 2 2 ))))) , 

= ^3Tj(a;i-l)a;i 2 (/i + /i(xi-l)x 2 + s(-2 + x 2 (l + x 2 -x 1 (3 + x 2 )))), (D.41) 

/o,s = s(x 2 - l)x 2 2 (l + (2x! - l)x 2 ) , (D.42) 

/2,i = -^pr( x i( l - s - + 5sxi + xi(3 - 2xi + s(6xi - 2))x 2 + 

+(x! - 1)(3 + s(4x! - 3))x 2 2 - 2(s - l)(x! - l) 2 x 2 3 + (D.43) 
-/i(x 2 + 2x 1 -l)(l + (x 1 -l)x 2 ))), 



fi,2 = x 2 (-x 1 (/i+(3 + /i)x 2 )(l + (x 1 -l)x 2 ) + s(2 + 2x 1 -4x 2 + 
Xi(7 + x x )x 2 + (X! - l)(7x! - 2)x 2 2 )) , 



(D.44) 

/ 3 ,i = (X! - 1)X! 2 (1 + (X! - l)x 2 ) , (D.45) 

A,3 = -x 1 (x 2 -l)x 2 2 (l + (x 1 -l)x 2 ), (D.46) 
h,2 = x 1 (l + x 1 (x 2 -2)-x 2 )x 2 (l + (x 1 -l)x 2 ), (D.47) 
Ao = /o,4 = 0. (D.48) 
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